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Abstract. This review is devoted to the analysis of the mutual consistency of the spin 
foam and canonical loop quantizations in three and four spacetime dimensions. In the 
i-H three-dimensional context, where the two approaches are in good agreement, we show how 

the canonical quantization a la Witten of Riemannian gravity with a positive cosmologi- 
q cal constant is related to the Turaev-Viro spin foam model, and how the Ponzano-Regge 

amplitudes are related to the physical scalar product of Riemannian loop quantum gravity 
( without cosmological constant. In the four-dimensional case, we recall a Lorentz-covariant 

W) formulation of loop quantum gravity using projected spin networks, compare it with the new 

spin foam models, and identify interesting relations and their pitfalls. Finally, we discuss 
ff) the properties which a spin foam model is expected to possess in order to be consistent with 

^ the canonical quantization, and suggest a new model illustrating these results. 
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1 Introduction 

Canonical loop quantization and spin foam models are two complementary approaches to quan- 
tum gravity [16, 21, 126, 136, 142]. The former is designed to describe a quantum space in the 
frozen time formalism, via a (kinematical) Hilbert space on which certain constraints encoding 
the dynamics of the theory are imposed. In contrast, the latter is supposed to be a version of 
the discretized path integral for gravity, or a sum over histories where every history represents 
a quantum spacetime. Each of the approaches has its own advantages and shortcomings. For ex- 
ample, the canonical approach is well grounded and relies on a solid mathematical construction, 
but at the same time in this framework all the dynamical questions look extremely complicated. 
On the other hand, spin foam models bring powerful geometric methods and suggest a simple 
way to recover the dynamics of the theory, but their derivations are much less robust than the 
corresponding canonical constructions. 

In this situation, it is of primary importance to consider spin foam and canonical quantizations 
in parallel, in order to be able to learn some lessons in one approach and use them in the second 
one. If such an exchange of ideas and results is successful, it would allow to overcome the 
problems arising in each particular quantization scheme and, if we are lucky, to find a model 
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for quantum gravity resulting from both approaches. This implies that the results derived from 
a spin foam model, which is claimed to be a possible candidate for quantum gravity, must be 
in agreement with the canonical quantization. If the agreement is found, this can be considered 
as a serious argument in favor of the model. If not, then either the spin foam model or the 
canonical quantization (if not both) should be reconsidered. 

Thus, the study of mutual consistency of the loop and spin foam quantizations is crucial for 
getting new insights and eventually a better understanding of the existing constructions. This 
issue will be the main subject of this review, as our goal is to present the current understanding 
of the relations between the results of these two approaches. In particular, we would like to show 
an amazing interplay between the canonical and covariant quantizations of three-dimensional 
gravity, possibly coupled to point-like particles, and to analyze the current situation in the 
four-dimensional case. 

Let us recall that a close relation between loops and spin foams follows already from the very 
basic features of these two approaches. Both the kinematical states in loop quantum gravity 
(LQG) and the boundary states between which one evaluates transition amplitudes in spin foam 
models are represented by spin networks, i.e. by graphs colored with group theoretic data such 
as representations and intertwiners. Moreover, in their seminal work [134], Reisenberger and 
Rovelli provided a precise link between the different physical quantities associated to these spin 
network states. Namely, they showed that an expansion in powers of a Hamiltonian operator of 
the physical scalar product between two kinematical states in LQG generates a sum over spin 
foams with weights given by matrix elements of this Hamiltonian. 

In principle, this result suggests a direct way to construct a spin foam model fully consistent 
with LQG. In three dimensions, although the spin foam models have been derived originally in a 
different way, their link with the canonical approach is established using the same set of ideas. In 
fact, there are two rather different frameworks in which this can be done. The first one relies on 
the work of Witten [150], who proposed the first complete and consistent quantization of three- 
dimensional Riemannian gravity with a positive cosmological constant. In this quantization, the 
path integral is represented by the scalar product between physical states. On the other hand, 
Turaev and Viro [146] proposed a state sum representation of the same path integral, which 
was the first mathematically well-defined spin foam model. The second framework is provided 
by Riemannian gravity with vanishing cosmological constant. The corresponding spin foam 
quantization was found by Ponzano and Regge [132], and then in [121] it was shown that the 
Ponzano-Regge amplitudes are reproduced by evaluating the physical scalar product between 
spin network kinematical states of the loop quantization. 

Thus, in the case of three-dimensional gravity (at least for Riemannian signature) there are 
already solid results confirming the agreement between the spin foam and canonical quantiza- 
tions. In contrast, the agreement in the four-dimensional case is less obvious. The difficulties 
associated with the Hamiltonian operator and the evaluation of its matrix elements have led 
people to reject the strategy which consists in finding the spin foam representation from the 
canonical side. Instead, the spin foam quantization has developed its own, completely indepen- 
dent methods. As a result, although at the qualitative level the models derived in this way 
produced a picture similar to LQG, quantitatively for a long time there was a striking disagree- 
ment. Since spin foam models represent a completely covariant quantization of gravity in the 
first order formulation, their boundary states are labeled by the group theoretic data of the 
unbroken gauge group in the tangent space, which is either the Lorentz group or Spin(4) de- 
pending on the signature. On the other hand, LQG is based on the SU(2) group obtained by 
a partial gauge-fixing at the classical level. A longstanding problem was to realize how a similar 
reduction to the compact subgroup can take place in spin foam models. 

The progress came from two sides. On one hand, it was understood how to represent LQG in 
a Lorentz-covariant form [15]. In particular, this brought into play the so-called projected spin 
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networks, which constitute a certain generalization of the usual spin networks. As it turned out, 
they are equally relevant for the spin foam approach since the boundary states of all existing 
spin foam models of four-dimensional gravity can be expressed as their linear combinations. On 
the other hand, some new spin foam models were introduced [66, 75], with the hope that they 
could potentially overcome various problems encountered with the previous ones. Furthermore, 
it was argued that the new models implement the reduction to the SU(2) gauge group mentioned 
above, which is a necessary step if one requires agreement with the canonical quantization. These 
findings suggest that we are now in a much better position than just a few years ago, and that 
the issue of the mutual consistency might be soon (if not already) resolved. 

In this review, we reconsider this issue in detail and provide an analysis of the relationship 
between the spin foam and canonical quantizations in three and four spacetime dimensions. We 
start by presenting the general ideas of the spin foam approach in Section 2. In particular, we 
recall its relation to the physical scalar product of the canonical loop quantization established 
in [134] and explain the main strategy used to derive spin foam models in four dimensions. 

In Section 3 we concentrate on Riemannian three-dimensional gravity. After briefly review- 
ing the classical aspects of pure gravity and gravity coupled to point particles, we first recall 
the quantization a la Witten [150], which relies on the Chern-Simons formulation of Rieman- 
nian gravity with a positive cosmological constant. In this context, we show how the physical 
Hilbert space can be constructed using the so-called combinatorial quantization and how the 
physical scalar product is related in a very precise way to the Turaev-Viro spin foam model. 
Then we turn to Riemannian gravity without cosmological constant. After reviewing its loop 
quantization and the Ponzano-Regge spin foam model, we show, following [121], how the spin 
foam representation of the path integral can be obtained from the kinematical Hilbert space of 
the canonical quantization by projecting states onto the physical Hilbert space. This provides 
an example of the concrete realization of the ideas of [134]. We also present some attempts to 
generalizing this result to the case of a positive cosmological constant. 

In Section 4 we study the four-dimensional case. First, we review the basic notions of LQG 
and its Lorentz-covariant reformulation. For this purpose, we introduce the projected spin 
networks, which turn out to play also an important role in spin foam models. Second, we re- 
call the canonical structure of the Plebanski formulation of general relativity which is at the 
heart of the spin foam approach in four dimensions. After this, we define the most relevant 
spin foam models studied in the literature (BC [32, 33], EPRL [66] and FK [75]) and present 
their main features. The EPRL model is then compared in detail to LQG. This comparison 
suggests some exciting relations, but at the same time a thorough analysis opens various prob- 
lems challenging some of the conclusions which have been taken for granted. In particular, 
it reveals a tension between the derivation of the new models and the canonical quantization. 
Taking seriously the lessons obtained in the canonical theory, we forget about the existing mod- 
els and discuss restrictions imposed on the spin foam quantization in four dimensions by the 
requirement of consistency with the canonical approach. In particular, we address the issues 
of the imposition of constraints, the path integral measure and the form of the vertex ampli- 
tude. 

Then, in Section 5, we describe a new model proposed recently in [88]. Although this model is 
three-dimensional, it is designed to test the implementation of the simplicity constraints appear- 
ing in the four-dimensional case. Its analysis shows that implementing the constraints following 
the method of either the BC or the EPRL model, one does not arrive to the expected vertex am- 
plitude. Instead, the correct result can be recovered by incorporating the ideas presented at the 
end of Section 4 and modifying the standard prescription for evaluating the vertex amplitudes 
from the simplex boundary states. In the construction of [88] however, the weak imposition 
of the simplicity constraints is an essential ingredient which enables to recover the expected 
amplitude. The final section contains our conclusions. 
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Some comments about our notations. The signature of spacetime will be encoded in a sign 
function which is a = 1 in the Riemannian case and a = — 1 in the Lorentzian case. The corre- 
sponding isometry algebras of the flat Euclidean or Minkowskian space, so(4) and so(3, 1), both 
admit su(2) as a subalgebra. Notations are such that //, v, . . . refer to spacetime indices, a,b, . . . 
to spatial indices, /,«/,... to so(4) or so(3, 1) indices, i, j, ... to su(2) indices, and (• •) (or [• •]) 
to the (anti-)symmetrization taken with weight 1/2. Tr and tr will denote traces over so(4) or 
50 (3, 1) and su(2) indices, respectively. The star * is the Hodge operator acting on bivectors (or 
equivalently in the adjoint representation of the isometry algebra) as {-kB) IJ = ~ e IJ klB kl . 
In the spin foam context we will use <t, t, t to label four-simplices, tetrahedra and triangles of 
a triangulation A; v, e, f for vertices, edges and faces of the dual triangulation A*; and re, t 
for nodes and links of the boundary graphs T. For the quantities which have a meaning for 
both, say A* and T, the indices may be used interchangeably. For example, the holonomy as- 
sociated to a face / bounding the link £, can be denoted either as gj or Note also that the 
intertwiners are considered in this review as generic tensors in the tensor product of represen- 
tations, not necessarily invariant. Instead, invariant tensors are usually specified explicitly as 
invariant intertwiners. Finally, we will assume that the d-dimensional spacetime manifold Ai is 
topologically SxR, where £ is a (d — l)-dimensional manifold without boundaries. 



2 Spin foams — general concepts 

In this introductory section, we expose the main concepts behind the spin foam approach to 
quantum gravity. We first give the definition of a spin foam and a spin foam model, then explain 
the qualitative relations between the spin foam and canonical approaches, and finally present 
the general strategy used to derive spin foam models in three and four spacetime dimensions. 

2.1 Structure of spin foam models 

A spin foam T = (A*,A/,Z e ) is defined as a two-dimensional cellular complex 1 A* colored 
with group representations Aj on its faces and intertwiners X e on its edges. One can think 
of this coloring as some kinematical information inherited from the boundary states, which 
are one-dimensional cellular complexes (graphs) obtained by intersecting the spin foam with 
codimension-one surfaces. These graphs are colored with representations Xe on their links and 
with intertwiners I n on their nodes. If we think of them as encoding quantum spatial geometry, 
the analogy allows us to see a spin foam as a representation of quantum spacetime. In this 
way, the sum over different spin foams can be viewed as a discrete path integral over possible 
geometries interpolating between given boundary information. Note that in this picture the 
dynamics is encoded at the vertices of the spin foam, which therefore carry no labels. 

A spin foam model is defined by an assignment of complex amplitudes Af (A/), A e (Xf De ,I e ) 
and A v (Xf Dv ,Z eDv ) to the faces, edges and vertices of T respectively. These amplitudes are 
local in the sense that Af depends only on the representation coloring the face /, A e on the 
representations coloring the faces / D e meeting at e and the intertwiner coloring e, and A v on 
the representations and intertwiners carried by the faces / D v and edges e D v meeting at the 
vertex v. Given two boundary states Si and 52, and a two-complex A* interpolating between 
their graphs, one can define the transition amplitude 

w(a*) = (s 1 \s 2 )a* = I d^ {Xf} f d / u {Xe} n^/(A/)n^( A /'^)ri^( A /'^)- w 

f e 



: A two-dimensional cellular complex is a combinatorial structure consisting of faces / meeting at edges e, in 
turn meeting at vertices v. 
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Figure 1. A spin foam defined on a two-complex interpolating between two boundary states Si and <S 2 - 
A section defines an intermediate state in the middle of the evolution. 

For the sake of generality, we have written this amplitude with arbitrary integration measures, 
d^| A/ } and d/j,^ Xe }, depending on the group-theoretic data coloring the two-complex. If the spin 
foam model is built with a compact Lie group, these continuous integrals become discrete sums 
over all possible colorings with representations and intertwiners of all the internal faces and 
edges of the two-complex. 

Two comments are in order. First of all, being similar to quantum amplitudes associated 
with Feynman diagrams in quantum field theory, the spin foam amplitude (2.1) may diverge. 
Given this analogy, it is natural to ask the following question: if the amplitude does diverge, can 
we understand the nature of this divergence and find a procedure to regularize it? The answer 
will turn out to crucially depend on the dimension of spacetime and on the particular model 
under consideration. 

Second, the amplitude (2.1) represents a contribution of only one particular two-complex. 
Recalling again the analogy with Feynman diagrams, it is natural to expect that the full tran- 
sition amplitude should be given by a sum over cellular complexes interpolating between the 
boundary graphs. The same expectation is suggested by the relation to canonical quantization 
presented in the next subsection. However, in the important case of topological theories, it turns 
out that the contribution of one spin foam is sufficient. In particular, we are going to see that 
this is the case for general relativity in three dimensions, where the amplitude (2.1) becomes 
independent of the choice of cellular decomposition. This was in fact the original motivation 
for introducing spin foam models from the mathematical point of view, since it enables one to 
compute topological invariants of three-manifolds. 

Now, in order to understand the geometric meaning of a spin foam, let us consider a simplicial 
decomposition A of the four- dimensional spacetime manifold Ai. Every such triangulation 
appears to be dual to a cellular two-complex A* whose vertices are all five-valent. Indeed, 
the dual two-complex can be constructed in the following way. To each four-simplex of A, we 
assign a point in its interior, which corresponds to a vertex v of A*. Then, to each of the 
five tetrahedra corresponding to the boundary of this four-simplex, we assign an edge e E A* 
which is connected to the vertex v. Finally, to each of the four triangles corresponding to 
the boundary of a tetrahedron, we assign a face / in A* whose boundary is the edge e. Spin 
foams with vertices of valence higher than five (or those whose boundary graphs have nodes of 
valence higher than four) can be obtained from discretizations which are more general than the 
simplicial ones considered above. The spin foam represented in Fig. 1 is defined on a simplicial 
triangulation of a three-dimensional manifold. 

Given this duality between spin foams and triangulations, one can view a spin foam model as 
a proposal for a regularized version of the path integral for gravity. Namely, given a manifold A4 
with boundaries Si and S2, and a diffeomorphism class of metrics q\ on Si and qi on S2, the 
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spin foam model is supposed to give a meaning to the transition amplitude 

((Ei,g 1 )|(E 2 ,?2))phys = / ^9 exp(iS), 

where the left-hand side denotes the physical scalar product in the yet to be defined physical 
Hilbert space of the theory. The integration is performed over the diffeomorphism class of 
metrics g on A4 which agree with the metrics q\ and q2 on the boundaries Si and £2, respectively. 
The idea is therefore to give a meaning to the physical inner product between boundary states 
of quantum geometry. 



2.2 Relation to canonical quantization 

We saw that the boundary states of spin foam models are represented by graphs colored with 
representations and intertwiners of a certain group. Such objects are known as spin networks 
and generically appear as the so-called kinematical states in the canonical loop quantization 
of first order gravity 2 . These states solve the constraints generating gauge transformations 
in the tangent space and spatial diffeomorphisms. However, it remains to solve the so-called 
Hamiltonian constraint which implements the dynamics of the theory. This can be done by 
mapping the kinematical states into the kernel of the quantum Hamiltonian operator H. As was 
first realized in [134], this directly leads, at least formally, to the spin foam picture. We recall 
here the main steps of this construction. 

The idea of Reisenberger and Rovelli is to construct the projection operator using the path 
integral representation of the ^-distribution, 

V H = J] = fvN exp [iH(N)), (2.2) 

where 



H(N) = J d 3 xN(x)H(x). 



If we forget for the moment about the mathematical obstructions to defining such an operator, 
we can think of this formula as a group averaging map which enables to define the physical 
inner product as 

where Si and S2 are two kinematical states. Using an expansion of the exponential in (2.2), this 
inner product becomes 



r I 00 -n 

(Si|7V?2>kin = jVN (S x Y,-(H(N)) n 




(2.3) 



where the right-hand side is defined as n actions of the Hamiltonian constraint on the state £2- 
A key property of the Hamiltonian constraint is that it acts only at the nodes of a graph, creating 
new links and nodes, and thereby changing the geometry encoded in the state. Therefore, we 
can rewrite (2.3) in the following form 



n=0 T n :Si^S 2 v 



2 Here we discuss only the general structure. The precise relationship between groups, representations and 
intertwiners coloring the graphs in the canonical and spin foam approaches will be discussed in the following 
sections. 
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where the (divergent) integral over the lapse function was formally included in the definition of 
the amplitude A v . In this formula, the second sum runs over the spin foams T n with n vertices 
interpolating between the two boundary states. The contribution of each T n corresponds to n 
actions of the Hamiltonian constraint on these states, each creating a vertex v and producing 
a factor given by A v . 

Given this construction, it is clear that the kinematical structure of canonical quantum grav- 
ity, together with a suitably regularized Hamiltonian constraint, should lead to a certain spin 
foam model. However, the difficulties in dealing with the quantum Hamiltonian constraint 
(see [45] for a recent overview) have forced people to address the problem from the other side. 
Namely, one usually starts by constructing a spin foam model along the lines explained in the 
following subsection, and then tries to compare it with the discrete Hamiltonian dynamics built 
from the canonical theory. Thus, on a general level, the canonical and covariant approaches 
should be considered as two completely different paths. On the other hand, if at some point 
the results obtained in the spin foam context and in the canonical quantization do converge 
in the same direction, this might be considered as a non-trivial consistency check of the two 
approaches. 

2.3 Strategy 

Spin foam models for quantum gravity were originally introduced in three spacetime dimensions. 
In this case, pure gravity is a topological field theory, i.e. it has no propagating local degrees 
of freedom (no gravitons). It is one of the reasons why it is possible to discretize the theory in 
a way which does not depend on the choice of discretization. Because of this crucial property, 
the amplitudes (2.1) do not depend on A*, and the sum over different two-complexes (for a fixed 
topology) becomes irrelevant 3 . The first two models for quantum Riemannian general relativity 
which were built along these lines are the Ponzano-Regge model [132] (which has a vanishing 
cosmological constant) and the Turaev-Viro model [146] (which has a positive cosmological 
constant). We will review this construction in the next section. 

In four spacetime dimensions, the situation is drastically different, since the theory is not 
topological anymore. However, it is possible to exploit the fact that, in any dimension, gravity 
can be written as a constrained topological field theory. The topological theory itself, known 
as BF theory and introduced for the first time in [97], is built using an so(4)- or so(3, l)-valued 
{d — 2)-form B and the curvature two- form F of a connection u. Its action is given by 



The equations of motion, DB = (where D is the covariant differential associated to oS) and 
F = 0, imply that all the solutions are locally pure gauge. General relativity (in dimensions 
d > 3) is obtained from BF theory by requiring that the B field is not arbitrary, but rather 
constructed from vierbein one- forms e 1 . In four dimensions this implies that it is constrained to 
satisfy the relation 



in which case (2.4) becomes the usual Hilbert-Palatini action. To achieve this, one can modify 
the BF action to obtain the Plebanski action [51, 58, 131] 



3 However, the question of summing over topologies is very interesting but still unsolved. The group field 
theory formulation of spin foam models offers a nice framework to study this question [72, 124]. 




(2.4) 
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The field has 20 independent components generating constraints which ensure that the B 
field comes from the wedge product of two tetrad fields. These are the so-called simplicity 
constraints, reducing the topological field theory to general relativity (plus additional sectors 
that we will discuss later on). 

Most of the spin foam models for four- dimensional quantum gravity rely on the Plebanski 
formulation and use the same strategy, but translated to the quantum level. Namely, instead 
of quantizing the phase space of the full Plebanski action, obtained by imposing the simplicity 
constraints at the classical level, one first discretizes the classical theory on a cellular complex, 
quantizes its topological part given by the discrete BF action, and finally imposes the simplicity 
constraints at the quantum level. The last step is implemented via certain conditions on the 
group-theoretic data associated to the spin foam quantization of BF theory, and is supposed to 
convert the trivial dynamics of a topological spin foam model into the dynamics of quantum 
general relativity. 

3 Three-dimensional quantum gravity 

Three-dimensional classical general relativity admits two equivalent first order formulations. In 
the BF formulation, the dynamical variables are a triad e and an h-valued spin connection u, 
where the Lie algebra f} of the gauge group is su(2) in Riemannian signature and su(l, 1) in 
the Lorentzian case. In the Chern-Simons formulation, the dynamical variable is a g-valued 
connection A, where g is a Lie algebra which is totally and uniquely determined by the signature 
of the spacetime and the sign of the cosmological constant (which will be denoted by A in what 
follows) according to the following table: 





A < 


A = 


A > 


a = 1 


80(3,1) 


isu(2) 


so (4) 


a = -1 


80(2,2) 


isu(l, 1) 


so(3,l) 



In spite of the apparent differences, the BF and Chern-Simons formulations are equivalent 
at the classical level. More precisely, the Chern-Simons connection A appears as a linear combi- 
nation of the triad one- form e and the spin connection cu. Furthermore, the gauge algebra of the 
Chern-Simons theory is nothing but the symmetry algebra of the associated BF formulation. 
This comes from the fact that a three-dimensional BF theory, associated to a gauge algebra f), 
admits a larger symmetry algebra whose nature depends on the sign of the cosmological constant 
according to the classification that we have given above. In all cases, q contains t) as a sub- 
algebra and the remaining symmetries describe translations in the coset space G/H, where G 
and H are universal coverings of the Lie groups associated to g and h, respectively. 

As we have already argued, the BF and Chern-Simons formulations posses the same space 
of classical solutions. However, in spite of this equivalence at the classical level, the two theo- 
ries naturally lead to different quantization schemes for three-dimensional gravity. If it is clear 
that both formulations can be quantized canonically and covariantly, it turns out that the co- 
variant quantization is much simpler with the BF formulation, whereas Chern-Simons theory 
offers a very powerful framework for the canonical quantization (regardless of the sign of the 
cosmological constant). To make this point more concrete, one can develop an analogy with the 
quantization of the harmonic oscillator in quantum mechanics. Indeed, the harmonic oscillator 
admits two well-known first order formulations. In the first one, the dynamical variables are 
the position q and the momentum p, whereas in the second one the dynamical variables are the 
chiral variables a and a', which are in turn complex linear combinations of q and p. The position 
and momentum variables are clearly more suited for a covariant quantization, whereas it is much 
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simpler and much more powerful to use the chiral ones for the canonical quantization. Thus, 
e and lv can be seen as being equivalent to q and p, whereas the Chern-Simons connection can 
be interpreted as being analogous to a. As a consequence, establishing a direct link between the 
covariant and canonical quantizations is a priori a non-trivial task. Nevertheless, this has been 
achieved for the Chern-Simons formulation of three-dimensional gravity, and led to amazing and 
deep relationships between topological invariants, knot invariants, conformal field theories and 
quantum gravity. In particular, these discoveries have opened a new way towards the under- 
standing of knot polynomials. On the other hand, the link between the covariant and canonical 
quantizations of BF theory leads to a clear understanding of the relation between LQG and spin 
foam models in three dimensions (at least in the case of a vanishing cosmological constant). 

We are going to present these aspects in much more details in this section. The first subsec- 
tion is devoted to the classical theory. We will present the two first order formulations mentioned 
above, perform their canonical analysis and discuss the inclusion of point particles. In the sec- 
ond subsection, we will focus on the canonical quantization of Chern-Simons theory using the 
so-called combinatorial quantization scheme. This method is very powerful and allows to un- 
derstand clearly how quantum groups appear in three-dimensional quantum gravity. Then, we 
will review the covariant quantization a la Witten and explain how it provides a bridge between 
topological and knot invariants on one side, and Riemannian quantum gravity with a positive cos- 
mological constant on the other side. The third subsection is devoted to the quantization of BF 
theory. We will mainly concentrate on the case where the gauge group is G = SU(2) and where 
there is no cosmological constant. In this framework, we will present the canonical loop quanti- 
zation and derive the spin foam representation as a realization of the Reisenberger-Rovelli path 
integral. This will enable us to construct the physical scalar product and to show the relation- 
ship between the canonical quantum theory and the Ponzano-Regge spin foam model. We will 
present some observations about the theory in the presence of a positive cosmological constant. 

3.1 Classical theory 

In the standard second order formulation, gravity is a theory of the metric of the spacetime 
manifold M and its dynamics is governed by the Einstein-Hilbert action 



where d 3 x y/\g\ is the volume form on M., A is the cosmological constant, R is the scalar cur- 

vature, and Gn is the three-dimensional Newton constant. As in four dimensions, the resulting 
Euler-Lagrange equations of motion are the well-known Einstein equations 



However, the solutions are very different from the solutions in the four-dimensional case. Indeed, 
in three dimensions, the components of the Riemann tensor are linear combinations of the 
components of the Ricci tensor, and therefore the Weyl tensor vanishes identically. But the Weyl 
tensor is responsible for the presence of local degrees of freedom, and in particular gravitational 
waves in four dimensions. As a consequence, there are no gravitational waves in three dimensions, 
there are no local degrees of freedom, all solutions are locally of constant curvature, and the only 
non-trivial degrees of freedom are topological. In other words, the spacetime manifold has to be 
topologically non-trivial for the theory to admit non-trivial and physically interesting solutions. 

The absence of local degrees of freedom gives a priori the feeling that three-dimensional 
gravity is too trivial to be interesting in itself, and cannot be a good toy model to test ideas and 
principles of classical and quantum gravity in four dimensions. But this viewpoint is too naive, 
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firstly because the quantization of three-dimensional gravity is far from being trivial, and in 
fact turns out to be mathematically extremely rich, and also because it has proven to be a very 
insightful toy model to understand important aspects of quantum gravity in four dimensions [52]. 
For instance, this is the context in which spin foam models have originally been introduced. 

In this subsection we perform the classical analysis of the two first order formulations of three- 
dimensional gravity mentioned previously, which we take as a starting point to explain the main 
quantization schemes. In a first part, we will present the BF formulation, and in a second part 
the Chern-Simons formulation. Then, we will study the inclusion of point-like particles, which 
are crucial ingredients to understand the link between quantum gravity and knot invariants like 
the Jones polynomial [98]. In the last part, we will perform the canonical analysis and construct 
the classical physical phase space in both formulations, with and without particles. 

3.1.1 The BF formulation 

In the BF formulation of (Riemannian) three-dimensional gravity, the basic field variables are 
a triad e and a spin connection cj. Both fields can be viewed as su(2)-valued one-forms in the 
sense that e = e^Tidx^ and ui = uj^Tidx^ , where {rj} forms a basis of su(2) defined by the 
commutation relations [t^Tj] = eij k Tj~, with the totally antisymmetric tensor such that 
£123 = +1- Indices are raised and lowered with the flat metric 5ij. In the spinorial (two- 
dimensional) representation, the su(2) generators are obtained from the Pauli matrices {a{\ via 
Ti = — io"j/2, and in the rest of this work we will often identify the abstract generators Tj with 
their two-dimensional representation. 

In the BF formulation, the triad variables are related to the spacetime metric by the relation 
9fiv = e^Sijei, and the spin connection is taken to be an independent field. When expressed in 
terms of these variables, the Einstein-Hilbert action reduces to the topological BF action 



where tr denotes the trace in the spinorial representation of su(2) and F = Duj = doo + wAw 
is the curvature two- form corresponding to the connection u. Varying this action with respect 
to the triad and the spin connection leads to the following two equations: 



The second equation is the torsion-free condition and is equivalent to the requirement that 
the triad e be covariantly constant with respect to uj. When the triad (viewed as a square 
three-dimensional matrix) is invertible, the torsion-free equation implies that u is the unique 
connection compatible with the triad which can be expressed explicitly in terms of e. Substi- 
tuting the expression w(e) into the BF action (3.3), one ends up with the initial second order 
Einstein-Hilbert action (3.1). This establishes the equivalence between the BF formulation and 
the second order formulation of three-dimensional gravity, when the triad is assumed to be 
invertible. 

Note that some of the classical equations (3.4) do not involve time derivatives and therefore 
should be considered as constraints on the dynamical variables. As we are going to see in the 
Hamiltonian framework, these constraints are first class and generate the gauge symmetries of 
the theory. Understanding these symmetries is essential in order to describe properly the physical 
solutions. However, we postpone the canonical analysis of BF theory after the introduction and 
the study of the Chern-Simons formulation of gravity for the sake of clarity and simplicity. 




(3.3) 
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3.1.2 The Chern Simons formulation 

To go from the BF formulation to the Chern-Simons formulation of three-dimensional gravity, 
one combines the triad and the spin connection into a single gauge field: the Chern-Simons 
connection A. To see how this works explicitly, it is necessary to introduce a family of Lie 
algebras qa labeled by the cosmological constant A. gA is an extension of the algebra f) = su(2), 
in the sense that it is generated by the three generators of su(2), denoted {r^}, supplemented 
by three extra generators {Pi}- They satisfy the following commutation relations: 

[Ti,Tj] = e i:j k T k , [Pi, Pj] = -sgn(A)£i/r fc! [n, Pj] = £ij k P k , (3.5) 

where sgn(x) is the sign function defined such that sgn(O) = 0. These relations define three 
different algebras characterized by the sign of the cosmological constant and given in the first 
line of the table on p. 9. These algebras admit a Killing form defined by the pairing 

(•,•): BAX0A — {P i ,T j ) = 5 ij , (3.6) 

all the other pairings being trivial 4 . 

These algebras have an immediate geometric interpretation as the isometric algebras of the 
three-dimensional spaces of constant scalar curvature. More precisely, isu(2) is the isometry 
algebra of the three-dimensional Euclidean space E 3 (of vanishing scalar curvature), so (4) is the 
isometry algebra of the three-dimensional sphere § 3 (of positive scalar curvature), and so (3, 1) is 
the isometry algebra of the three-dimensional hyperbolic space H 3 (of negative scalar curvature) . 
Locally, all the solutions to the Einstein equations (3.2) look like one of these three spaces, de- 
pending on the sign of the cosmological constant. It is therefore not a surprise to see that the cor- 
responding isometry algebras will play an important role in three-dimensional gravity. This role 
is particularly important and explicit in the Chern-Simons formulation. From now on, in order to 
simplify the notations, we will omit the index A in the Lie algebras 0a when there is no ambiguity. 

The Chern-Simons connection is defined by 

A = wVi + y^Pi, 

where the cosmological length t c is a function of the cosmological constant A given by £ c (0) = 1 
and £ C (A) = |A| -1 / 2 . The gauge field A is a one-form with values in the Lie algebra g. It is 
shown in [148] that the first order action for three-dimensional gravity can then be rewritten as 
a Chern-Simons action 

S CS [A(e,u)] = ^-f^xe^ [{A, } d v A p ) + ^(A^[A u ,A p ])^ , (3.7) 

where (• , •) is the bilinear form (3.6) and k = icG^ 1 is called the level of the Chern-Simons 
action. Using (3.5), it is easy to check that this action is indeed equivalent to (3.3) up to 
a boundary term for dA4 ^ 0, which does not affect the equations of motion. 

Before going any further, let us discuss in more detail the change of variables to go from the 
BF variables e and oj to the Chern-Simons gauge field A. As we explained in the introduction to 
this section, this is not a simple change of variables, but rather a kind of "complexification" of 
the initial fields. Indeed, everything works as if we had extended the field (in the mathematical 
sense) of rotations to a "complex" field which involves the rotations but also the translations. 
In other words, we have doubled the number of dimensions of the initial field. Then, we have 
combined these two types of variables, e and u>, into a single "complexified" connection A whose 



4 Notice that we use the notation Tr and tr exclusively for the bilinear pairing of so (4) or su(2) elements 
respectively. 
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rotational part is the spin connection, and whose translational part (in the complex direction) 
is the triad. In this sense, A is analogous to the chiral variable a for the harmonic oscillator, 
this later being related to the phase space variables q and p by the complexification a = q + ip. 

Varying the Chern-Simons action (3.7) with respect to the gauge field results in a flatness 
condition on the 3-valued Chern-Simons connection A. As expected, this flatness condition 
on A combines the torsion- free and constant curvature conditions for the spin connection uj: 



F^[A] = 



F^ v [uj] = -A[e M ,e„], 
T^ v [e,uj] = 0. 



The consequences of this unification of the variables e and uj into a single gauge field A are 
manifold. First, the Euler-Lagrange equations are encoded into a single, simple and beautiful 
equation, which is the requirement of flatness for the connection A. Second, the gauge symmetry 
group is now appearing explicitly, whereas in the BF formulation a part of these symmetries was 
somehow hidden. Third, it makes obvious that the theory admits no local degrees of freedom 
since it is always possible to find a gauge transformation that sends locally every flat connection 
to the trivial connection. 

As the Chern-Simons formulation of three-dimensional gravity is a gauge theory with local 
symmetry group G, its action admits an infinite-dimensional symmetry group given by Q = 
C°° (M. , G) . This is the group of G- valued smooth functions on M. , and its action on the 
connection reads as 

VflSe A^ A3 =g- 1 Ag + g- l dg. (3.8) 

The invariance of the Chern-Simons action (3.7) with respect to the infinitesimal version of 
the transformations (3.8) is an immediate consequence of the invariance of the bilinear form 
(• , •) under the adjoint action of G. To see how these transformations encode diffeomorphisms, 
let us consider an infinitesimal version of the transformation (3.8) generated by <p € 0. If we 
decompose the Lie algebra element 4> into its translational and rotational parts according to 
4> = u l Pi + uVi, we can use the commutation relations (3.5) to show that e and uj transform 
according to the following rules: 

be = £ c (du + [uj, u]) + [e, v], 5uj = dv + [uj, v] + y[e, u]. 

Setting u = £^e^ and v = ^u„, we can then express these transformations in terms of the Lie 
derivatives C% along the vector field £ = ^d^. This leads to 

Sep = C^e^ + i u T^ u [e,uj], Suj^ = L^uj^ + ^(F^u) + A[e M , e„]). (3.9) 

As a result, when the gauge fields satisfy the equations of motion, the infinitesimal gauge trans- 
formations (3.9) are in fact given by diffeomorphisms, which establishes the on-shell equivalence 
between infinitesimal diffeomorphisms and infinitesimal Chern-Simons gauge transformations. 
However, this equivalence applies only to gauge transformations and diffeomorphisms that are 
connected to the identity, whereas the status of large (i.e. not infinitesimally generated) gauge 
transformations is more subtle [89, 114]. 

In fact, the Chern-Simons action is not invariant under the large gauge transformations (3.8). 
Instead, one has the relation 

Scs[A^-S cs [A]=w[M,g) = ^- I (g~ x dg A g^dg A g~ l dg). (3.10) 

6vr J M 

The quantity on the right hand side, which is known as the Wess-Zumino-Witten term, is 
a priori not null. It is a topological invariant of the principal G-bundle over j\A, and as such 
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it is a constant which does not modify the classical Euler-Lagrange equations. In this sense, it 
is fair to claim that the Chern-Simons action is invariant under gauge transformations. Note 
that when the gauge group is compact, it is possible to show that the topological invariant is 
= 2nNk, where N is an integer and k is the level of the Chern-Simons action. 

3.1.3 Inclusion of particles and Wilson loops 

In three-dimensional spacetime, there is no gravitational interaction between massive particles. 
Indeed, from a classical (non-relativistic) point of view, it is possible to show that massive 
particles do not create any gravitational field in three dimensions. From the point of view of 
general relativity, particles do not deform the spacetime manifold. Locally, in the vicinity of the 
particle, spacetime remains flat (if A = 0) or of constant positive (if A > 0) or negative curvature 
(if A < 0). This property contrasts drastically with what happens in four dimensions, where the 
matter fields have in general a strong effect on the geometry of spacetime. Therefore, at first 
sight, it might seem useless to study particles in three dimensions since they have no effect. 

However, it turns out that a particle in three dimensions has a non-local effect on space- 
time [60, 61]. It does not modify the geometry but the topology of spacetime. More precisely, 
at any time, a massive (non-spinning) particle creates a conical singularity at the place where 
it is located, so that space looks like a cone around the particle. The angle a of the cone is 
proportional to the mass m of the particle in Planck units. In the presence of several particles, 
space looks like a multi-cone and its dynamics is far from being trivial [111]. For instance, 
there is a close and deep relationship between the properties of the diffusion of particles and 
the so-called exotic statistics [26]. Furthermore, when the particles possess a non-trivial spin, 
the topology becomes even more complicated. In particular, time-helical structures appear and 
lead to the presence of closed timelike curves [91]. 

These examples show that in the end the dynamics of point particles in three dimensions is 
quite complicated, even if there are just modifications of spacetime topology, and not geometry. 
Unfortunately, we cannot review all these features in the present work. Instead, we are going 
to focus only on the aspects which are directly related to our subject. The reason why we are 
interested in studying particles in three-dimensional quantum gravity is that they are intimately 
related to the Jones polynomials. To see this, we first need to understand how to describe 
a point particle coupled to gravity at the classical level. This can be done using the so-called 
algebraic action. Then the coupling to the gravitational field reduces to gauging the global 
symmetries of this action with the Chern-Simons gauge field. 

The algebraic formalism is adapted to describe a free point particle evolving in any homo- 
geneous three-dimensional spacetime. In the Riemannian three-dimensional case, the homoge- 
neous manifolds are the Euclidean space E 3 , the sphere § 3 , and the hyperbolic space H 3 . They 
correspond to the maximally symmetric solutions to Einstein equations with vanishing, posi- 
tive or negative cosmological constant, respectively. Each of these spaces can be obtained as the 
coset G/H, where H = SU(2) and G is the associated isometry group. In particular, G = ISU(2) 
for E 3 , G = Spin(4) for S 3 and G = SL(2, C) for H 3 . As the action of H is transitive, any element 
g € G can be written as the product 5 g = qp of an element q 6 G/H with an element p 6 H. 
This is summarized by the following map: 

G — > G/H X H, g = qp — ► (q,p). 

In the algebraic formulation, we consider g G G as a variable describing the dynamics of the 
point particle evolving in the corresponding homogeneous space. We identify q € G/H with the 
position of the particle in spacetime, whereas the element p £ H is viewed as its momentum. 



J Note that this factorization is not necessarily unique. 
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In this context, relativistic particles at rest are classified by the unitary irreducible represen- 
tations of the isometry group G. For any value of A, these representations are labeled by a couple 
(m, s) of real parameters fixing the values of the two Casimir operators of (the enveloping algebra 
of) q in the given representation: 

Cg } = PVF4 - sgn(A)rVi = m 2 - sgn(A)s 2 , = FS* = ms. 

Even if the spacetime is taken to be Riemannian here, we will interpret m as the mass and s as 
the spin of the particle. When G is compact, which corresponds to G = Spin(4) in the present 
case, the two parameters are half-integers. If G is the non-compact group ISU(2) or SL(2,C), 
the mass m is any real number and the spin s is a half-integer. 

The algebraic action for a relativistic particle of mass m and spin s evolving in the homoge- 
neous space G/H is given by 

S[g] = / {x(m,s),g~ 1 dg), 

where x( m ) s ) = mr i + S P £ characterizes the particle at rest and 7 denotes its world line. 
Such an action has been introduced in [27] and studied in the context of Chern-Simons theory 
in [59, 64, 150]. We refer to these articles for more details. Here we will just recall the main 
results. 

• The algebraic action is a first order action for the particle since it involves both position q 
and momentum p through the variable g. 

• The equations of motion are easy to obtain and are simply given by 

= 0. (3.11) 

The general solution is g(t) = go//(t), where go is a constant and /j,(t) is any function in 
the maximal two-dimensional Cartan torus of G generated by exp(ri) and exp(Pi). These 
equations imply that q satisfies the geodesic equations on G/H. 

• The action is clearly invariant under the global left action of the isometry group on itself, 
g 1 — > g g = gg for any g £ G. Physically, this symmetry represents the invariance of the 
theory under a change of inertial frame. It implies the existence of Noether charges which 
are nothing but the momentum Pj of the particle and its total angular momentum ] i . 
These conserved quantities are shown to be Dirac observables. They satisfy the Poisson 
algebra 

{jiJj} = e ij k h, = -sgn(A)£i/j fe , {ji,Pj} = £i/p fc , 

as well as the quadratic Casimir relations 

p 2 — sgn(A)j 2 = m 2 — sgn(A)s 2 , p • j = ms. (3-12) 

This is the reason why the algebraic action correctly captures the dynamics of the free 
particle of mass m and spin s evolving in the homogeneous space G/H. 

• Besides, the action is invariant under the local two-dimensional (infinitesimal) gauge sym- 
metry g 1 — > g(ari + (3 Pi). This symmetry is closely related to the invariance of the action 
under time reparametrizations. As a consequence, the particle is determined at rest by its 
mass and its spin, and dynamically by the group element g up to this local symmetry. In 



-ldg 
dt 
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other words, the configuration spaces for the particles in three dimensions are in one to 
one correspondence with the coadjoint orbits 



of the Lie algebra q. In this formulation, g is the group element which sends the dynamical 
particle characterized by g to its rest frame. Furthermore, the coadjoint orbits of q are 
canonically endowed with a symplectic structure given by the Kirillov-Poisson bracket 
which coincides with the symplectic structure following from the Hamiltonian analysis of 
the algebraic action. 

Coupling the relativistic particle to the gravitational field simply amounts to gauging the 
global symmetry (i.e. the left action of the group G) of the free particle action by the Chern- 
Simons gauge field A. This is possible because the connection A is valued in the Lie algebra g 
of the same group G. As a consequence, the particle is minimally coupled to gravity and the 
action of the coupled system reads 



Note that one can always add or absorb the factor fc/47r in front of the action by a simple 
redefinition of the mass m and the spin s. It is also easy to rewrite this action in the BF 
formulation. To this end, one replaces A by its expression in terms of e and u and, up to 
a boundary term, the action becomes 



where g appears only in the definition of the momentum p l and the total angular momentum j*, 
which are the Noether charges of the free particle satisfying the quadratic Casimir relations (3.12) 
and defining its mass m and spin s. 

The equations of motion for the degrees of freedom g of the particle are easily obtained by 
replacing the derivative dg in (3.11) by the covariant derivative Dg = dg + Atg, where At is the 
component of the connection in the direction of 7. At is not a dynamical variable of the theory 
in the sense that the coupled action does not involve time derivatives of At. It is a Lagrange 
multiplier which can be set to zero with a gauge choice. As a result, the particle coupled to 
gravity evolves as a free particle, which is consistent with the fact that particles are not sensitive 
to the gravitational field in three-dimensional spacetimes. 

3.1.4 The classical phase space 

On manifolds of topology M = Sxl, one can give a Hamiltonian formulation of the theory. For 
simplicity, we first focus on the case where £ is an oriented two-surface of arbitrary genus. The 
case of a surface with punctures representing massive spinning particles is a straightforward 
generalization which will be discussed in a second time. A detailed discussion can be found 
in [52] and in references therein. 

Canonical analysis of pure gravity 

Decomposing the gauge field A = A^dt + A a dx a into a time component Aq and a gauge field 
Ay, = A a dx a on the spatial surface, we can rewrite the Chern-Simons action (3.7) as 



m,s — 



{g = gx{m,s)g 1 I g E G) C 



(3.13) 
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where e ab = e 0ab . This implies that the phase space coincides with the space of G-connections 
on the surface E, which we denote by ,A(G, E). The canonical Poisson bracket on this phase 
space is given by 

2-7T 

{AUx),A b j(y)} = T e ab ^j^j}5 2 (x-y), (3.14) 

where S {n, Pj}ij=o,i,2 are the generators of the Lie algebra q and S 2 (x — y) is the delta 
distribution on E. The time components A® of the gauge field act as Lagrange multipliers which 
impose the six primary constraints 

F 1 = e ab Fl b [A] « 0. 

It is easy to check that these constraints are first class and that the system admits no more 
constraints. In terms of the smeared curvature F(a) = J* s dPxonF 1 defined for any a = a/£j G 
G°°(E,g), their algebra is given by 

{F(ai),F(a 2 )} = —F([ ai ,a 2 ]), 

and the Hamiltonian action of these constraints on the Chern-Simons connection, 

{F(a) , A} = Da = da + [A, a] , 

reproduces as expected the infinitesimal gauge transformations. 

When expressed in terms of the BF variables e and u, the only non-trivial Poisson brackets 
in (3.14) are the ones which pair the components of the triad and spin connection, 

{el(x),ui(y)} = T 8^e ab 5 2 (x - y), (3.15) 

so that the triad e and the connection u> are canonically conjugated variables. Moreover, the 
first class constraints can be grouped into the two sets 

F = e ab F ab [oj] + Ae ab [e a , e b ] « 0, T = e ab T ab [e,oj] « 0. (3.16) 



The physical phase space 

To give a simple description of the physical phase space, it is advantageous to work with the 
Chern-Simons formulation of the theory. Let us recall that solutions to the constraints span the 
infinite-dimensional affine space of flat G-connections on E, which we denote by 

T(G, S) = {A £ A{G, E) | F[A] = 0}. 

This space inherits the Poisson bracket (3.14) from the Chern-Simons action and the action of 
the gauge symmetry (3.8). The physical phase space, denoted V(G, E), is the moduli space of 
flat G-connections modulo gauge transformations on the spatial surface S: 

V{G,T) =.F(G,E)/0e, fe = C°°(E,G). 

It inherits a symplectic structure from the Poisson bracket on J-(G, E) and, remarkably, is of 
finite dimension. More specifically, the physical phase space V(G, E) can be parametrized by 
the holonomies along curves on E and is isomorphic to the space Hom(7Ti(E), G)/G, where the 
quotient is taken with respect to the action of G by simultaneous conjugation. The physical 
observables are, by definition, functions on V(G, E). A basis can be constructed using the notion 
of spin networks on E. Alternatively, one can work with conjugation- invariant functions of the 
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holonomies along a set of curves on £ representing the elements of its fundamental group 7Ti(X). 
The simplest elements of this type are the Wilson loops 

wf ) (A)=trp)( W (A))], 

which are given by the trace of the holonomy ge(A) of the connection along the closed loop t 
evaluated in a finite-dimensional representation n^' ■ G — > End {H^') of the group G. The 
fact that the representation is finite-dimensional ensures that the Wilson loop is well-defined. 
The Poisson bracket between two such observables was first described by Goldman [90]. 

It is easy to see that the dimension of V(G, X) is given by (2g— 2) dim G, where g is the genus of 
the surface. It is in general very difficult to exhibit an explicit basis of the physical phase space. 
We know how to do it for the torus and the genus-two surface, but simple parametrizations for 
higher genus surfaces have not yet been found. As a consequence, a concrete description of the 
physical phase space is not yet available and this is an obstacle for its direct quantization, i.e. 
after the implementation of all constraints. Due to this, the only possible strategy will consist 
in quantizing before implementing the constraints. But before focusing on the quantum aspects, 
let us briefly describe how one can include point particles in this scenario. 

The physical phase space with point particles 

To account for the presence of point particles, let us consider a Riemann surface S punctured 
with p particles at the points {z±, . . . ,z p }. The particle at the point Zi has a mass mi and 
a spin Si in the sense described in the previous subsection. We know that the connection A is 
flat everywhere, except at the locations z% of the particles where it is singular. The singularity is 
completely encoded in the conjugacy class of the holonomy g^. (A) of the connection along a small 
closed loop £i around the particle (small means that the loop surrounds only the particle Zj). 
We have shown that each particle z% comes with a phase space C mi)Si (3.13), which is a coadjoint 
orbit of g endowed with the Kirillov-Poisson bracket. Furthermore, the action of the gauge 
group £/£ = C°°(S, G) on the space A(E, G) of regular connections on the surface £ extends on 
the coadjoint orbits (by the gauge principle) and acts by the adjoint action on each C muSi . This 
action is Poissonian, which means that it is compatible with the Kirillov-Poisson bracket. For 
the sake of simplicity, we keep the notation for the gauge group on the punctured surface S. 
As a consequence, the generalization of the phase space to the case of punctured surfaces is now 
immediate and is given by 

P(G, S) = ({A G A(G, S) | F[A] = 0} x C rnuSl x • • • x C rrip , Sp ) /fe. (3.17) 

This phase space is the cartesian product between the space of flat connections and the p 
coadjoint orbits, up to the simultaneous action of the gauge group It is the action of Q^, 
which "couples" the coadjoint orbits with the space of flat connections. As in the case without 
particles, the phase space is related to the space Hom(7ri(£), G) where S is now the punctured 
surface. If one denotes by li, % = 1, . . . ,p, the p generators of Hom(7Ti(£), G) corresponding to 
the loops around the punctures, then V(G, S) is isomorphic to 

{<j> e Hom(7n(E), G) | (j>(£i) E O mi;Si }/G, (3.18) 

where O mjS C G are now conjugacy classes 6 of G and the quotient is taken again with respect 
to the action of G by simultaneous conjugation. 

This completes the description of the classical theory. However, we emphasize that we did 
not propose any explicit parametrization of the physical phase space (or equivalently the space 

6 Notice that the elements of the conjugacy classes mi , Si are defined as the "exponentiated" elements of the 
coadjoint orbits C mi , Si defined in (3.13). 
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of classical solutions). Such a parametrization is very difficult to obtain in general (for any 
punctured Riemann surface) and is not very convenient for the quantization. For these reasons, 
we ended the Hamiltonian analysis with the abstract definition (3.18), which is the starting 
point of the combinatorial quantization. Instead, in loop quantum gravity one starts with the 
definition (3.17), written in terms of the BF variables. Of course, as we have already emphasized, 
the two definitions are totally equivalent, the latter being the infinite-dimensional version of the 
former coset. 

3.2 Quantum gravity in the Chern— Simons formulation 

To build the quantum theory, we are going to perform the quantization before implementing 
the gauge symmetries or, in other words, before imposing the first class constraints generating 
these symmetries. Thus, the main difficulty in the quantization process is to understand the 
symmetries at the quantum level. Exactly the same issue arises in four-dimensional loop quan- 
tum gravity, where the implementation of the quantum scalar constraint is certainly the most 
important unresolved issue. In three dimensions, we do know how this can be done at least in 
some cases. 

The first complete quantization of three-dimensional gravity was done by Witten [150], who 
found amazing and deep relationships between quantum gravity and the theory of knots and 
manifold invariants. This discovery has opened a new way towards the understanding of topolog- 
ical invariants and offered in particular a simple and beautiful framework to unify the different 
knot invariants like the Jones [98] and the HOMFLY (Hoste-Ocneanu-Millett-Freyd-Lickorish- 
Yetter) [85] polynomials. Even the more "mysterious" Vassiliev invariants [147] are incorporated 
into this approach (see [103] for a general review of these aspects). Unfortunately, we will not 
consider these aspects here in detail. It is nonetheless important to emphasize that any other 
approach to quantum gravity in three dimensions must reproduce at least the results obtained 
by Witten. In this sense, we have a good test for any candidate of quantum gravity in three 
dimensions. 

For the present discussion, the most interesting point in Witten's work is the clear relation- 
ship between the canonical and the covariant quantizations of three-dimensional gravity. It is 
indeed possible to perform both quantizations in the case of Riemannian gravity with a positive 
cosmological constant. The reason for this is that it corresponds to the only case in which the 
gauge group G of the Chern-Simons theory is compact 7 . The construction of a bridge between 
the covariant and canonical quantizations was crucial for Witten to establish the link between 
knot invariants and quantum gravity. After him, Turaev and Viro [146], and Reshetikhin and 
Turaev [135] made this link even more precise, because they settled the basis for the construction 
of the physical Hilbert space, on one hand, and proposed a state sum formulation of the path 
integral, which is known as the Turaev- Viro model, on the other hand. In their construction, 
the link between the canonical and the covariant quantizations is very clear. This can therefore 
help us to clarify the relationship between the two quantization schemes in other approaches like 
loop quantum gravity and spin foams in three and four dimensions. For all these reasons, we 
think that it is really useful to recall the basis and the main results of the quantization initiated 
by Witten and established more rigorously by Reshetikhin and Turaev. 

This section is structured as follows. First, we review the main ideas behind Witten's con- 
struction. We will see how particles are naturally associated to knots in spacetime, and how the 
computation of the path integral of Riemannian gravity with a positive cosmological constant 
coupled to particles leads to knot polynomials. In fact, the techniques involved mix both the 

7 The other Riemannian cases (vanishing and negative cosmological constant) are much more complicated. It 
remains possible to perform the canonical quantization and to construct the physical Hilbert space [50, 112], but 
the covariant quantization is still poorly understood. 
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covariant and canonical quantizations, and require a precise knowledge of the physical Hilbert 
space. The second part will be devoted to the construction of the physical Hilbert space using 
the combinatorial quantization, which is a Hamiltonian quantization of Chern-Simons theory 
where quantum groups play a central role. This method is powerful because it is very general. 
In particular, it works for all punctured Riemann surfaces and with all the gauge groups G of 
three-dimensional gravity, the only exception being SO(2,2), which corresponds to Lorentzian 
gravity with a negative cosmological constant. The third part will deal with the covariant quan- 
tization. We will show that it is possible to properly formulate the path integral of Riemannian 
gravity with a positive cosmological constant in the presence of point particles as a state sum 
model. This state sum, which generalizes the Turaev-Viro model, does not suffer from any 
divergences and is therefore mathematically well-defined. Furthermore, it is shown to give (or 
to be given by) the scalar product between physical states of quantum gravity. This property 
is known as the equivalence between Witten-Reshetikhin-Turaev and Turaev-Viro invariants. 
We will conclude with a summary and a discussion. 

3.2.1 Quantization a la Witten and topological invariants 

The goal of this quantization is to give a meaning and to compute explicitly the path integral 
for three-dimensional Riemannian gravity with a positive cosmological constant, coupled to 
point particles Mi, . . . , M p . Each particle Mj is characterized by its mass m,, its spin Sj and its 
worldline 7» corresponding to a closed loop in a Riemannian spacetime. Using the Chern-Simons 
formulation, the path integral can be formally written as 



where S[A, gj] is the Spin(4) Chern-Simons action coupled to the group elements gj which are the 
first order dynamical variables of the particles. Due to fact that Spin(4) ~ SU(2) x SU(2), the 
partition function can be written as the product of an SU(2) Chern-Simons partition function 
with its complex conjugate 8 . Therefore, it is sufficient to consider the partition function (3.19) 
with the gauge group SU(2). In this case, "half-particles" at rest are now characterized by only 
one parameter, which corresponds to a choice of a unitary irreducible representation of SU(2). 
We will denote it by rrii 6 N/2 and abusively refer to it as the mass of the "half-particle". The 
"half-particle" will also be abusively referred to as particle. 

Pure gravity and quantization of the level 

Let us start with pure gravity. While studying the symmetries of the Chern-Simons action, 
we have mentioned the presence of an anomaly (3.10) in the large gauge transformations. This 
anomaly has no effect in the classical theory, but its consequences at the quantum level are 
important. Indeed, if one requires that the formal measure in the path integral (3.19) be 
gauge invariant under the action of the gauge group, then the following identity holds for any 



where the winding number w[M, g] introduced in (3.10) is non-trivial when we the gauge trans- 
formations are not connected to the identity. For this relation to make sense, w[A4,g] must be 

8 The fact that we have the complex conjugate is due to the expression of the non-degenerate bilinear symmetric 
form which defines the Spin(4) Chern-Simons action. Indeed, when decomposing the Chern-Simons action into 
self-dual and anti self-dual parts, the associated levels come with an opposite sign. 




(3.19) 



g G C°°(7W,SU(2)): 
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an integer modulo 2ir. As the group SU(2) is compact, for any g and every closed manifold A4, 
w[Ai,g] = 2irNk where N is an integer. Since N can take in particular the value 1, k must be 
an integer as well. 

Let us emphasize that the path integral quantization is the most suitable framework to see 
that the level k is discrete. We are not aware of any argument for the discreetness of k coming 
from the canonical quantization. This is an important input from the covariant quantization 
and, as we will see in the remaining of this work, it is a good illustration of the complementarity 
between the covariant and canonical approaches. 



Integrating out the particle degrees of freedom 

Now we turn to gravity coupled to point particles. In order to compute the partition func- 
tion (3.19), it is easier to first integrate out the degrees of freedom of the particles. The integrals 
over the variables gj decouple. They are in fact all identical and given by 

W^[A] = j Pgexp^y^H^-Mg + g- 1 ^, (3.20) 

where m £ N/2 is the representation of SU(2) corresponding to the mass of the particle, 7 is 
a closed loop in the manifold A4, and x( m ) £ su (2) is a given element in a Cartan subalgebra 
u(l) C 5u(2). As the full theory (gravity coupled to particles) is invariant under the action of 
the gauge group, it is natural to require that the measure f Vg in (3.20) be also invariant under 
the gauge transformations. They reduce in the present case to the right action of an element 
g £ C°°(7,SU(2)) on the group variable g. An immediate consequence of the requirement of 
invariance is that 

w\ m) [A 9 } = w\ m) [A]. (3.21) 

Therefore, as it was suggested by the notation, [A] is in fact the Wilson loop of the connec- 
tion A along the closed loop 7 evaluated in the representation of spin m. To be more precise, 
equation (3.21) implies that W^"^L4] is proportional to the Wilson loop [64, 150], but one can 
always choose a normalization of the measure j T>g in such a way that the proportionality factor 
is trivial. As a consequence, the partition function of gravity coupled to point particles can be ex- 
pressed as the expectation value of a product of Wilson loop observables in Chern-Simons theory: 



Z(M;M 1 ,...,M P ) = J VA {^[W^[A]^ exp (iS cs [A]) 



Complementarity of the covariant and canonical quantizations 

Now we have all the ingredients to see how the canonical and covariant quantizations are related 
to one another. In particular, we will show that a complete description of the physical Hilbert 
space inherited from the canonical quantization of the theory provides a machinery to compute 
explicitly the path integral. 

Following [150], let us decompose the spacetime manifold M. into two pieces A4i and Ai 2 
sharing the same boundary £ = dM± = dM.2- Furthermore, let us assume that the common 
boundary has the topology of a two-sphere § 2 . When particles are present in the theory, their 
worldlines cross this boundary, which can then be interpreted as a punctured two-sphere. 

In the previous subsections, we have seen that any punctured surface is endowed with a certain 
symplectic structure. Therefore, the canonical quantization of the theory on these punctured 
surfaces associates a Hilbert space H to § 2 . In the next subsection we provide a precise descrip- 
tion of this Hilbert space. For the moment, we are only going to recall some of its properties: 
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• When there are no punctures on § 2 , H is one-dimensional. 

• When there is one puncture on § 2 , 7i is zero-dimensional. This follows from the fact that 
on a closed manifold without boundaries the wordlines of particles are closed loops and 
the number of punctures of § 2 is necessarily even. It is therefore only with an even number 
of punctures that the dimension of the Hilbert space is non-zero. 

• When there are two punctures on § 2 , 7i is one-dimensional. 

• When there are more than two punctures, the dimension of the Hilbert space will depend 
non-trivially on the representations carried by the punctures, but will always remain finite- 
dimensional 9 . For instance, when there are four punctures colored with the spin m = 1/2, 
T~L is two-dimensional. 

The path integral is directly related to the Hilbert space T~L in the sense that it can be written 
as a scalar product 

Z(M;M 1 ,...,M p ) = (v 1 ,v 2 ) (3.22) 

between two vectors V\ and v 2 in More precisely, v\ is defined as an element of the dual 
space T~L* , which we have implicitly identified with %. Formally, v\ and v 2 correspond to the 
vectors given by the path integral (3.19) restricted to the manifolds Mi and M 2 . This is 
a standard property of any topological quantum field theory. 

Furthermore, let us assume that there are no punctures. Then, Ti is one-dimensional and 
the path integral Z{M) is still given by the scalar product as in (3.22). Let wi, w 2 be the 
corresponding vectors for M = § 3 , which are a priori different from v\ and v 2 . Then one can 
use the "marvelous" property [150] 

{v 1 ,v 2 ){wi,w 2 ) = {vi,w 2 ){w 1 ,v 2 ), 

which follows trivially from dim % = 1 . If one translates this identity in the language of partition 
functions, one immediately obtains the factorization property 

Z(M)Z&) = Z(M 1 )Z(M>) <=* Z(M)_Z(M l )Z(M 2 ) 



Z(8 3 ) Z(8 3 ) Z(S 3 ) ' 

where Mi is the connected sum of Mi with a three-ball, which is illustrated in Fig. 2. 

It is remarkable how simple arguments coming from both the canonical and the covariant 
quantizations allow to derive such a non-trivial property. With the same construction, we can 
extend this property to the case where particles Mi, . . . , M p are present. If one assumes that 
the loops 7j associated to the particles Mj are not knotted one to the other, then one has 

Z(M; Mi, . . . ,M P ) = fr Z(M;Mj) 
Z(§ 3 ) 11 Z(<B 3 ) ' 

Proving these factorization properties is certainly enough to illustrate the complementarity 
between the covariant and the canonical quantization schemes. The path integral gives the geo- 
metric viewpoint, whereas the physical Hilbert space brings the algebraic tools. However, we 
want to go further and show how to establish the relationship with the Jones polynomials. 

In order to do so, let us consider the situation depicted in Fig. 3. It shows three knots in §> 3 
which coincide in the "southern hemisphere" but differ in the "northern hemisphere". The 



9 There is a nice way to understand physically the finiteness of the Hilbert space [150]. It is easy to show 
that the classical phase space has (generically) a finite volume computed with the symplectic volume form when 
the gauge group is compact. As the Hilbert space is a discretization of the phase space into quantum cells of 
microvolume ti, it is natural to expect that it has a finite dimension. 
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Figure 3. Examples of knots in S 3 . They coincide in the "southern hemisphere" of § 3 but differ in the 
"northern hemisphere" . 



shape of the knots in the southern hemisphere is not important for the moment. According to 
the previous discussion, the knots in the northern hemisphere of § 3 define three vectors ui, V2 
and v% in the Hilbert space H associated to the sphere with four punctures. If we assume for 
simplicity that the knots are all colored with the 1/2 representation, then dim% = 2. As a 
consequence, v\, V2 and V3 are not independent, and there are three real numbers a 1 , i = 1,2, 3, 
such that a l Vi = 0. The portion of the knots lying in the southern hemisphere defines a vector 
w in the (dual of the) Hilbert space H such that 

(vi,w) = Z(§ 3 ,M i ), 

where Mi denotes the particle whose worldline is the knot obtained by gluing the northern part 
of S 3 associated to v\ with the southern part associated to w. As a consequence, we obtain the 
following relation: 

a 1 ^(S 3 , Mi) + a 2 Z(<B 3 , M 2 ) + a 3 ^(S 3 , M 3 ) = 0. 

This relation is important because it is closely related to the skein relation defining the Jones 
polynomial. To see that this is indeed the case, we need to compute the explicit values of the 
coefficients a 1 , and in order to do so the concrete description of the physical Hilbert space T-L 
which will be developed in the next subsection is required. Using this description, we find that 

qZ(S 3 ,M 1 ) + (q 1 / 2 -q- 1 / 2 )Z(S 3 ,M 2 )-q- 1 Z(S 3 ,M 3 ) =0, (3.23) 

where q = exp(2'm/(k + 2)) is a root of unity because the level k is an integer. Before going 
further, let us point out that the values of a 1 depend on the choice of a framing corresponding to 
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K + K- K Q 

Figure 4. The three "crossings" appearing in the definition of the skein relation. 

a thickening of the knot into a ribbon. There is no canonical choice of framing, but two different 
choices lead to scalar products which differ only by a phase. In our example, a canonical framing 
exists (it is the trivial one), and with this choice (3.23) is exactly the skein relation 

qJ(q,K + ) + (q 1 / 2 -q- 1 / 2 )j(q,K )-q- 1 J(q,K-) = 0, 

which defines the Jones polynomials J(q,K) for any knot K embedded in S 3 , up to a global 
normalization [98]. In this relation, J(q,K) is a Laurent polynomial in the variable q, and 
K + , ET_ and Kq are three different knots which coincide with a knot K everywhere but in 
a small region around a crossing, where they differ according to the representation of Fig. 4. 
In order to evaluate the Jones polynomial of a knot, one needs to give it an orientation, but in 
the end the evaluation does not depend on this choice. As a consequence, the path integral of 
three-dimensional Riemannian gravity with a positive cosmological constant coupled to a point 
particle is closely related to the Jones polynomial of the knot associated to the worldline of 
the particle. One can make this statement more precise, but this would require an extended 
discussion about the role of the framing, and we want to avoid this for the sake of simplicity. 

This construction can be generalized to compute the path integral of the Chern-Simons 
theory with compact gauge group in any closed three-dimensional spacetime using the surgery 
representation of three-manifolds [150]. When the gauge group is SU(2) and the knot is colored 
with a spin-m representation, the path integral gives the generalized Jones polynomial. When 
the gauge group is SU(iV), it can be related to the so-called HOMFLY polynomial, which is 
a two-parameter generalization of the Jones polynomial. 

3.2.2 The combinatorial quantization — emergence of quantum groups 

This part is devoted to the construction of the Hilbert space % associated to a punctured 
Riemann surface S. To our knowledge, the most powerful technique to perform the canonical 
quantization of Chern-Simons theory is the combinatorial quantization. It is based on the 
seminal work of Fock and Rosly [70], and has been used in [2, 3, 4] to quantize Chern-Simons 
theory with the gauge group G = SU(2), in [50] when G is the Lorentz group and in [113, 114, 115, 
116] when G is the Euclidean group. These three cases correspond respectively to Riemannian 
gravity with a positive, negative, or vanishing cosmological constant. The Lorentzian case with 
A = has been studied as well, but we restrict ourselves only to the Riemannian case. 

The combinatorial quantization relies on the fact that the quantization turns classical sym- 
metry groups into quantum groups. In this approach, quantum groups are not introduced 
artificially by hand, but appear naturally in a very beautiful way during the quantization pro- 
cess. Let us start by presenting the quantum groups which are at the heart of this combinatorial 
quantization. The classical (universal coverings of the) symmetry groups of Riemannian gravity 
are the Euclidean group ISU(2), the Lorentz group SL(2,C), or the rotational group Spin(4) de- 
pending on the value of the scalar curvature. These groups are turned into quantum groups that 
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can all be described as Drinfeld doubles [53]: DSU(2), DU g (su(2)) for q £ R, and DU g (su(2)) 
with q a root of unity, respectively, where U 9 (g) denotes a deformation of the associated en- 
veloping algebra U(g) of a Lie algebra g. In fact, DUq(su(2)) for q real or root of unity can 
be shown to be a deformation of either SL(2,C) or Spin(4), and therefore these groups are also 
denoted by U g (sl(2, C)) and U g (so(4)). The deformation parameter q is the same as the one 
which appears in the skein relation (3.23) and is related to the level k of the Chern-Simons 
theory. When q is a root of unity, one has q = exp(2i7r/(/c + 2)), whereas for q real the relation 
is known only at first order in 1/k and reads q = 1 + 2ir/k + o(l/A;). Thus, the quantization 
leads to two important physical effects: the level k becomes an integer and classical groups are 
turned into quantum groups. As a result, the mass and the spin of the particles M, are, after 
the quantization process, given by unitary irreducible representations of the quantum groups 
instead of classical groups. In this sense, the mass and the spin become quantized. 

In order to make all the previous statements more precise, let us assume for simplicity that X 
is a sphere with p punctures colored with SU(2) representations rrii for % = 1, . . . ,p. The starting 
point of the combinatorial quantization is the definition (3.18) of the physical phase space, which 
reduces to 



a G O mi ® •••<E> O mp | l\a t = 1 /G, (3.24) 

i=i > 

where O m C G is a conjugacy class in SU(2) characterized by the spin m. Usually, conjugacy 
classes are labeled by an angle 9 G [0, tt], but in the quantum theory the angle is discrete and 
related to m by 

6 m = 2irm/k, (3.25) 

where k is the level. This relation suggests that m < k/2, which has an immediate interpretation 
in the quantum theory as we will see very soon. If we denote by t\ one generator of su(2), then 
a representative element of O m is given by h(9 m ) = exp(# m Ti) and the conjugacy class is 
explicitly defined by 

O m = {h€ SU(2) | h = gh(Q m )g- x , Vg G SU(2)} c SU(2). (3.26) 

It is clear from this definition that O m ~ S 2 and that any of its elements h satisfy XjQ 1 ) = 
Xj{h{0m))i where XjQ 1 ) is the character of h evaluated in the representation j. Thus, the 
space (3.24) is a finite-dimensional symplectic manifold whose Poisson bracket is inherited from 
the canonical analysis of Chern-Simons theory. The equation 

v 

F = Y[ ai = l (3.27) 

i=l 

expresses the flatness of the connection, and G acts by simultaneous conjugation, i.e. a% i — > 
g~ x a>ig for any i G [l,p] and any g G SU(2). Furthermore, F is the Hamiltonian generator of 
this symmetry. 

The combinatorial quantization proceeds by quantizing the space of connections of the punc- 
tured surface S, 

A(E) = O mi x • • • x O mp 

and then by imposing the flatness condition (3.27) at the quantum level. Any element aj in the 
space of connections can be viewed as the holonomy of the Chern-Simons connection starting 
at a given base point and turning around the puncture i. The problem is that such a space does 
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not possess a regular Poisson structure, and the one inherited from the canonical analysis is ill- 
defined due to the following reason: Although one knows how to compute the Poisson bracket 
between two holonomies associated to paths which cross each other, the Poisson bracket between 
holonomies associated to paths which end at the same point but do not cross is not defined. 
This is precisely our case since the paths associated to the elements of A(T,) do not cross, they 
just end or start at the same base point. As a consequence, «4(S) is not naturally endowed with 
a Poisson bracket and one needs to introduce a certain regularization. This is exactly what Fock 
and Rosly did in a very clever way [70]. They have shown that if one endows the space 
with a linear ordering (which amounts to saying that a\ comes first, a 2 second, and so on ... ), 
then it is possible to endow it with a quadratic Poisson-Lie structure schematically given by 
a formula of the type 

{ai 1, 1 <g> o 2 }fr = -r{ra\ <£> a 2 + a\ ® a 2 V), 
k 

where r is called a classical r-matrix (V is its transpose) and satisfies the so-called classical 
Yang-Baxter equation as a consequence of the Jacobi identity. In (3.28), a\ and a 2 are viewed 
as group elements 10 (in fact, they are holonomies), and r is an element of (8> 0. Similar 
Poisson brackets hold for any couples (aj, aj) of A(E). This prescription possesses the following 
important features: 

• First of all, it reduces to the standard Poisson bracket on the physical phase space (3.24) 
after imposing the flatness condition, or equivalently the gauge invariance. 

• Second, its quantization is well-known and leads to the appearance of quantum groups. In 
particular, the r-matrix becomes the classical limit of a quantum i?-matrix which satisfies 
in turn the quantum Yang-Baxter equation. This is a clear sign of the quantum group 
structure behind this construction. 

Indeed, when one proceeds to the quantization of the (non-physical) phase space A(E) ac- 
cording to the standard rules of quantum mechanics, the resulting quantum algebra, called the 
graph algebra £(£), is generated by the quantization hi of the classical variables a% and is defined 
by commutation relations of the form 11 

R{a x (g> (g> a 2 ) = (a 2 1)12(1 ® a^lT 1 , 

where R is the quantum .R-matrix. As we already said, the presence of R is a strong indication 
of the presence of quantum groups. This is in fact what Alekseev has shown in a beautiful 
theorem [1] stating that 

£(£) ~ U ff (*u(2))«* 

i.e. that £(£) is isomorphic to p copies of the quantum group U g (su(2)). This theorem (which 
can be generalized to any Riemann surface) allows to construct the kinematical states of the 
theory. They are given by unitary irreducible representations of £(£) which are obtained from 
the tensor product of p unitary irreducible representations ji of U g (su(2)), 

p 

H (j 1 ,...,j p ) = (g>H( ji \ 

8=1 

10 To make the Poisson-Lie structure more concrete, it is useful to evaluate the holonomies a; in finite- 
dimensional representations of SU(2) so that the variables of the phase space become the matrix elements 
of 7r ij) (ai). 

"More precisely, the graph algebra is generated by the quantization of the matrix elements 7r^(aj) which can 
be denoted 7P 3 '(aj), but with j labelling a representation of the quantum group instead of the classical group. 



(3.28) 



(3.29) 
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where H^) is the module of the representation labeled by j G [0, k/2]. The cut-off on the 
number of representations is a result of the quantization and is totally consistent with the 
classical relation (3.25) that we have obtained above. As in the classical case, the dimension of 
the representation of the quantum group is dj = 2j + 1 . 

Note that Hq is not the Hilbert space of physical states. Rather, it is analogous to the 
kinematical Hilbert space in loop quantum gravity. For this reason, we have called its elements 
kinematical states. To extract the physical Hilbert space out of T-Lq, one still has to impose 
two types of constraints. The first one is the quantum analogue of ai G O mi and fixes the 
representations 12 in (3.29) to j% = mi. Therefore, among all the Hilbert spaces Hq, only one of 
them is admissible. The second constraint is the quantum flatness condition F = 0. To impose 
this constraint, one has to understand how the action of the symmetry group G translates at 
the quantum level. As we know, the classical action of G is turned into the action of U 9 (su(2)) 
on Ho viewed as a tensor product of U q (su(2)) modules. As a consequence, the physical Hilbert 
space is simply given by 

U = lnv(U {mi) n {mp) ) , 

where the notation Inv denotes the subspace invariant under the action of U 9 (su(2)). 

Due to the U g (su(2)) recoupling theory, it is clear that dim'H is equal to one when there are 
no punctures, and to one or zero when there are three punctures at most. When there are more 
than four punctures, the dimension becomes a non-trivial function of the spins rrij. In the case 
in which there are four punctures colored with the spin 1/2 representation, one has 

(^(i/2))®4 = 2n (o )(B m n) e n m 

and dim% = 2 as was claimed previously. In order to justify the results that we have used 
to compute the skein relation (3.23), we now give a more precise description of T~L and show 
how the vectors in T~L are related to the path integral on the sphere with boundaries. These 
vectors and some of their properties can be given a nice pictorial representation as in Fig. 5. 
The interpretation of this graphical representation is the same as in the classical case, the only 
important difference being that now there is a non-trivial braiding. Using this representation, it 
is immediate to relate the three vectors V{ defined by the path integral in Fig. 3 to the elements 
of %. This relation is given in Fig. 6. The vectors v% are viewed as elements of T~L in an obvious 
way. Indeed, one identifies the projection of their graph onto a two-dimensional surface with 
a quantum spin network. This projection is not canonical and depends on the choice of the 
framing, but the "physical results" are independent of these choices. Once this identification is 
made, the calculation of the coefficients a 1 (up to a global factor) is straightforward and one 
ends up with the skein relation. 

This example is a good illustration of the interplay between the covariant and canonical 
quantizations of Chern-Simons theory with a compact gauge group. This construction can be 

12 It is possible to show that for any representation j, the elements 

W Q) =Tr w)(a i ) = Tr w)(^. )j 

where the quantum trace Tr q is defined from the ribbon element fi £ U q (su(2)), are in the center of the graph 
algebra £(£)• Therefore, these elements form an Abelian subalgebra and they are scalar operators in Ho- The 
quantum version of the constraint £ O mi fixes the value of each element Wf to 

where the quantum number [x] q is defined in (3.31). In fact, is the quantum analogue of the character Xj( m i), 
and it is easy to show that the quantum constraint fixes ji = m,\ 
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basis of physical states 




and 



quantum trace 



^ = Tr 9 (l) = 9 V2 + r 1/2 



twist 



= 



Figure 5. Graphical representation of the physical states on the sphere with four punctures labeled 
with j — 1/2. On the top, we have given two different basis. The one on the left is described in terms of 
intertwiners (j € {0, 1}), while the one on the right has two different vectors (one of them being defined 
using the quantum R- matrix). In the middle, we have represented the quantum trace. At the bottom, 
we have given some properties related to the framing. The links have to be understood as ribbons, and 
therefore the twist has a non-trivial effect on the evaluation of the quantum scalar product. The phase 
is given by = — q 



-3/2 



path integral representation 



vectors in 7i 







Figure 6. A map from the path integral representation of the physical states to the elements of the 
physical Hilbcrt space R,. 

generalized to any arbitrary manifold Ai using its surgery representation. On the other hand, 
the generalization to non-compact gauge groups is still an open problem [149]. 



3.2.3 The state sum quantization — spin foams 



The previous quantization, a la Witten-Reshetikhin-Turaev, relies technically on the canonical 
quantization. Even if the use of the path integral was necessary to establish the relation with 
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knot invariants, physical amplitudes were explicitly computed using the scalar product on the 
physical Hilbert space. Therefore, this quantization appears to be more canonical than path 
integral-like. 

After the seminal paper of Witten, the path integral quantization of Chern-Simons theory or 
three-dimensional gravity has been extensively studied. The perturbative quantization of Chern- 
Simons theory was initiated soon after the work of Witten and led to a better understanding 
of the Vassiliev invariants. The non-perturbative quantization was also studied and led to the 
first mathematically well-defined spin foam model, the so-called Turaev-Viro model, which we 
will now briefly review. 

The original idea was to construct some manifold and knot invariants in three-dimensional 
spacetime in a totally covariant way, by taking into account an important indication coming from 
the canonical quantization, namely the fact that quantum groups play a fundamental role. The 
construction is as follows. Let M be a closed three-dimensional spacetime, A a triangulation 
of M, and A* its dual triangulation. Elements of A* consists of a set of vertices v (dual to 
tetrahedra r in A), edges e (dual to triangles t in A) and faces / (dual to links I in A). To each 
face / 6 A* (or equivalently to each link t of the triangulation), we assign a unitary irreducible 
representation of U g (su(2)) labeled by a spin jj. Then one can define the following quantity 

- {-wrkT E IT ra, + i],n{ % %%}. (3-30) 

where no denotes the number of vertices of A, the sum is taken over all possible assignments 
of irreducible representations of U g (su(2)) to the faces, and the two products are taken over all 
the faces and vertices of the triangulation, respectively. The labels j\ , . . . , j% are associated to 
the six links I bounding a tetrahedron r £ A. Finally the weights associated to the faces are 
the quantum dimensions 

[x] q = ^El!, (3.31) 
q-q 1 

with x = 2jf + 1, and the weights associated to the vertices are the totally symmetrized quan- 
tum {6j} symbols. The quantity defined by (3.30) is finite and was shown to be independent 
on the choice of the triangulation. Therefore, it gives rise to a well-defined topological invariant 
of the manifold Ai and we can write that 



Z TY (A*) = Z TY {M). 

It was realized later on, after the introduction of the Ponzano-Regge model and the study of 
its asymptotics, that the Turaev-Viro model is in fact a quantization of Riemannian three- 
dimensional gravity with a positive cosmological constant. As the action for gravity is the 
product of two uncoupled Chern-Simons actions with opposite levels k and —k, one can formally 
write that 

Z Ty (M) = \Z WRT (M)\ 2 , 

where ^vvrt(A'I) is the canonical evaluation a la Witten-Reshetikhin-Turaev of the path in- 
tegral, as described in the previous subsection. This relation has been rigorously proved very 
recently [145], and it establishes a precise link between the Turaev-Viro model and the physi- 
cal scalar product of the canonical quantization. Even if this relation seems natural from the 
physical point of view, the rigorous mathematical proof is much more involved. 

The construction of the covariant path integral can be generalized to accommodate for the 
presence of particles [36, 83, 84]. Let us recall that particles in Riemannian gravity with a positive 
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cosmological constant are classified by the unitary irreducible representations of so(4) and are 
therefore characterized by two half-integers. When the particles are assumed to be without spin, 
they are classified by their masses only (which label the so-called simple representations, see 
Section 4.3.1) and are characterized by one half-integer denoted m. After coupling to gravity, 
m becomes a label for the (simple) representation of the quantum group U 9 (so(4)) and as such 
it is bounded by k/2. 

In order to include massive particles in the framework developed by Turaev and Viro, we 
consider once again a triangulation A of Ai in which particles live on the links (one-cells) I G A 
and have trajectories describing closed loops 7, in A. The generalized Turaev- Viro amplitude, 
which now depends on the data Mj (the mass mi of the i th particle and the closed loop 7$), is 
given by 



where A^ denotes the set of faces dual to the links I where the particles live, and A^ de- 
notes the set of faces dual to the links I where there are no particles. It has been shown that 
-Ztv(A; Mi, . . . , M p ) is a knot invariant which does not depend on the choice of the triangula- 
tion A for M, and also that 



This provides a clear relation between the covariant and canonical quantizations and demon- 
strates their mutual consistency in the present case. 

3.2.4 Discussion 

Before going to the loop and spin foam quantizations of three-dimensional gravity, we would like 
to end the review of the quantization of Chern-Simons theory with a short discussion. First, 
let us emphasize once again the important role played by quantum groups in the canonical and 
covariant quantizations. They appear in a very natural way in the combinatorial quantization 
and were used to define the Turaev- Viro state sum. The question whether they are fundamental 
symmetries of quantum gravity or simply mathematical tools to construct the physical Hilbert 
space is still open. Indeed, quantum groups can be viewed as an intermediate step towards the 
construction of physical states, which are at the end of the day invariant under the action of 
quantum groups. In this sense, we do not see their effect on the physical states. 

Here we have considered exclusively Riemannian gravity with a positive cosmological con- 
stant. The reason is that it is the only case relevant for gravity which involves a compact group. 
All the other cases (Riemannian or Lorentzian) are formulated using non-compact gauge groups 
and they are much more difficult to deal with. It is nonetheless possible to perform the canonical 
quantizations for Riemannian gravity with a vanishing or negative cosmological constant. How- 
ever, the spin foam quantizations of these cases lead to state sums that need to be regularized. 
Usually, the regularization breaks the topological invariance and this makes these models less 
interesting. 

Thus, we conclude by emphasizing that Riemannian three-dimensional gravity with a positive 
cosmological constant admits a complete quantization, both from the canonical and the covariant 
point of view. It is therefore a good ground to test the ideas and the techniques of loop quantum 
gravity. In particular, LQG should reproduce the results obtained by Witten to claim that it is 
a serious candidate for quantizing gravity. We will see in the next section if this is indeed the case. 





Z TY (A*;M 1 ,...,M p ) = \Z WRT {M;M 1 ,...,M p )\ 2 . 
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3.3 Quantum gravity in the BF formulation 

Loop quantum gravity has been introduced as an attempt to quantize gravity in four dimensions. 
Soon after the foundations of LQG were established, it was applied to the quantization of three- 
dimensional gravity in order to test its ideas and techniques [20, 24, 80, 110, 121, 143]. Although 
many aspects of three-dimensional quantum gravity can be recovered from the LQG point of 
view, the results presented in the previous subsection are still out of reach. For instance, in 
the Riemannian case without cosmological constant it is possible to perform the construction 
of the physical Hilbert space [121] and to see that the Drinfeld double of SU(2) (the relevant 
quantum group in the combinatorial quantization scheme) plays a central role, even if it is 
somehow hidden [118]. On the other hand, the construction of the physical Hilbert space in the 
presence of a cosmological constant is still an open problem in LQG. As a consequence, in the 
context of LQG, one does not see the emergence of quantum groups as deformations of classical 
gauge symmetry algebras, and the link with knot invariants and the Jones polynomial remains 
a mystery. Note however that the preliminary results of [122] strongly suggest that the Jones 
polynomials could appear in LQG, but many things remain to be done in order to make this 
statement more precise. 

The story with spin foam models is a bit different. In contrast to LQG, they were first 
introduced in three dimensions, and then generalized to four dimensions. The first spin foam 
model was proposed by Ponzano and Regge [132] as an attempt to compute non-perturbatively 
the path integral of Riemannian gravity with vanishing cosmological constant. However, the 
original model suffers from divergences and requires a suitable regularization. Although it is 
possible to do this without modifying the original theory, a simple and nice regularization is 
provided by the inclusion of the cosmological constant and leads to the Turaev-Viro state sum 
model [146]. The cosmological constant compactifies the symmetry gauge group already in the 
classical theory, which has the effect of curing the divergences 13 . Thus, in three dimensions we 
have at least two well-defined spin foam models. 

This section is devoted to reviewing some aspects of loop quantum gravity and spin foam mod- 
els in three dimensions. We would like to focus on the main results and the main open questions 
arising in this context. Both quantizations are based on the BF formulation of three-dimensional 
gravity instead of the Chern-Simons one. We start with a brief overview of the Ponzano-Regge 
model. In the second part, we review the loop quantization approach. Then, we make a contact 
between the canonical and covariant quantizations. Finally, we explain the preliminary results 
obtained in the LQG framework for the case of a positive cosmological constant. 



3.3.1 Ponzano Regge model 

The idea leading to the spin foam quantization of gravity is to use a discretization of the 
spacetime manifold in order to define a regularization of the formal partition function 

Z 3d = J VeVuj exp ^i J tr (e A F[u])j . (3.32) 

To this end, we choose an arbitrary cellular decomposition A of the three-manifold M, together 
with its dual two-complex A*. Since the triad e l a is a one- form, it is natural to integrate it along 
the one-cells £ of A (dual to / £ A*) to form the su(2)-valued elements Xt. The connection 
of the gauge group G = SU(2) is discretized by computing its holonomy g e along the edges e 
of A*, so that the discrete curvature is associated to the faces and is given by the product of 



1 Note that this procedure has been used as an attempt to regularize four-dimensional spin foam models as 
well [69, 93, 123]. 
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holonomies along the edges lying on the boundary of /: 

9f = ll9e 

eCf 

With these variables, the discretization of the action appearing in (3.32) can be written as 

S[X f ,g f ]= £ **(Xf9f), 
/eA* 

and the partition function becomes 

Z PR (A*) = ( II / dX f I f II / d 9e ] exp [ i £ K(Xf9f) ) , (3-33) 

where dXf is the Lebesgue measure on 5u(2) ~ M 3 and dg e the Haar measure on SU(2). The 
integral over Xf can be trivially performed to find 

Z PR (A*) = ( IT / d 9e) IT <%)> (3.34) 

where the delta distribution over SU(2) imposes the flatness of the connection 14 . Using the 
Peter- Weyl decomposition, one can write 

*p*(a*) = e n (aw + 1) ( n L m d ^ ) n n4 (3.35) 

where the sum is taken over all possible SU(2) representations j labeling the set of faces / G A*. 
For arbitrary cellular decompositions A, let us call F the number of faces meeting at every edge 
e G A*. In (3.35) we will have an integral over g e of F products of its representation matrices. 
Such an integral can be evaluated as follows: 

f dg e D&>Gfe) • • • D^)( 5e ) = £ 

where J*= (ji, • • • ,Jf) denotes the collection of spins coloring the faces meeting at the edge e, 
D^(fl'f) is the spin-j representation matrix of the group element g#, and the sum runs over an 
orthonormal basis in the space InVgy^ = Inv^g^Lj ® • • • <g> 'Hgu( 2 )) °^ invariant intertwiners 
between the representations j. Then, all the intertwiners meeting at a vertex v can be contracted 
to define a vertex amplitude A v (j f Dv ,i eDv ). Finally, the partition function can be written as 
a sum over colorings decorating the dual two-complex of the amplitudes constructed from A v : 

Z PR (A*) = EE II ("i) 2 ^/ + 1) II Mjf,ie). (3.36) 
j-+f i-s>e /eA* u£A* 

To clarify the meaning of this formula, let us assume that the cellular decomposition A is 
simplicial (i.e. three-dimensional). In this case, vertices in A* are four-valent, while edges are 
three- valent. The vertex amplitude is therefore given by a contraction of four (normalized) 



14 Technically, the integration over the variables Xf leads to a delta function over SO(3) and not SU(2), but 
we will ignore this minor complication in the following. In fact, the model can be slightly modified in order to 



immediately yield the SU(2) delta function [83, 107]. 
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three-valent intertwiners, which is the so-called totally symmetric {Gj} symbol, and (3.36) takes 
the form [132] 

z PR (a*) = e n (-i) 2 ^/ + 1) n { £ 1 1 } ' (3 - 37) 

where the labels j\ , . . . , j§ are associated to the six links bounding the tetrahedron dual to the 
vertex v G A*, and we absorbed all possible sign factors associated to the edges into the normali- 
zation of the {6j} symbol (see [36]). Notice that the sum over intertwiners has now disappeared 
since there is a unique (up to normalization) three-valent intertwiner, i.e. dim (l nv gu(2)) = ^ 
Some important comments concerning the result (3.37) are in order: 

• An extremely important observation is that the partition function (3.37) does not depend 
on the choice of discretization A, which is a manifestation of the topological nature of the 
quantum theory. More precisely, given two different cellular decompositions A and A and 
using the Biedenharn-Elliot relation, one can prove the following relation: 

S - n °Z PR (A*) = s-"°Z PR (A*), (3.38) 

where s = ^2j(2j + 1) is a divergent factor and no is the number of zero-cells in A, i.e. the 
number of bubbles of A*. Thus, (3.37) defines a formal (due to divergences) three-manifold 
invariant and, as a result, one can formally write that 

Z PR (A*) = Z PR (M). 

In particular, this implies that an appropriate choice of a single two-complex is sufficient 
to provide the full transition amplitude and the sum over triangulations should be omitted 
in its definition. 

• Although the sum over different two-complexes in the Ponzano-Regge model is redundant, 
it is meaningful to ask the question of the sum over topologies. This question has been 
partially studied in [82] and [109] using the group field theory formalism for spin foam 
models. It was shown that the sum over triangulations in three-dimensional quantum 
gravity can be tamed rigorously and paves the way for the renormalization program in 
group field theory [74]. 

• The divergence of the sum over spins in Z PR (A*) and in the definition of the factor s 
makes the relation (3.38) rather formal. In fact, this divergence has a physical origin and 
arises due to a remnant of the gauge symmetry at the discrete level [83]. Indeed, the spin 
foam state sum is a version of the path integral where no gauge-fixing was done. On the 
other hand, the classical gauge group in the present case is the non-compact (Euclidean) 
Lie group ISU(2) and the divergences appear as a manifestation of its infinite volume. 
They can be removed by partially fixing the gauge freedom remaining after discretization, 
which amounts to eliminating the bubbles from A*. Another way to understand these 
divergences is to see that they result from the ill-defined product of distributions. This 
product can be regularized by replacing in (3.34) the ^-distributions over the group SU(2) 
by the heat kernel. This allows to explicitly identify the divergent part of the amplitudes 
and to relate it to the bubbles in A* [46, 47, 48]. Then the physical convergent part turns 
out to produce the well-known Reidemeister torsion. In fact, this was already expected 
from the work of Witten [149], and was proved in a less general case in [37]. As we are 
going to see in what follows, the canonical quantization enables one to avoid the infrared 
divergencies as well. 
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Finally, let us briefly discuss the inclusion of particles into the model. In the language of 
the quantization a la Witten, the particles correspond to Wilson loop observables in the Chern- 
Simons theory with the non-compact gauge group ISU(2). They are supposed to have closed 
trajectories describing knots 7« in spacetime. For the particles with masses m% and with no spin, 
the partition function of the coupled system is formally given by [83, 84] 

Z PR (A*;M 1 ,...,M P ) = Y J II (-1)^(2.7/ + 1) 

x II (-l) 2jf Xj f (m f ) JJ | jl 11), (3.39) 

where Mi = (mj,7j) is the data of each particle, A^ are the faces in the cellular decomposition 
that are crossed by the trajectories ji, and A*^ the faces that are not crossed. A face / 6 A*, 
initially colored by a representation jt and crossed by a particle of mass m, inherits a new 
color rrif which is precisely the mass of the particle mj = m. As a result, its weight Wf in the 
state sum is no longer (2jf + 1) but 

sin((27f + l)mf) 

Wf = X Jf (m f ) = — ^ '—11. 3.40 

This factor has a clear group-theoretic interpretation in terms of characters. Although the state 
sum (3.39) is still formal, in some cases it can be regularized. 

Note that the inclusion of particles makes it clear that Riemannian quantum gravity without 
cosmological constant is related to the quantum double DSU(2), since the particles are shown 
to be classified by unitary irreducible representations of DSU(2) and not by those of ISU(2). 
This fact can be illustrated in two ways. The first one is that the masses of the particles are 
bounded by 2ir, as can be seen from the expression (3.40), where the mass is defined mod 2ir 
only. The physical origin of this bound can be traced back to the conical singularity created 
by the particle after its coupling to gravity. The deficit angle of the cone coincides (in some 
units) with the mass of the particle, and therefore the mass cannot exceed the value 2ir. And 
we know that the unitary irreducible representations of DSU(2) are labeled by a couple (m, s) 
with m < 2ir, whereas this is not the case for representations of ISU(2). The second argument 
is based on the way in which interactions between particles are described, since it can be shown 
to follow the rules of DSU(2) recoupling theory [101]. This latter is drastically different from 
the recoupling theory of ISU(2). 

The coupling to particles offers the possibility to build a quantum field theory over a quantum 
background, which is often referred to as a third quantization. Let us explain how this comes 
about. It is not possible to fully quantize a matter (or gauge) field coupled to gravity using 
the standard techniques. Indeed, the standard approach consists in first quantizing the field on 
a flat background, and then quantizing perturbatively the gravitational field around the flat 
metric. We know that this does not work because of the non-renormalizability of perturbative 
gravity. Another idea would be to quantize the gravitational field first, while keeping the matter 
(or gauge) degrees of freedom classical. In four-dimensional gravity, such an idea seems to be 
unrealistic, because there is no regime in which spacetime is quantum-mechanical while matter 
remains classical. In three dimensions however, this idea does make sense because matter does 
not really interact with gravity. This was realized explicitly and first done in the context of 
spin foam models [77, 78], for matter fields that are massive point particles. The result is very 
interesting. Integrating out the gravitational degrees of freedom while keeping the particles 
classical turns the classical flat geometry into a non-commutative geometry. In other words, 
quantizing particles coupled to gravity, is equivalent to quantizing particles without gravity 
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but leaving in a non-commutative spacetime. Afterwards, it was realized [99, 118, 119] that 
the emergence of a non-commutative space has a very simple physical interpretation. Indeed, 
functions on this non-commutative space form a non-commutative algebra which is covariant 
under the action of the quantum double DSU(2). More precisely, the classical isometry algebra 
ISU(2) of the flat Euclidean space E 3 is deformed at the Planck scale and becomes the Drinfeld 
double DSU(2). This result is very exciting because it is an explicit realization of the old idea 
that quantum gravity turns classical geometry into quantum geometry. It can be generalized 
to the case of Riemannian gravity with a non-vanishing cosmological constant. When A > 0, 
the classical isometry algebra so (4) is quantized into DUq(su(2)) with q a root of unity, and 
when A < 0, the classical isometry algebra so(3, 1) is deformed into DU 9 (5u(2)) with q real. 
As expected, these quantum groups are the ones appearing in the context of the combinatorial 
quantization. The associated non-commutative spaces are the quantum analogues of the three- 
sphere § 3 (for positive A) and of the three- hyperbolic space H 3 (for negative A). 

3.3.2 Canonical loop quantization with A = 

Now let us look at the canonical loop quantization of three-dimensional Riemannian gravity. 
The starting point is the classical phase space defined by the Poisson brackets (3.15) together 
with the constraints (3.16). We do not consider particles for the moment and we assume that 
there is no cosmological constant, i.e. A = 0. In contrast to the combinatorial quantization, we 
start with an infinite-dimensional phase space and take its quotient with an infinite-dimensional 
gauge group generated by the infinite set of constraints (3.16). The result is a finite-dimensional 
phase space, as follows from the classical analysis. This strategy is chosen in order to mimic the 
situation in four dimensions. In this way, one can use three-dimensional gravity as a toy model 
to test the ideas of the loop quantization. 

As in the case of the combinatorial quantization, the canonical quantization proceeds in two 
steps. The first one is to promote the classical variables to an operator algebra and to determine 
its unitary irreducible representations, which defines the space of quantum states. In loop 
quantum gravity, the state space is obtained by completion of the so-called space of cylindrical 
functions, which itself is obtained by taking the projective limit of the spaces Cyl(r) = 1? (G Lr ) 
associated with particular graphs T. Cyl(r) is the space of functions of the G- valued holonomies 
assigned to the links t of the graph, with the scalar product defined by 



where f\ and fi represent the functions on the group and Lr denotes the cardinal of the set of 
oriented links of T. Note that the topological nature of the theory in three dimensions allows one 
to work on a single graph, as long as it is sufficiently refined to resolve the topology of S. The 
resulting quantum theory will be independent of the choice of the graph. Therefore, in three 
dimensions, the projective limit is redundant. The second step is to construct the kinematical 
and physical Hilbert spaces. This can be done by promoting the constraints (3.16) to operators 
acting on the space of cylindrical functions I? (G Lr ). Schematically, kinematical states span 
the kernel of the quantum operator associated to the discretized version of the torsion T(x), and 
physical states are kinematical states that are in the kernel of the operator corresponding to the 
curvature F(x). 

The kinematical Hilbert space 

The basic discrete variables of loop quantum gravity are the holonomies ge of the connection u 
along links t of T and the "fluxes" Ei of the triad field along the links £* dual to the links I. 




(3.41) 
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They are defined by 




where exfj denotes the path ordered exponential 15 . On the space Cyl(r) the holonomy acts by 
multiplication and the flux by derivation. These actions are completely characterized by the 
relations 



9t> > 9i = 9e>ge, E\, > 9e = 4>i 9e 1 T t ge 2 , (3.42) 

where £p\ = KG is the three-dimensional Planck length, £% is the part of the oriented link I before 
the intersection I n (£')*, and £2 is the remaining part. This formula is a direct quantization of 
the Poisson bracket (3.15). It follows that the operators Eg act as vector fields on the space 
of cylindrical functions and the action (3.42) provides a unitary representation of the algebra 
of quantum operators on Cyl(r). Note that the flux operator is not gauge covariant as it is 
written, however its action on a holonomy is gauge covariant and therefore it is well-defined. 

The kinematical Hilbert space T^kin(r) is obtained as the set of solutions to the quantum 
Gauss constraint. They appear as elements of I? (G Lr ) that are invariant under the action of 
the gauge group G at the nodes n of the graph: 

H kin (T) = L 2 (G L r)/G N r. 

Due to the left and right invariance of the Haar measure on G, the norm (3.41) is compatible 
with the action of the gauge group and therefore induces a norm on %kin(r)j whereas the 
representation defined in (3.42) provides a unitary representation of kinematical operators on 
this space. 

In our case, G = SU(2), and a dense basis of T-Ly m (T) is provided by the spin network 
functions, obtained by assigning a representation of G to each link £ and an intertwiner to each 
node n. More precisely, harmonic analysis over the group SU(2) enables one to expand any 
cylindrical function in terms of unitary irreducible representations as 

With this decomposition of the states, it is easy to see that imposing the Gauss constraint 
amounts to requiring that the intertwiners, attached to the nodes of the graph and used to 
contract the indices of the representation matrices, be SU(2)-invariant: 




An orthonormal basis of the gauge- invariant Hilbert space %ki n (r) is then given by the so-called 
spin network states 

S { Tjd9i\ ^®V^{ 91 ) • ® »(3 (3.43) 



15 Notice that the fluxes defined in this way are not gauge-covariant under the action of gauge transformations 
of the group G. This can be fixed by defining their smearing with an additional holonomy dependance [73, 144], 
but this is however not important for the general discussion since we will only be interested in the action of the 
flux operators on functionals of the holonomies. 
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where the dot denotes the contraction between the indices of the intertwiners and the represen- 
tation matrices. 

The construction of the kinematical Hilbert space can be easily generalized to account for 
the presence of particles [120]. Let us recall that particles in Euclidean space are classified by 
unitary irreducible representations of the isometry group ISU(2). Thus, the physical Hilbert 
space for a massive particle of mass m and spin s is nothing but the module %( m,s ) of the 
associated representation of ISU(2), which is defined by [102] 

n (m,s) = G L 2 (SU(2)) | <j>{xh{e)) = e i6s ^{x), Vx G SU(2)}, 

where h(8) G SU(2) is an element of a given maximal Cartan torus of SU(2). This latter is 
isomorphic to U(l) and can be chosen such that h(0) = e 0T1 , as in (3.26). As %( m ' s ) is a subvector 
space of L 2 (SU(2)), it is naturally endowed with the L 2 measure on SU(2) and therefore inherits 
a Hilbert space structure. Let us remark that %( m '°) ~ L 2 (SU(2)/U(1)), which makes the study 
of particles without spin a bit simpler. Applying the Peter-Weyl decomposition theorem, it is 
straightforward to show that (j) G %( m,s ) can be decomposed in terms of SU(2) representations 
as follows: 

j n=-j 

Therefore, the set {D^j, j G N/2, n G gives an orthogonal basis of T-L^ m ' s \ Elements of 

this set are the point particle analogues of spin networks. 

Now we have all the ingredients to construct the kinematical Hilbert space for gravity coupled 

to point particles. In order to do so, we choose a graph T, sufficiently refined to resolve the 

(p) 

topology of S, with particles sitting at its nodes. We distinguish the nodes Np where particles 
are present and its complement Np . Then the space of cylindrical functions for the coupled 
system is defined by the tensor product 

Cyl(r ; Mi, . . . , M p ) = L 2 (G Lr ) ® ^(g) H (mi ' s A 

between the space L 2 (G Lr ) encoding the gravitational degrees of freedom and the spaces %( mi > Si ) 

(p) 

encoding the degrees of freedom of the particles attached to the nodes Np . At this stage, gravity 
and the particles are not coupled yet. The coupling appears after imposing the constraints. In 
particular, imposing the Gauss constraint leads to the kinematical Hilbert space 

n kin (T; M u ...,M p ) = Cyl(r ; Mi, ... , M p )/G m r. 

The part of the symmetry group acting on the nodes with particles acts simultaneously on 
cylindrical functions and particles states. It is clear that the requirement of the invariance 
under this simultaneous action is responsible for the coupling between the gravitational degrees 
of freedom and the degrees of freedom of the particles. 



Physical scalar product and transition amplitudes 

The kinematical Hilbert space %kin carries an action of the quantum curvature constraint, which 
still needs to be imposed in order to obtain the physical Hilbert space % p h ys - For this purpose, 
we are going to construct a sum over histories (or spin foam) representation of the so-called 
"projector" onto physical states 

= n = / Va ex p (i / tr » ( 3 - 44 ) 
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where a G su(2). For the moment this "projector" is mathematically ill-defined, even if its 
physical content is clear. The fact that it is written as an infinite product of delta distributions 
implies that one has to introduce a regularization to give it a meaning, but also that it cannot 
be a "projector" because V F does not make sense. We will assume in this section that £ is 
a closed Riemann surface of genus g with no punctures. 

The projector can be defined using the Gelfand-Naimark-Segal (GNS) construction (see [19] 
for examples of the GNS construction in LQG). The idea is to find a positive linear functional w 
on the space Cyl(S) of cylindrical functions, which is called GNS state and defines the physical 
scalar product between two elements ^1,^2 £ Cyl(S) as follows: 

(*l,^2>phys = tS7(*T* 2 ). 
This state should enforce the projection onto the physical space, which can be written formally as 

(*i,fH^)ph, s = o. 

Two ingredients are needed in order to make this construction concrete and rigorous: 

• A cellular decomposition of the surface £ which supports the graph T. This decomposition 
is in fact arbitrary, and in the end the physical scalar product will not depend on it. 

• A regularization of the curvature operator F[u]. The classical curvature F[co] is not a cylin- 
drical function and, when promoted to a quantum operator F[u], its action by multiplica- 
tion does not leave the space of cylindrical functions invariant. For this reason, we need 
to introduce a regularization of the quantum curvature constraint such that its action on 
Cyl(S) is well-defined. The projection onto the flat connections enforces the condition 
that the holonomy of the connection around any cell be trivial. It is therefore natural to 
trade the formal definition (3.44) for the discrete operator 

V = HS{g c ), (3.45) 

c 

where g c £ SU(2) is the holonomy around the cell c and the product runs over the finite 
set of cells of the cellular decomposition 16 . This regularized "projector" onto the physical 
states is a well-defined distribution because, in contrast to (3.44), the product is finite. It 
is therefore a good candidate for the GNS state. As a result, the physical scalar product 
between two cylindrical functions and ^2 is defined by 

(*l,*2)phyB = tS7(*T* 2 ) = (*l,7 ? * 2 )km- 

Nonetheless, this is not the end of the story because one has to verify that the product of delta 
distributions does not introduce divergences and that w is a positive state. This indeed turns 
out to be the case. To address the issue of the divergences, it is useful to use the following bound 
on the physical scalar product: 

<*l,*2>ph y s<C]T(2j + l) 2 - 2 S, 
3 

where C is a constant depending on the states. This bound can easily be obtained by decompo- 
sing the delta distributions into SU(2) representations. As a result, the physical scalar product 
is well-defined when the genus g < 2, but it requires an additional regularization when g = 



16 Another regularization was recently introduced in [44]. It is particularly interesting because it mimics the 
situation in four dimensions. 



Spin Foams and Canonical Quantization 



39 




.72 




/ 



' 31 



\ 



Figure 7. Two spin networks states and "J^ which differ only in a small neighborhood. 

and g = 1 . The case of the sphere g = is easy because the divergence is due to a redundancy 
in the product of delta distributions in V, which is a consequence of the discrete analogue of the 
Bianchi identity. One can show that eliminating a single arbitrary cell holonomy from the pro- 
duct in (3.45) makes the "projector" well-defined. The case of the torus (g = 1) is more subtle 
and needs an explicit description of the physical Hilbert space in order to be dealt with [121]. 

As a result, one has a well-defined physical scalar product between cylindrical functions. 
However, one does not have an explicit description of the physical Hilbert space yet. To obtain 
such a description, we need to quotient the kinematical Hilbert space endowed with the physical 
scalar product by the null states. The null states have a vanishing physical norm and form an 
ideal. Their presence is a consequence of the fact that the invariance of £ under diffeomorphisms 
was not taken into account. It is possible to show that the physical scalar product as defined 
above is invariant under diffeomorphisms and therefore compatible with the quotient of %kin by 
the ideal of null states. This quotient defines the physical Hilbert space where a concrete basis 
can be identified. 

It is now straightforward to see the relation between the physical scalar product and the 
Ponzano-Regge model [121] 17 . To illustrate this link, let us consider two spin network states ^1 
and ^2 associated with two graphs Ti and T2 which differ only in a small region, as illustrated 
in Fig. 7. In this small region, T\ has one vertex, whereas T2 has three vertices. Outside of 
this region, the spin networks (i.e. the graphs and the colorings) are exactly the same. A direct 
calculation shows that 



with the notations of the picture. This factor precisely agrees with the one from (3.37) predicted 
by the Ponzano-Regge model! 

The generalization of this construction to the case in which particles are present can be 
done easily since on both sides we already know how the particles are incorporated. Indeed, it 
just amounts to replacing (3.45) by a new projector which enforces that the holonomy of the 
connection around a particle is in the conjugacy class O m associated to its mass m [120]. At 
points where there are no particles, the connection remains flat. Therefore, the computation 
of the physical scalar product leads to the same expression (3.46), where now the dimension dj 
associated to the face crossed by a particle is replaced by the character Xj( m )- 

3.3.3 Positive cosmological constant 

According to the previous discussion, three-dimensional Riemannian gravity with A = is an 
example of a theory in which there is perfect agreement between LQG and spin foam quantization 
techniques. However, the situation is still not quite satisfactory because the Ponzano-Regge 




(3.46) 



17 Note that the relation between the Ponzano-Regge boundary states and spin networks was realized for the 
first time in [71]. 
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model is in general mathematically ill-defined (for closed spacetimes without boundaries) and, 
as a result, the link between the canonical and covariant quantizations stays at a formal level. 
For instance, the evaluation of the path integral with particles associated to closed spacetime 
loops leads necessarily to divergences due to the non-compactness of the gauge group. 

This provides a strong motivation for the study of the relation between LQG and spin foam 
models in three-dimensional Riemannian gravity with a positive cosmological constant, where 
the gauge group is compact and no divergences are expected. However, in this case, the link 
between spin foam amplitudes and kinematical states of LQG is far from being trivial. The 
latter are constructed from representations of the classical group SU(2), whereas the former are 
defined from representations of the quantum group U g (su(2)) (when there are no particles). 
Understanding how the structure of quantum groups emerges from the implementation of the 
curvature constraint can be viewed as a very non-trivial test for LQG. A first step towards a 
resolution of this problem has been done in [122, 127] and we would like to present here the 
main results of this work. 

The idea is to generalize what has been done in the case with A = and to regularize 
the curvature constraint in (3.16) in order to construct a "projector" onto physical states and 
the physical scalar product. Then, we expect to recover the Turaev-Viro amplitudes from 
the physical scalar product between kinematical states. When the cosmological constant is 
non- vanishing, the two constraints (3.16) are equivalent to the requirement that the SU(2) 
connections A± defined by 



be flat. These two connections are in fact the self-dual and anti self-dual components of the 
associated Spin(4) Chern-Simons connection. 

At the kinematical level, the Gauss constraint has already been imposed, and therefore it is 
enough to impose that only one of the connections (3.47) be flat in order to "project" onto the 
physical states, as clearly follows from the following relation 



We will denote generically A one of these connections. If one mimics the construction for the 
flat case (A = 0) presented in the previous subsection, one should impose that the holonomy of 
one of these connections be trivial around each cell. However, contrary to what happens in the 
flat case, the holonomy now depends on the non-commuting variables e and oj, and therefore 
one encounters ordering problems when promoting the classical variables to quantum operators. 
These ordering problems are quite difficult to solve, and to our knowledge no solution giving 
the right physical scalar product has been obtained. 

In [122], the quantization of the (one-parameter family of) classical (kinematical) observab- 
les <? 7 [A], defined as the holonomy of the connection A along the path 7, has been explored. This 
operator is the building block for the construction of the "projector" onto the physical states. 
Its action on the vacuum creates a Wilson line excitation, i.e. it acts simply by multiplication 
by the holonomy of oj along the path, namely 



where g 7 [A] and h~ [oj] denote the holonomies of the Chern-Simons connection A and of the spin 
connection oj, respectively. If however the holonomy g^[A] acts on a spin network, one expects 
a less trivial result. This is because the operator e in g^[A] acts as a derivative operator with 
respect to the components of the connection. As a result, one expects additional contributions 
when the path 7 in g^[A] is self-intersecting or the spin-network contains nodes on (or links 
transversal to) 7. 




(3.47) 



F[A ± ] = (F[w]+AeAe)±\/AT[e,w]. 



57 [A] > I = [oj] , 
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The simplest non-trivial example is the action on a transversal Wilson loop in the fundamental 
representation. To define the quantization of 57 [-A], we would like to quantize each term in its 
series expansion in powers of y/A. But the problem is that quantization of these terms is 
potentially ill-defined due to the non-commutativity of the operators associated to e. A similar 
problem has recently been investigated in [139], where the authors provided a new derivation of 
the expectation values of holonomies in Chern-Simons theory. In the analysis of [139], the same 
sort of ordering ambiguities arises due to the replacement of holonomy functionals under the 
path integral by a complicated functional differential operator. It was shown that the known 
result can be recovered once a mathematically-preferred ordering, dictated by the so-called 
Duflo isomorphism, is adopted. Therefore, following this example, one can also make use of this 
mathematical insight. But in our case the Duflo map does not do all the job. In fact, since the 
ambiguities in the quantization of the Chern-Simons holonomy arise due to the presence of non- 
linear terms in the e field, a second piece of information has to be taken into account, namely 
the quantum action of flux operators in LQG. Combining these two elements leads to a well- 
defined quantization for each term in the perturbative expansion in y/X. Moreover, the series 
can be summed and the result can be expressed in a closed form, leading to algebraic structures 
remarkably similar to those appearing in the skein relation defining the Jones polynomial. 

More precisely, if one concentrates on a single intersection (a crossing) between the path 
defining the holonomy g-y[A] and a transversal spin-network link colored by the fundamental 
representation j = 1/2, one obtains a formula which looks identical to the skein relation (3.23) 
with q = exp(i£p\/£ c ) where £p\ is the Planck length and l c = A _1//2 the cosmological length. Un- 
fortunately, despite the strict resemblance with the Jones polynomial, this result is not sufficient 
yet to recover these polynomials from LQG 18 and should be considered only as a preliminary 
step in this direction. It suggests that it is not impossible to reproduce Witten's results from the 
loop quantization, but the precise way to achieve this still remains unknown. The solution of 
the problem would require a precise regularization of the action of the quantum holonomy <? 7 [A] 
on the holonomy /i 7 [cj] (which is the building block of kinematical states) associated to the 
same path 7, and it is easy to imagine that such a regularization procedure is full of ordering 
ambiguities. 

4 Four-dimensional quantum gravity 

In this section, we focus on the case of four-dimensional gravity. We first start with a brief 
reminder of the basic elements of LQG followed by the introduction of the so-called projected spin 
networks. These states, being an extension of the usual spin network states, allow to put LQG in 
a covariant form and provide the habitat for the boundary states of four-dimensional spin foam 
models. Then we recall some facts about the canonical analysis of the Plebanski formulation 
of gravity and give an overview of various proposals for four-dimensional spin foam models 
available in the literature. These proposals are then compared with the canonical quantization 
and finally the latter is used to deduce some generic restrictions on the form of spin foam 
amplitudes imposed by consistency with the canonical approach. 

18 Note that the quantum parameter q obtained here is slightly different from the one obtained in the Chern- 
Simons quantization, where q = exp(27ri/(fc + 2)). In this last relation, the level k oc lp\/l c is shifted by the 
SU(2) Coxeter number (which takes the value N when the gauge group is SU(iV)). The shift in the level is due 
to a one-loop quantum correction to the path integral, and this correction is exact in the sense that there are no 
higher order contributions. However, one does not see the emergence of this shift in the canonical quantization, 
neither in LQG nor in the combinatorial quantization, where q is in fact determined only at first order in l/k 
when k goes to infinity. Even in [139] where the authors computed Wilson loop expectation values using the LQG 
techniques supplied by the Duflo map, one recovers the Jones polynomials only up to the shift in the definition of 
the quantum parameter. Understanding how to recover the shift in the context of the Hamiltonian quantization 
is certainly a very interesting question, and some arguments have already been given by [104] . 
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4.1 Canonical approach — spin networks and projected spin networks 
4.1.1 Loop quantum gravity in four dimensions 

The canonical loop quantization in 3 + 1 spacetime dimensions follows the same ideas which were 
presented in Section 3.3.2. However, there are several important differences compared to the 
three-dimensional Riemannian case considered there. The first one is that in the Lorentzian sig- 
nature the gauge group in the tangent space is non-compact. This gives rise to certain difficulties 
in rigorously defining various mathematical structures important for the loop quantization. In 
particular, a naive generalization of spin networks to a non-compact group is supplied with 
serious problems [79]. The second difference is that the canonical analysis of the first order 
Hilbert-Palatini formulation of general relativity, or its Plebanski reformulation presented in 
the next subsection, leads to second class constraints which cannot be solved in a covariant 
way. This fact makes the canonical structure of the theory quite involved since the symplec- 
tic structure relevant for quantization is given now by Dirac brackets [10]. Or equivalently, 
one can solve explicitly the second class constraints, which leads to a simple symplectic struc- 
ture [18, 39, 54, 87], but the price to pay is a very complicated form of the first class constraints 
and the loss of explicit covariance. 

Fortunately, there is a simple way to avoid both of these problems. It involves two ingredients. 
First, one changes a bit the starting point. It is now given by the so-called Hoist action [96] 

S M [e,u) = ~ [ euKLe 1 Ae J a(f kl [uj] + ± * F KL [u,]\ , 
1 JM V 7 / 

which differs from the usual Hilbert-Palatini action by the presence of the second term. However, 
this term does not affect the equations of motion and at the classical level the two actions are 
completely equivalent. The coupling constant 7 in front of the second term, which is known 
as the Immirzi parameter, is therefore classically irrelevant. The second ingredient, which is 
responsible for all the simplifications, is to impose a partial gauge-fixing. Namely, one requires 
that e° = 0. This "time gauge" fixes the boost part of the gauge freedom in the tangent space, 
which therefore reduces to the SU(2) gauge group thereby solving the non-compactness issue of 
the Lorentz group. Since the Lorentz covariance is broken anyway, one can also solve explicitly 
the second class constraints. As a result, one finds a canonical formulation very similar to the 
one we had in the three-dimensional Riemannian case: the phase space is parametrized by the 
densitized inverse triad and an SU(2) Ashtekar-Barbero (AB) connection 

E? = \e abc e ljk e{e k c , A\ = r a (E) - 7 u,°\ (4.1) 

where T l a (E) = ^e l jkUJa k is the Christoffel connection compatible with the triad. Note that their 
commutation relation 

{^ a {x),E b J (y)}= 1 5)5 b a 5 2 (x-y) (4.2) 

involves the Immirzi parameter and this fact leads to important physical consequences. On this 
phase space, one should impose a set of first class constraints generating the gauge symmetries: 
three Gauss constraints Gi generate SU(2) transformations, three vector constraints H a are 
responsible for spatial diffeomorphisms, and the Hamiltonian constraint H encodes the time 
repar ametrizations . 

The loop quantization proceeds then as in Section 3.3.2. Assuming that the holonomies of 
the AB connection give rise to well-defined operators in the Hilbert space of quantum gravity, 
one constructs a kinematical Hilbert space %kin with the scalar product induced by the same 
formula (3.41) defined on cylindrical functions. The orthonormal basis in %kin is then provided 
by the usual SU(2) spin networks. 
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On this Hilbert space one still has to impose the vector and Hamiltonian constraints. The 
conceptual meaning of the former is clear - it amounts to considering the equivalence classes 
of embedded spin networks up to three-dimensional spatial diffeomorphisms. As a result, only 
some of their topological properties become important 19 . On the other hand, the imposition of 
the Hamiltonian constraint is much more complicated and has not been properly implemented 
so far. The difficulties start from the form of its classical expression which, in contrast to all 
other constraints, is non-polynomial. Due to this, its quantization is very involved and supplied 
with an infinite set of ambiguities. Even if one chooses a favorite quantization, finding the 
kernel of the corresponding constraint operator seems to be an impossible task. 

Despite these problems, at the kinematical level LQG furnishes a nice picture of quantum 
space where the states are represented by SU(2) spin networks. Moreover, it allows for com- 
putation of spectra of various geometric operators such as area [22, 137], volume [23, 137] and 
length [40, 141]. In particular, these results provide a geometric meaning to the coloring of spin 
networks: the representations attached to the links give the area of transverse surfaces, whereas 
the intertwiners assigned to the nodes encode information about the volume of the surrounding 
region. 



4.1.2 Projected spin networks 

However, the structure provided by LQG is not sufficient for our purpose, which is to relate 
the canonical quantization to the spin foam approach. The problem is that, due to the partial 
gauge-fixing, it operates with an SU(2) subgroup of the initial Lorentz gauge group, whereas 
none of the existing spin foam models make such a reduction. Therefore, we need to understand 
how the results presented above can be reformulated in a Lorentz-covariant way. 

This can be done using the so-called projected spin networks [106] which provide a bridge 
connecting the canonical loop quantization with the spin foam approach. These objects can 
be introduced for any group G and its compact subgroup H, and can be viewed as the usual 
spin networks for the group G with a special choice of intertwiners depending on an element 
of the factor space X = G/H. Namely, let H x be the stationary subgroup of x G X and 7~Lq^ 
be the representation space of an irreducible representation A of G. This representation can be 
decomposed with respect to the subgroup H x as the direct sum 

4 A) = ©^ (4-3) 

j 

and we use p and rrij to label basis elements in and 71$ , respectively. With these 

definitions, the class of intertwiners giving rise to projected spin networks is given by 

k=l rnj 1 ...mj h fc=l 

where g x is a representative of x G X in G, A and j are a set of representations of G and H , and 
/ L \ 



G Inv 07$ 



Ok) 

k=l 



is an invariant i?-intertwiner. Thus, besides the L representations Afc coupled by the intertwiner, 
it is parametrized by L representations jk of the subgroup H appearing in the decomposition 



1 Nevertheless, even at this level there seem to remain some unsolved issues which, in particular, lead to the 
fact that the resulting Hilbert space is still non-separable [92]. 
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Figure 8. Structure of a projected spin network. The link i carries a representation \g of G, which 

(I s ) ( i t ) 

is projected at its source node onto H H and at its target node onto H H . The node n is labeled 

by an H- invariant intertwiner i n . The projected spin network state depends on the group elements gg 

associated to the links and on the normals x n associated to the nodes. 



of Xk and by an invariant .ff-intertwiner between these jk- Note that by itself the intertwiner I 
is not invariant. Although it does stay invariant under transformations from the subgroup H x , 
under general G-transformations it transforms in a covariant way 




In other words, the tensor field X{x) is invariant provided the transformation acts not only on 
the tensor indices, but also on the argument x. Such an invariance property was called in [9, 11] 
relaxed closure condition. 

Combining the intertwiners (4.4) with the group elements in representations \g as in (3.43), 
one obtains a projected spin network function 

S {r,xss,m\9^n] =®T>M(gt) • ®I {lnM {x n ). (4.6) 

l n 

It should be clear that the projected spin networks carry the following coloring: each link is 
colored by a representation A^ of G and by two representations of H, and jh associated with 
the source and target nodes of the link, and each node comes with an /^-invariant intertwiner i n . 

The idea is that the matrix element D^^(^) is "projected" from the left to Hj^ ^ and from the 

tit) 

right to % H l (see Fig. 8). Since the subgroup H is supposed to be compact, this projection 
selects a finite matrix block in the original (infinite-dimensional if G is non-compact) matrix. 
Such a construction allows to avoid the problems arising for the usual spin networks defined 
for a non-compact gauge group [79] since all traces become effectively finite dimensional. In 
particular, the projected spin networks are well-defined on the trivial group elements ge = I. 

It is important to remember that the projected spin networks depend not only on the group 
elements gi G G, but also on the elements of the factor space x n £ X. Due to this additional 
dependence, and despite the projection to a subgroup, the functions (4.6) are G-invariant 

where £ s and 1% denote the source and target nodes of the link I. This equality can be easily 
verified using (4.5). In fact, the projected spin networks provide an orthonormal basis in the 
space of all gauge invariant functions of a G-connection and an X-valued scalar field [13]: the 
connection appears through its holonomies along the links of the graph T and the scalar field 
through its values at the nodes. The corresponding Hilbert space structure is induced by the 
scalar product defined as usual by [106] 

r L r 

(*r,/il*r,/ a ) = / Y[ d ^f^9i,---,gL r ,x 1 ,...,x Nr )f 2 (gi,...,gL r ,xi,...,x Nr ), (4.7) 
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where ^r,/ is a (generalized) cylindrical function defined by a graph T and a function / : 
g l t x x n f — y C) and the gcalar p roc i uc t (4.7) can be checked to be independent on the choice 
of x n £ X due to the left and right invariance of the Haar measure. 

In the case of four-dimensional gravity we are interested in, G should be taken to be either 
the Lorentz group or Spin(4), depending on the signature, and H = SU(2). The projected spin 
networks appear then as a habitat for kinematical states in the Lorentz-covariant approach to 
the loop quantization (see, for example, [16]). Moreover, we will see that they also represent 
the boundary states of all spin foam models of gravity which have been introduced up to now. 
Thus, they provide the most direct link between the canonical and spin foam quantizations. 



4.1.3 Lorentz-covariant form of loop quantum gravity 

The advantage of projected spin networks can be seen immediately from the fact that they allow 
to formulate LQG in a Lorentz-covariant way. To see how this comes about, we should better 
understand the origin of their arguments, the connection and the X-valued scalar field. The 
latter has a very clear geometric meaning: it defines the normal to the foliation used to perform 
the 3 + 1 decomposition of spacetime. In particular, choosing this normal to be time-directed, 
one reproduces the time gauge used in the canonical formulation of the Hoist action leading to 
LQG. 

The choice of the connection defining the holonomy in the first argument of (4.6) is a more 
subtle issue. The most evident possibility would be to take the usual spin connection a/" 7 . 
However, due to the second class constraints mentioned in the beginning of this section, it has 
very non-trivial commutation relations. As a result, the holonomies constructed from lo ij do not 
diagonalize the action of the area operator [17], which is in drastic contrast with the situation 
in LQG. Due to this, the spin connection does not seem to be suitable for the loop quantization. 

Fortunately, one can construct other candidates for the connection to be used in holonomy 
operators and in projected spin networks. In particular, in [7] using the approach based on the 
Dirac brackets and in [87] by explicitly solving the second class constraint (see also [54]), it was 
found that there exists a Lorentz connection A IJ which provides a generalization of the SU(2) 
AB connection (4.1) in the following sense: 

• In the time gauge x = 0, h reduces to the AB connection 
a jj = fo, 1 = 0, J = i, 

a time gauge 1 k A%, I = i, J = j. 



• Its commutation relations generalize the symplectic structure (4.2) in a natural way: 

{ A 1 /, A« L } D = 0. { A 1 /, e bcd e^e L d } D = 7 6 b a I IJ ' KL (x) 5 3 (x - y), (4.8) 

where lj/i{x) = S^ K dj^ — 2ax S^x^-i is the projector on the SU(2)a; part of the Lorentz 

group, so that in the time gauge x 1 = Sq one recovers again (4.2). In particular, the area 
operator evaluated on holonomies of this connection leads exactly to the same spectrum, 
given by the Casimir of SU(2) [7, 86], as in LQG. 

The explicit expression of A through the spin connection and the normal to spacelike slices is 
given by 

= lKL(l ~ l^a L + 2(1 + ^dax'l (4.9) 
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Moreover, in [15] it was shown that the loop quantization of the symplectic structure (4.8) 
using projected spin networks gives a quantum theory equivalent to LQG but formulated in 
a Lorentz-covariant form. The clue to this result is a constraint satisfied by the connection A IJ : 

(x [J 6^x l] )aK l = ax^dax 1 ^, (4.10) 

as can be easily verified using (4.9). The operator acting on the connection on the left-hand 
side is the projector on the boost part of the Lorentz algebra orthogonal to su(2) x . Therefore, 
for a constant x , this relation implies that only the su(2) x part of the connection is non-trivial 
and its holonomy ge(A) belongs to the subgroup SU(2)a;. Since by definition such operators do 
not mix various terms in the decomposition (4.3), one obtains the following result 20 



Thus, the projected spin networks evaluated on the connection A are non-vanishing only for 
those SU(2) spin labels, which are the same on both ends of every link, jl = 3i = Hi do 
not depend on the representation labels A^, and coincide with the usual SU(2) spin networks 
evaluated on the AB connection. Furthermore, the constraints (4.10) also affect the scalar 
product, which is not given by the simple formula (4.7), but should involve the <5-functions of 
constraints in the measure. The correct kinematical scalar product is then given by 

r L r _ 

(*r,/il*r,/ 2 )kin = / TT [d^feffe/ g£g x „))h(gt,x n )f 2 (gi,x n ), (4.12) 

where Snig) is the 5- function on the group with support on the subgroup H such that 
J G dgdH(g) = fjjd-9- It is clear that for the states (4.11), the formula (4.12) reduces to the 
usual scalar product on the SU(2) spin networks. 

This construction identifies the kinematical Hilbert spaces of LQG and Lorentz-covariant 
quantization. Besides, it provides an important lesson: the projected spin networks form an 
extended state space which should be reduced to the kinematical Hilbert space by imposing 
certain constraints on their arguments, i.e. on the connection entering the definition of the 
quantum holonomy operators. In fact, these constraints (in particular (4.10)) are nothing else 
but the secondary second class constraints conjugated to the simplicity constraints, playing the 
prominent role in the spin foam approach. We derive them in the next subsection by studying 
the Plebanski action. 



4.2 Plebanski formulation and its canonical analysis 

Most approaches to four-dimensional spin foam models take as a starting point the formulation 
of general relativity as a constrained topological field theory. As we already mentioned in 
Subsection 2.3, this can be done by considering the Plebanski modification of the BF action (2.4). 
The Plebanski action with non-vanishing cosmological constant is given by 

Spi[u,B,X]=J (rr(B A F) + | Tr(*B A B) + ^ Tr(A • B A B)\ . (4.13) 

In this expression, the density tensor \ IJKL = \vA\ KI A is symmetric under the exchange of the 
pairs [IJ] and [KL], and satisfies the tracelessness condition £ijklX ijkl = 0. It generates 
the 20 simplicity constraints ensuring that the 36 components of B\^ v reduce to 16 components 

20 We will use h to denote holonomies of the subgroup H C G. 
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associated to a tetrad field e^. To understand how this comes about, let us look at the equations 
of motion derived from (4.13). They are given by 



DB = dB + [u,B] = 0, 

F u + l -(Ke IJKL + \ IJ 
B lJ AB KL = aVe UKL 



L 



where 



1 



}IJ A T}KL 



V = ^e IJKL B 1J AB 1 

and a is as usual the sign distinguishing the Riemannian and Lorentzian signatures. The last 
equation of motion is known as the simplicity constraint. In the non-degenerate case, i.e. when 
V 7^ 0, one can show that this condition implies that the B field can be written as 

• topological sector (I 1 * 1 ): B IJ = ±e 7 A e J , 

• gravitational sector (IF* 3 ): B IJ = ±\e IJ KL e K A e L . 

Plugging the solution for the gravitational sector back into the Plebanski action, one recovers the 
usual Hilbert-Palatini formulation. This demonstrates that, when Plebanski theory is restricted 
to this sector, the two formulations are equivalent at the classical level. 

It is useful to note that, again for V / 0, the simplicity constraints can be equivalently 
rewritten as 

B u A b kl = aV£ UKL ^ EijklB^B^ = aVe^pa. (4.14) 

Of these two equivalent forms, it is the second one that will be important for our purposes. It 
can be obtained from an action similar to (4.13) which takes the form 



Sp\[u,B,X] 



J d A x [rr{B^F pa ) + ~ TrfrB^B^j + -X^Tr^B^B^) 



(4.15) 



The difference is the nature of the Lagrange multiplier field A. Whereas before it lived in the 
tangent space, here it is a spacetime tensor, satisfying the same symmetry properties as X IJKL . 
The canonical analysis of this theory was carried out for the first time in [49], and further 
developed in [12], whereas the canonical structure of the action (4.13) has been elucidated 
in [14]. The canonical analysis of (4.13) turns out to be simpler than that of (4.15), and not 
surprisingly can be obtained from the latter by solving some of the auxiliary constraints. Here 
we present the results of the original analysis of [49]. 

To perform the canonical analysis of the action (4.15), it is convenient to add new non- 
dynamical variables [i a and 7r a in order to enforce the vanishing of the momenta conjugated 
to Bq g . This can be done by adding to the action the following term 

J d 4 x(Tr(7r a d B 0a ) + Tr(/x a ^)) . 

The variables u)q, X and B$ a are non-dynamical (they appear without their time derivatives), 
but the fact that B$ a appears quadratically in the action prevents us from treating it as a true 
Lagrange multiplier, and this is the reason for which we have to use this additional term. The 
action can then be written in its canonical form as 

S P i[o;,B > A,7r,/i]= I dt f d 3 z(T¥(P a d o; a ) +Tr(Tr a d B 0a ) - H T ), (4.16) 
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where we have introduced P a = e B bc and Ht is the total Hamiltonian given by a linear 
combination of primary constraints. The results of the Hamiltonian analysis can be summarized 
as follows. The total list of first class constraints is given by 

Q = D a P a -[P a ,B 0a }^0 : 

K a = £ a6c Tr(7T 6 P C ) « 0, 

/Co = Tr(7r a ,B 0a ) « 0, 

V a = TiiP b d a uj b ) - d b Tr(u a P b ) + Tr(« b d a B 0b ) - <9 6 Tr(^P 0a ) « 0, 
P = s abc Tr(B 0a F bc ) + ATr(*P a P a) « 0. 

Among these first class constraints, the first three are primary and the last two are secondary 
arising from the requirement that lC a and /Co be preserved in time. The constraint Q is the 
generators of the Lorentz gauge transformations, while the remaining constraints are related to 
the generators of the diffeomorphisms. 

Apart from the first class constraints, the phase space variables of the Plebanski theory are 
also subject to second class constraints. These are given by 



k)\ = Tr(B ob Tr a ) - i 6 a b lC = 0, 


(4.17a) 


$(P,vr) afe = Tr(p( a 7r 6 ))=0, 


(4.17b) 


<S>(B,B) ab = Tr(*B 0a B ob ) = 0, 


(4.17c) 


PA = Tr(*5 06 P a ) - 8%V = 0, 


(4.17d) 


${P,P) ab = Tr(*P a P 6 ) = 0, 


(4.17e) 


^ab = e «i(«Tr(*D c P 6 )s od ) = 0. 


(4.17f) 



One can recognize that the second class constraints $>(B,B), $(P,P) and $(P, P) are the 
various components of the simplicity constraints (4.14). Moreover, one can solve explicitly the 
first four constraints, which removes 14 out of 18 components of Bq^ and nfj. The remaining four 
components represent the usual lapse J\f and shift M a and their conjugated momenta, which are 
treated in this formulation as dynamical variables. In particular, the first class constraints fC a 
and /Co are responsible for their transformations under diffeomorphisms. Reinterpreting J\f 
and M a as Lagrange multipliers, one establishes a direct contact with the canonical formulation 
of the Hilbert-Palatini action [12l 21 . 

In contrast, the remaining second class constraints, $(P, P) a and ^ , cannot be solved in 
a Lorentz-covariant way. The meaning of the former is to say that Pf, is expressed in terms 
of the triad Ef and the normal x introduced in the previous subsection. On the other hand, 
the latter is a restriction on the connection co^ J which allows to express six of its components in 
terms of the triad and the normal. In fact, it is a secondary constraint obtained by stabilization 
procedure from $(P, P) ab . They are conjugated to each other and together form a second class 
system. Thus, when it is said that the simplicity constraints are of second class, it is important to 
remember that this is due to the presence of these secondary constraints which do not commute 
with the primary ones 22 . 

The presence of the second class constraints modifies the symplectic structure of the theory: 
the usual Poisson bracket on the phase space has to be replaced by the appropriate Dirac 
bracket which is simply the symplectic structure induced on the constraint surface. The Dirac 

21 Note that these results apply only in the gravitational sector. In the topological sector the canonical structure 
changes drastically. But the two sectors stay disconnected from each other. 

22 For example, in the topological sector there are also secondary simplicity constraints, which however commute 
with the original simplicity. As a result, all constraints appear as first class [105]. 
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brackets of the variables parametrizing the phase space have been computed in [12]. These 
results in particular imply that P a and uj a are not canonically conjugated anymore and the 
spin connection becomes non-commutative. As we will see in the next subsection, instead of 
quantizing the resulting complicated symplectic structure, in the spin foam approach one follows 
a different strategy which ignores the presence of the secondary constraints. We will discuss 
implications of this approach in Section 4.4, whereas in Section 4.5 we will try to formulate 
a spin foam quantization taking into account the canonical structure we have just described. 

4.3 Four-dimensional spin foam models 

In this subsection we present the spin foam models of four-dimensional gravity which have 
received the most attention. These are the Barrett-Crane (BC) model [32, 33], the Engle- 
Pereira-Rovelli-Livine (EPRL) model [66] , and the Freidel-Krasnov (FK) model [75] . We outline 
their derivations, give the resulting constructions, and discuss their relations to the canonical 
quantization. In particular, we express their boundary states in terms of the projected spin 
networks introduced in Section 4.1.2. 

All these models are based on the strategy outlined in Section 2.3. Namely, they are ob- 
tained by implementing a quantum version of discretized simplicity constraints of the Plebanski 
formulation in the partition function of four-dimensional BF theory. The latter can be found in 
the same way as its three-dimensional analogue in Section 3.3.1. The only difference is that in- 
stead of discretizing the triad one-form as in three dimensions, one has to discretize the bivector 
fields Bj^. This can be done by integrating the bivectors along the two-cells (triangles) of A 
dual to the faces / G A*, in order to obtain Lie algebra elements Bj J . Then one can proceed as 
in the three-dimensional case and integrate over the discretized B field in the partition function 
to obtain the analogue of (3.34): 




(4.18) 



where the group G is now either SL(2, C) or Spin(4). Using the Peter- Weyl decomposition and 
integrating out the group elements, one arrives at a representation of the partition function 
similar to (3.36), where the vertex amplitude is given by the evaluation on a flat connection of 
the spin network dual to the boundary of a four-simplex and represented by the pentagon graph 




(4.19) 



This result also tells us that the space of the boundary states in quantum BF theory is spanned 
by the usual spin networks of the gauge group G. 

The next and the main step is to properly discretize, quantize and incorporate the simplicity 
constraints. This is where different spin foam models put forward different proposals and start 
to diverge from each other. The problem is that in general there is no unique, well-defined 
procedure to make such a "quantum reduction" because from the very beginning it goes against 
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the principles of Dirac's quantization, which state that the second class constraints should be 
implemented prior to quantizing the initial symplectic structure. Nevertheless, keeping this 
difficulty in mind, but proceeding with the suggested strategy, one can build interesting models 
and some of them turn out to be closely related to the results of the canonical approach. 

4.3.1 BC model 

The BC model was the first spin foam model of four-dimensional gravity. The first step towards 
its definition is the discretization of the quadratic simplicity constraints (4.14). In the discrete 
setting, they can be split into three family of constraints, depending upon the position of the 
triangles defining the discretized B field. The various constraints are given by 

• diagonal simplicity: £jjKLBj J B^ L = 0, if two triangles are the same, 

• cross simplicity: sukl^Y Bf L = 0> if two triangles share an edge, 

• volume simplicity: eijklBY B^, L = ±V V , if two triangles meet only at a vertex, 

where V v is the volume of the dual four-simplex. It is convenient to distinguish between these 
three types of constraints because they play rather different roles in the spin foam quantization. 
Besides, note that the bivectors Bf possess a clear geometric interpretation: their norm is equal 
to the area of the triangle dual to the face /, and their tensor structure encodes its orientation 
in spacetime. Given this interpretation, the bivectors should satisfy an additional constraint, 
which originates from the fact that the triangles forming the boundary of a tetrahedron r cannot 
be arbitrary. They are restricted to fulfill the following condition: 



which is known as the so-called closure constraint. One can show [108] that this constraint 
together with the diagonal and cross simplicity implies the volume simplicity. Due to this, the 
volume simplicity is often ignored and we remain with the problem of implementing the other 
constraints in the boundary states and in the partition function of the topological BF theory. 

To translate the constraints to the quantum level, one identifies the classical bivectors Bj J 
with the generators T IJ of the Lie algebra of the gauge group G. This identification is justified by 
the observation that in BF theory the B field is canonically conjugated to the spin connection. 
As a result, in the boundary Bilbert space it is represented by a derivative operator whose 
action on the spin networks is equivalent to an insertion of the generator of the Lie algebra in 
the representation associated with the given face. This can be clearly seen from the action of 
the smeared B field on the holonomy in representation A. This action can be schematically 
represented as 



where S is a two-dimensional surface transverse to the link I = l\ U £2- Thus, the quantization 
map 




(4.20) 



fCT 





Bf ^ T f = T (A ^ 



(4.21) 



used to quantize the simplicity constraints, is a simple consequence of the symplectic structure 
of BF theory, in accordance with our strategy which implies the ignorance of the second class 
constraints in the process of quantization. 
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Once this identification has been made, the simplicity constraints are imposed as conditions 

d-5 (rXf) = (4.22) 

on the boundary states of the BF theory, which has the effect of fixing the representations and 
the intertwiners coloring the graphs. The easiest constraint to impose is the diagonal simplicity, 
which involves only the bivectors referring to one face. As a result, it restricts the allowed 
representations A/. More precisely, using the map (4.21), one finds that it leads to the condition 
that the second quadratic Casimir operator of G should vanish: 

cg\\ f )=T f *T f = 0. (4.23) 

The representations satisfying this condition are commonly called simple and will be denoted 
by X s . To display more explicitly what the simple representations are, we distinguish the cases 
in which G is either Spin(4) or SL(2,C). For the former, the unitary irreducible representations 
are labeled by two half-integers, A = and for the Lorentz group, representations of the 

principal series are given by A = (n,p) with n 6 Z/2, p £ 1, We use the principal series only 
because they are the only ones which appear in the Plancherel formula for the Lorentz group. 
The two quadratic Casimir operators are given, respectively, by 23 

C Sn(4)(i + >.T ) = 2j+(i+ + 1) + 2j-(j- + 1), Cg? (2>c) (n,p) = n 2 - p 2 - 1, 

C Sp!n(4)(i + '^) = 2 ^ + (i + + !) - 2 -Ttr + 1), C^ (2C) (n,p) = 2np. 

Thus, the solution to the condition (4.23) is provided by the representations with j + = j~ in 
the Riemannian case and n = (or p = 0) in the Lorentzian case. 

The next step is to impose the cross simplicity constraint. It involves different faces dual to 
triangles sharing an edge, and these triangles are in the boundary of a tetrahedron r. Since a 
tetrahedron is represented in the quantum theory by an intertwining operator, the imposition of 
the cross simplicity requires that we look for intertwiners that are compatible with the coupling 
of simple representations. More precisely, since it is always possible to decompose an intertwiner 
with arbitrary valency into the sum of three- valent intertwiners, the cross simplicity is equivalent 
to the condition that the intermediate representations of the decomposition be simple as well. 
A solution to this condition, called BC intertwiner, was proposed in [33] and was shown to be 
unique in [133]. 

Finally, it remains to implement the closure constraint, which simply amounts to imposing 
the invariance of the intertwiners with respect to the group G. Since however the original 
boundary states of BF theory represented by G spin networks already possess this invariance, 
there is actually nothing to impose. As a result, we arrive at the boundary state space of the 
BC model, consisting of spin networks for the group G with representations restricted to be 
simple and intertwiners fixed to coincide with the BC intertwiner. We can denote such states 

(r,A») 

The restrictions on representations and intertwiners obtained by imposing the simplicity 
constraints should now be inserted into the spin foam representation of the partition function for 
BF theory. The result is supposed to give a model for four-dimensional quantum gravity. Since 
the constraints implemented in this way do not change the summand, but only the set of group 
theoretic data which one sums over, the vertex amplitude is given by the same formula (4.19) 
as in BF theory, but with arguments satisfying the above requirements. Denoting by T a this 
pentagon graph (let us remind that it is dual to the boundary of a four-simplex), we therefore 
can write 

J^ C (X')=S^ ) ]S\. (4.24) 



23 Note that the Riemannian and Lorentzian cases are related by the formal identification 2j + 1 = n ± ip. 
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This defines the BC spin foam model. Strictly speaking, to complete its definition one should 
also provide the face and edge amplitudes. From the procedure outlined above, one finds that 
A e = 1 and Af = d\ f . However, one can expect on a general ground that this choice is 
modified by the implementation of the simplicity constraints. See in particular the discussion 
in Section 4.5.1 where this choice is related to a non-trivial measure in the canonical path in- 
tegral [41]. Although several proposals for these lower dimensional amplitudes have been put 
forward, none of them was sufficiently convincing. Thus, it is still considered as an impor- 
tant open issue, especially because the convergence properties of spin foam transition functions 
strongly depend on the choice of Af and A e [57, 129]. 

The BC model has been extensively studied in the literature. In particular, for some special 
choice of the face and edge amplitudes it has been shown to be finite [57, 125] and to possess 
an interpretation in terms of a group field theory [130]. It also was generalized in order to allow 
for the presence of timelike bivectors [128] where however the finiteness properties get lost. 

At the same time, no relation with the canonical quantization has been established, except 
one remarkable fact. It has been noticed in [15, 106] that the boundary states of the BC model 
can be expressed in terms of the projected spin networks. Indeed, using the representation for 
the BC intertwiner found in [76], it is easy to show that 



i.e. the boundary states are obtained from general projected spin networks by fixing the SU(2) 
representations one projects onto to be the trivial ones jf = j\ = 0, which also forces the SU(2) 
intertwiners to be trivial i n = I, and by integrating over all factor space elements living at the 
nodes 24 . The integral is necessary in order to ensure the invariance of the intertwiners, which 
should be contrasted with the covariance property (4.5) valid without integration. Due to this, 
one can rewrite the BC vertex amplitude (4.24) in the following way: 



This integral representation was very useful in the asymptotic analysis of the vertex amplitude 
and plays an important role in attempts to overcome the shortcomings of the BC model. Al- 
though the relation (4.25) suggests a possible connection with the canonical quantization, where 
the projected spin networks have been introduced for the first time, it does not explain how to 
reproduce the restrictions on the labels on the canonical side. As a result, the lack of a clear 
connection with the loop approach called for modifications of the BC model [11, 68]. 

In fact, this was not the only reason. A lot of criticism has been raised concerning the 
ability of the BC model to reproduce correctly quantum simplicial geometry. It seems that the 
asymptotic (i.e. the large spin limit) behavior of the vertex amplitude [25, 38, 81] does not lead 
to the exponential of the Regge action, but rather to crippled geometric configurations. Another 
serious issue was the inability of the BC model to reproduce the graviton propagator [6] . Most 
of these problems seem to be related to the so-called ultralocality property of the BC model, 
which originates in the uniqueness of the BC intertwiner and results in the fact that the simplices 
do not talk to each other. This feature can be traced back to the way in which the simplicity 
constraints have been implemented. Thus, if one wishes to define a new model which is not 

24 In fact, in [76] the notion of simple spin networks representing the boundary states of the BC model has been 
generalized to any dimension d > 4. In particular, for a general gauge group G, a representation A is called simple 
if its decomposition (4.3) into irreducible representations of the maximal compact subgroup H contains the trivial 
representation. Remarkably, such simple spin networks are given by the same formula (4.25) in all dimensions, 
provided G is the corresponding group Spin(d). 




(4.25) 
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plagued by the same problems, the constraints have to be imposed in an alternative way . Such 
alternatives have been suggested by the EPRL and FK models, which are presented in the rest 
of this subsection. 

4.3.2 EPRL model 

The EPRL model is a straightforward generalization of the BC model based on three new inputs: 

• linearization of the simplicity constraints, 

• inclusion of the Immirzi parameter, 

• imposition of the simplicity constraints in a weak sense. 

Let us present these new ideas and their effects on the spin foam quantization. 

As we know from the classical analysis, the Plebanski formulation is equivalent to four- 
dimensional gravity only in one of the sectors of solutions to the simplicity constraints. Besides 
the gravitational sector, there is a topological sector which may "interfere" with the physical one 
at the quantum level. The BC model does not distinguish between the two sectors and therefore 
one cannot be sure that one describes pure gravity. On the other hand, there is a simple way 
to resolve this problem. It relies on the geometric meaning of the diagonal and cross simplicity 
constraints. What they simply mean is that the four triangles described by bivectors Bf, which 
belong to the same tetrahedron r D /, lie in the same hyperplane. This is equivalent to the 
condition that there exists a vector x T , normal to the tetrahedron, such that [11, 68] 

{*B f ) IJ {x T )j = 0. (4.26) 

This linear version of the simplicity constraints turns out to be stronger than the original 
one (4.14). It is satisfied only by the gravitational sector of solutions for the B field. 
The topological sector 1^ can also be described in a similar way, but requires a different linear 
constraint which can be obtained from (4.26) by removing the star operator. The new models, 
EPRL and FK, are both based on the linearization (4.26). One of the consequences of this is 
that from the very beginning, the formulation features the normals x T . They turn out to coin- 
cide with the elements x n G X appearing as arguments of the projected spin networks (recall 
that the tetrahedra r G A are dual to the spin foam edges e G A* which in turn correspond to 
the nodes n G T of the boundary spin networks). This fact facilitates the comparison with the 
canonical framework. 

The second point is to incorporate the Immirzi parameter into the quantization scheme. 
This is, of course, motivated by the hope to find a spin foam model consistent with LQG, whose 
results depend on this (classically irrelevant) parameter in a non-trivial way. This goal can be 
easily achieved by replacing the B field in the BF part of the action by its Hoist generalization 
B+ , y~ 1 *B. Although this replacement does not affect the partition function of the unconstrained 
theory, it does change the quantization of the bivectors Bf. As a result, the previous quantization 
map (4.21) takes now a 7-dependent form, 

B f + l*B f M> T f & B f ^ (T f - 1 *2>) , (4.27) 

provided that 7 2 7^ a. Using this map in the simplicity constraints, one obtains a quantum theory 
which also depends non-trivially on the Immirzi parameter. What remains to be understood is 
how the constraints are actually imposed. 

25 One should however mention some recent developments in the non-commutative flux formulation of group 
field theory pointing towards a possible reconsideration of the BC model [31]. We refer the reader to this last 
paper for a discussion of various arguments that have been put forward in this context. 



54 



S. Alexandrov, M. Geiller and K. Noui 



We saw that the problems of the BC model were related to an over-imposition of the simplicity 
constraints. Mathematically this is explained by the fact that, after quantization (4.21), they 
become non-commuting operators. In other words, they represent a second class system, and 
imposing them in a strong way as in (4.22), one actually requires that all their commutators be 
vanishing as well. To avoid this problem, it has been suggested [68] to replace (4.22) by a weaker 
condition where one requires only the vanishing of their matrix elements: 

(V|l>|*) = 0. 

Moreover, the initial quadratic simplicity constraints are supposed to be replaced by their li- 
nearized version (4.26) and quantized using the 7-dependent map (4.27). In fact, two other 
methods have also been suggested: using the so-called master constraint [66] and using the 
coherent states [56] playing an important role in the FK model. All three methods lead to the 
same construction and we present only the final result encoded in the following two equations: 



Several comments concerning this result are in order: 

• In all three methods, not all the simplicity constraints are treated on an equal footing. In 
fact, the diagonal simplicity is still imposed in a strong way as in (4.22). The distinguishing 
role of this constraint is justified by the fact that it lies in the center of the algebra formed 
by the quantum constraint operators, i.e. it commutes with all the other constraints, and 
can therefore be considered as first class. It is this constraint that, after using (4.27), leads 
to the first equation (4.28a). 

• The second equation (4.28b) is a consequence of the cross simplicity. The linearized version 
of this constraint (4.26) depends explicitly on the normal to the tetrahedron r which can be 
considered as an element x T £ X. In its quantum version, this dependence appears through 
the Casimir operator Ch, since the subgroup H = SU(2) appearing in this equation is the 
stabilizer of x T . (This is why we put two indices on j T f.) 

• The presence of the SU(2) subgroup and the normal x T immediately implies that the solu- 
tion to the quantum simplicity constraints will lie in the class of projected spin networks. 
The SU(2) representations j T / are then identified with the representations jf and j\ at- 
tached to the ends of links of the boundary graph and appearing in the construction of 
Section 4.1.2. Equation (4.28b) fixes these representations in terms of A/ and thereby 
constrains the intertwiner associated with the tetrahedron r. 

• Strictly speaking, the two equations (4.28) do not have common solutions neither in the 
Riemannian nor in the Lorentzian case. It is assumed however that they may acquire 
quantum corrections due to ordering ambiguities in the Casimir operators. Moreover, 
these corrections are assumed to be such that the equations do have non-trivial solutions. 
What is presented below corresponds to one particular choice of such corrections, but it 
should be mentioned that there are also other interesting possibilities [8, 62]. 

One can show that it is possible to adjust the quantum corrections to the Casimir operators 
such that the equations (4.28) have the following solution in the Riemannian and Lorentzian 
case, respectively, 




(4.28a) 



cg ) (\ f )-2 1 C H (j Tf ) = 0. 



(4.28b) 





(4.29) 
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We see that the representations Xf of the group G are completely determined by the SU(2) 
representations j T t. An important consequence of this fact is that even in the Lorentzian case 
the spectrum of allowed representations becomes discrete, despite the non-compactness of G = 
SL(2, C). As a result, the same pattern holds for various geometric operators [62]. On the other 
hand, in the Riemannian case there is a clash between the discrete nature of representation labels 
and the continuous Immirzi parameter. Indeed, for X(j) defined in (4.29), it belongs to [N/2] 2 
only if the Immirzi parameter is a rational number. This is a very unexpected result which has 
no clear justification coming from the canonical quantization 26 . 

The restrictions on representations (4.29) imply that the boundary states of the EPRL model 
coincide with the following particular class of projected spin networks: 

This formula represents the kinematical setup of the EPRL model. In principle, one should also 
impose the closure constraint (4.20), which is not satisfied automatically anymore because of 
the explicit dependence on x T . This can be easily achieved by integrating over the normals as 
in (4.25). This is what has been done in the original version of this model [66]. However, as 
was argued in [9, 11] and will be justified below in Section 4.5.2, the closure constraint must 
be relaxed and the integration over x T should be omitted. And it seems that now there is 
a convergence towards this viewpoint, see for instance [30, 42, 138]. 

The dynamics of the model is encoded in the EPRL vertex amplitude which is as usual 
obtained by restricting the labels in the BF vertex (4.19) according to (4.29). Thus, it is 
completely determined by our kinematical result and appears as the simplex boundary state 
(with the closure constraint implemented) evaluated on a flat connection. Equivalently, it can 
be written as the following integral: 

A7 RL (J, •) = / II d 9n Sptosat) *o] • (4-31) 

J n 

We postpone a more detailed discussion of the relations between this model and the canonical 
quantization to Section 4.4. 



4.3.3 FK model 

Finally, we present the FK spin foam model. It is based on the same set of ideas which underlie 
the EPRL model, namely, the linearization of the simplicity constraints, the inclusion of 7 and 
the imposition of the constraints in a weak sense, but this last step is realized in a different way. 
Instead of imposing the simplicity on the boundary states, the FK model uses a coherent state 
decomposition in the path integral of BF theory and imposes the constraints on the semi-classical 
bivectors representing the coherent states. Let us explain the main steps of this construction by 
focusing on the Euclidean theory. 

We start by recalling some facts about coherent states on a group. For the SU(2) group, 
these states are defined by acting with an SU(2) element on the highest weight state \j,j) in 
a representation of spin j 

\h,j)=B^(h)\j,j). 

26 Note, however, that a similar phenomenon occurs in three-dimensional Riemannian gravity, where the level k 
(related to the cosmological constant) of the Chern-Simons theory is discrete at the quantum level. It is clear 
from the covariant point of view that k must be discrete but, to our knowledge, no explanation coming from the 
canonical point of view exists. In the Lorentzian regime, there are no restrictions on the values of k. 
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They form an overcomplete basis and can be used to write a decomposition of the identity as 
Ij=dj [ dn\j,n)(j,n\, (4.32) 



where we took into account the possibility to reduce the integral to the quotient space § 2 = 
SU(2)/U(1). The advantage of these states is that they possess a geometric interpretation: the 
state \j, n) describes a vector in M 3 of length j and of direction n, given by the action of h n on 
a unit reference vector, 

(j,n\J\j,n) = jn. (4.33) 

Furthermore, they minimize the uncertainty of the quadratic operator J 2 and therefore can be 
considered as semi-classical states. To get the coherent states for the Spin(4) group, which we 
are really interested in, it is sufficient to pair two SU(2) coherent states 

|A,n) = \j + ,n + ) ® |j + ,n + ), 

where A = (j + ,j~) and n = (n + ,n~). It is clear that they satisfy the decomposition of the 
identity analogous to (4.32). 

Let us now consider the partition function of BF theory obtained by expanding (4.18) in the 
sum over representations, but before performing the integration over the holonomies g e . Since 
an edge e G A* connects two four-simplices sharing the tetrahedron r G A dual to the edge, 
the holonomies can be represented as a product of two "half- holonomies" , g e = g aT g~} , where 
g GT goes from the center of the four-simplex a to the center of the tetrahedron r. Inserting the 
decomposition of the identity in terms of the Spin(4) coherent states between the two group 
elements of this representation for all g e , one finds that the partition function is given by 

Z BF (A*) = Y.Tl d *f [ II d ^ I II d * f dn rf \{ (X fl n Tf \g^g aTl \X f ,n T , f ). (4.34) 

*/ f (t,«t) J (r,f) (a,f) 

The idea implemented in the FK model is that the simplicity constraints have to be imposed 
directly in (4.34) as certain restrictions on the representations Xf and the normals n T f labeling 
the Spin(4) coherent states. 

In order to translate the simplicity conditions in terms of these data, one uses the semi- 
classical nature of the coherent states. The idea is that the state |A/, n r j) encodes the geometry 
of a semi-classical tetrahedron in the sense that Xf represents the area of the triangle dual to 
the face /, whereas n T f is the normal to the triangle as viewed from the tetrahedron r. To 
extract this information, one should evaluate the expectation value of the appropriate quantum 
operator. Since classically the tetrahedron is described by the discretized B field, it is natural 
to expect that one should consider its quantization Bf which is provided by the map (4.27). As 
a result, we have to consider the following semi-classical bivectors: 



XfcU/) = <A/,Hr/| ( TlJ - l -* TlJ ) |A/,n T/ ). 



The simplicity constraints are then imposed by requiring that this bivector is simple, i.e. it can 
be represented as (the Hodge dual of) wedge product of two one-forms, and therefore it should 
satisfy (*X( A/ . nT/ )) /J (x r )j = for all / C r and some normal vector x T . Using (4.33), this 
condition can be easily solved and leads to a restriction on Xf, 
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and to a restriction on n T f, 

/ |i-i| \ 
M 55 ^/^, "WV/ e J' with u xt=3x t (9x t ) i ( 4 -36) 

where ^ is as usual a representative of s £ X in Spin(4), the U(l) elements Ha, = e*^ " 3 take care 

of the fact that the are defined up to a phase, and e is either 1 or the matrix *J such that 

/ie = e/i for all /i € SU(2), depending on whether 7 is smaller or larger than 1. These restrictions 
are then substituted in the partition function (4.34) so that instead of a four-dimensional integral 
over n T f one is left with a two-dimensional integral over n T f only. 

Notice that the first condition (4.35) agrees with the restriction on representations (4.29) 
appearing in the EPRL model and can be considered as a consequence of (4.28a). At the same 
time, the second condition (4.36) clearly distinguishes the two cases 7 < 1 and 7 > 1. It turns 
out that the model behaves very differently in these two ranges of the Immirzi parameter [75]. 
In the first case, it literally coincides with the EPRL model. In particular, the boundary states 
and the vertex amplitude are given by (4.30) and (4.31), respectively. In the case 7 > 1, the 
presence of the matrix e in (4.36) has drastic consequences. The boundary states are spanned 
by all projected spin networks with representations A/ satisfying the constraint (4.35). In other 
words, the simplicity constraints do not lead to any restrictions on the intertwiners. What they 
fix is the weight of each intertwiner in the partition function. More precisely, the partition 
function of the FK model can be represented as the following sum: 

13} 

^k 7>1 =ee e n^/^n^wnKctin^cj/.jr/.ir), 

jf ir jrf=jf f T Tf V 

where Af and A e are the face and edge amplitudes given by very simple factors in terms of 
dimensions of the corresponding representations, A v is the vertex amplitude given by the formula 
similar to (4.31) but with the arbitrary representations j T f one projects onto, and finally 

r j = ( r jjk \ 2 = m m 

\yj-30) (2j-k)\(2j + k + l)\ 

is the factor fixing the weight of the projected spin network intertwiners. 

Since the FK model for 7 < 1 is identical to the EPRL model, it does not require a separate 
discussion. In contrast, the FK construction for 7 > 1 suggests an alternative spin foam model, 
which however has not been investigated as deeply as its EPRL cousin. In particular, its relation 
to the canonical quantization remains unclear. On the other hand, there are several evidences 
in favor of this model compared to EPRL [8] one. In particular, in the limit 7 — > 00, the EPRL 
model reduces to the original BC model, whereas the FK model is truly different. But as is well 
known, the BC model suffers from various problems. Furthermore, in [75] it was shown that 
this model can also be reproduced in the FK framework, but involves an erroneous step, which 
makes its results more than suspicious. Due to these reasons, for an Immirzi parameter 7 > 1, 
the FK model might be considered as most promising one. 

4.4 EPRL vs. LQG 

4.4.1 Identification of kinematical states 

Let us take a closer look at the results of the EPRL model. One observes that in both Euclidean 
and Lorentzian signature, the boundary states, represented by the particular class of projected 
spin networks (4.30), are labeled by graphs colored by one spin per link and one invariant SU(2) 
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intertwiner per node. But this is precisely the coloring carried by the usual SU(2) spin networks! 
This fact suggests to identify the space of boundary states with the kinematical Hilbert space 
of LQG [66]. 

This identification can be mathematically formalized as a certain bijection map r from %kin 
to the space K, spanned by the states (4.30). It is given by the following integral formula [63]: 

T- St T3 fi i-> S^ P J^ L) [g e ,x n } = / Y[[&h t K(g~l ■ g e ■ g h e )) S {rJ ^ [h e ], (4.37) 

K ' J ' ' JSU(21 „ V * > 



where the kernel is 



" o U ( 2, ) 



'SU(2) 

and Xg(s) = TrD( A )(g) denotes the character of representation A. This map lifts SU(2)-invariant 
states to Lorentz-invariant ones. However, these Lorentz-invariant states are not arbitrary. The 
image of t can be characterized by the property that the reduction of t[S] to the SU(2) subgroup 
coincides with the original SU(2) state S. In other words, the inverse map is given the following 
simple restriction 

r- 1 : 5™ h+ 5 (r ^N =5™[^,x (0) ], (4.38) 

where avg) = ^ s ^ ne norma l vector in the time gauge. 

This construction however raises some issues. First of all, in the Euclidean case the map is 
defined on the whole %kin only for very special values of the Immirzi parameter. Indeed, let 

1 + 7 V 



|l-7l 



where p, q are two coprime integers. Then (4.29) implies that A S [N/2] 2 and therefore the 
map (4.37) is well-defined only for [16] 



1 , 

-m{p + q), 7 < 1, 
1 , 

-m{p-q), 7 > 1, 



These values run over all half-integers only for 7 = 2n + 1, n G N. From the canonical point of 
view, it is not clear at all why the Immirzi parameter should be restricted to one of these values. 
Thus, in Euclidean signature the identification seems to be problematic. 

Fortunately, in the physically interesting case of Lorentzian signature, there are no restrictions 
on the Immirzi parameter and the map (4.37) is always well-defined. But here one encounters 
another problem: the norm of the states on K, induced from the scalar product on the space of all 
projected spin networks is divergent. This divergence arises due to the delta function 5{p — p') of 
continuous representation labels and can be traced back to the non-compactness of the Lorentz 
group. It has been suggested in [138] that the relevant scalar product on the constrained state 
space K, is obtained by replacing 8\\> = S nn /5(p — p') by Sjf since A is determined by the discrete 
variable j. After such a replacement, r becomes a unitary map. One immediately recognizes that 
this change of the scalar product is exactly the same as the passage from (4.7) to (4.12) in the 
canonical approach and therefore it is properly justified. However, it is important to note that 
it is equivalent to a certain change in the integration measure over holonomies which acquires 
a 5-function of the secondary second class constraints. As we will see in the next subsection, 
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if this modification is done not only in the scalar product, but also in all other places in a 
consistent way, it gives rise to interesting new results for the spin foam quantization. 

Thus, it seems that, at least in the Lorentzian signature, the EPRL model reproduces the 
kinematical structure of LQG. In fact, not only the states are identified, but there is also a way 
of defining quantum operators of area and volume in the spin foam context such that their 
spectra coincide with those of the canonical quantization [62] . 

However, at the dynamical level the situation is far from being clear. What has been exten- 
sively studied is the asymptotics of the EPRL vertex amplitude at large spins [34, 35, 55, 117]. It 
has been shown that it is dominated by non-degenerate simplicial geometries whose contribution 
scales as the exponential of the Regge action. Although this fact is promising, it is not sufficient 
to claim that the model has the correct quasiclassical limit. There were also some attempts to 
relate this vertex amplitude to a Hamiltonian operator in LQG [5, 43], but this issue remains 
largely unexplored. 



4.4.2 Problems 

It is natural to compare the map r (4.37) to the construction of the covariant representation of 
LQG of Section 4.1.3. More precisely, the inverse map (4.38) seems to be essentially the same 
as the equality (4.11) reducing the projected spin networks to the usual spin network states. 
However, between these two there is a crucial difference: the map r is a formal identification, 
whereas (4.11) is a consequence of special properties of the connection on which these states are 
evaluated. This observation raises the question of why the holonomies entering in the definition 
of the boundary states are restricted to belong to an SU(2) subgroup. In the case of covariant 
LQG, this follows from the property (4.10). What is a similar reason in the EPRL model? 

First of all, it is impossible that this is the same reason, namely that the holonomies are 
defined with respect to a connection satisfying the same constraints. The EPRL model is derived 
such that its boundary states are functionals of the usual spin connection oj whose holonomies do 
not satisfy any restrictions and span the whole Lorentz group. Another possibility would be that 
it is the projection entering the definition of projected spin networks that is responsible for the 
restriction to the subgroup. Indeed, the spin connection satisfies a very appealing property [8]: 

n U) . W W .^U) =AV\ (4.39) 

where A a is the su(2)-valued Ashtekar-Barbero connection (4.1), the upper index in parenthesis 
indicates an irreducible representation, -k^ denotes the projector to Hq^m) m the decomposi- 
tion (4.3), and A is supposed to be related to j (in the Lorentzian case) as X(j) = (j,~f(j + 1)). 
This is a slight modification of the condition (4.29). It also provides a solution to the con- 
straints (4.28). Moreover it is an exact solution for the second equation, i.e. it does not require 
any additional large-spin approximation or ordering fitting in (4.28b) [62]. However, the pro- 
perty (4.39) is not sufficient for our purpose since we need its exponentiated version 

n U) . dW(j,( u )) • tt u) = B^{h £ (A)). (4.40) 

In other words, the EPRL model projects at nodes only, while LQG requires to project contin- 
uously at all points of the graph. It is clear that these two projections do nto coincide and the 
equality (4.40) is generally not true. In the absence of a justification for the reduction to the 
SU(2) subgroup in the spin foam context, it is impossible to claim that at the kinematical level 
the complete agreement with the canonical quantization has been already achieved 27 . 



27 Note also the following difference: the covariant form of LQG can be formulated using the projected spin 
networks with arbitrary A^, whereas the EPRL model requires a specific map X(J). See p. 71 for further discussions 
on this issue and on how it can be resolved in the spin foam context. 
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Furthermore, going beyond the kinematical level, one reveals another, more serious problem of 
the EPRL model. To see this, it is again indispensable to make a comparison with the canonical 
approach. Namely, let us compare the imposition of the simplicity constraints in the EPRL 
model and in the canonical quantization. The main difference between the two quantization 
procedures is the one which is at the heart of the spin foam approach: whereas the canonical 
quantization implies that the second class constraints should be imposed at the classical level, 
either by solving them explicitly or by introducing the Dirac bracket, in spin foam models the 
constraints are implemented at the quantum level via certain conditions on an auxiliary state 
space. As a consequence, in the two approaches one quantizes completely different symplectic 
structures: in the canonical quantization it is the symplectic structure induced on the constraint 
surface, and in the spin foam approach it is the symplectic structure of BF theory. As a result, 
the simplicity constraints are realized in completely different ways: either as operators identically 
vanishing on all states, or as some non-trivial operators annihilating only a certain "physical" 
subspace. 

As has been shown in [8], the second quantization procedure, employed in most spin foam 
models, leads to internal inconsistencies. For example, the Hamiltonian turns out to be not well- 
defined on the "physical" subspace given by the constraints and its eigenstates cannot belong 
to this subspace. This last fact, in particular, implies that the kernel of the Hamiltonian, which 
provides the real physical states of the theory, has a vanishing intersection with the boundary 
state space determined in this way. 

Moreover, the spin foam quantization completely ignores the presence of the secondary con- 
straints which are responsible for the fact that the simplicity constraints are of second class. 
Although they do follow from the primary simplicity constraints imposed at all times, this holds 
only in the quasiclassical approximation. At the full quantum level, they are allowed to have 
non-vanishing fluctuations. At the same time, in the canonical approach they are completely 
frozen, which is ensured by inserting their delta function in the path integral measure. In par- 
ticular, this is the reason why it appears in the scalar product (4.12) and, as will be shown in 
the next subsection, this insertion might have important effects on the quantum dynamics. 

Thus, some of the quantization steps leading to the EPRL model or to the natural iden- 
tifications with the LQG structures, when critically examined, show a deep difference with the 
canonical quantization, which can be traced back to the standard spin foam strategy. This strat- 
egy can be summarized by the recipe "first quantize, then constrain" meaning that one first 
quantizes the unconstrained theory, while the constraints are incorporated only at the quantum 
level. Given the above observations, the results of the EPRL model, and any other model derived 
using this strategy, should be taken with great care. 

4.5 Lessons from the canonical quantization 

Let us now try to understand what a spin foam model would look like if it was derived using the 
standard rules of Dirac's quantization. In other words, we would like to address the following 
question: which properties can we expect to be satisfied by spin foam amplitudes and their 
boundary states on the basis of the canonical analysis? 

4.5.1 Measure 

We start our discussion with the problem of the measure in the (discretized) path integral. 
Deriving spin foam models, one usually chooses a trivial measure which, after discretization, is 
given either by the Haar measure for variables taking values in a Lie group, or the Lebesgue 
measure for algebra-valued fields. On the other hand, the canonical path integral is known to 
have a non-trivial measure which might strongly affect the results of quantization. 



Spin Foams and Canonical Quantization 



61 



Let us recall the general form of the canonical path integral representing the partition function 
of a dynamical system with coordinates q, momenta p, primary and secondary second class 
constraints, <j) and ip, respectively, and a Hamiltonian H. It is given by 

Z = fvqVp V|detA| 8{4>)5{ijj) exp U f dt (p a q a - H)Y (4.41) 

where A is the Dirac matrix of second class constraints. We assume that there are no first class 
constraints since they are not essential for our discussion. If they are present, the Hamiltonian 
should be replaced by i7 e ff = Hn — {@,^}+ where SI is the BRST charge, Hq is an exten- 
sion of H commuting with f2, and O is a gauge-fixing fermion. Besides, one should add the 
Lagrange multipliers and ghosts to the kinetic term and the integration measure. If the gauge- 
fixing fermion is appropriately chosen (see, for example, [18] for the case of the Hilbert-Palatini 
formulation), the ghosts can be integrated out and one obtains the usual Faddeev-Popov path 
integral. 

One can distinguish two types of contributions to the measure in the path integral (4.41): 
determinants of Faddeev-Popov type and 5- functions of constraints. We consider the implica- 
tions of these contributions separately. The reason for this is that there is a trick [95] which 
allows to rewrite the partition function in the Lagrangian form 

Z = f VqVpV\fi(q,p)exp(i J dt (p a q a - H - A0)Y (4.42) 

where A is the Lagrange multiplier for the primary constraint, n(q,p) is a regular local measure, 
and the expression in the exponential is the Lagrangian in its first order formulation. In passing 
from (4.41) to (4.42), the ^-function of the secondary constraint has been traded for an additional 
local regular contribution which together with det A| gives the measure p(q,p)- Thus, this 
allows to avoid the contributions of the second type altogether. 

The local measure n(q,p) has been computed for the Plebanski formulation in [49, 65] and 
for the Hoist action in [65]. The result is given by 28 

m = y V, /iHoist = V 3 V, (4.43) 

where V and V are three- and four-dimensional volumes. Note that, due to the presence of the 
spatial volume, the measure is not covariant. This is in fact a familiar phenomenon which can 
be traced back to the fact that the diffeomorphisms coincide with the gauge transformations 
generated by the constraints of the canonical formulation only on-shell. The full path integral is 
by construction invariant under the gauge symmetries and in [94] it has been verified that the 
measure factors (4.43) are indeed consistent with this invariance. 

The implications of such additional factors on spin foam models have been considered in [41]. 
There are two possible approaches to take into account their contributions. The first one is to 
insert them in the initial path integral, as required by the canonical approach, and to repeat all 
the steps leading to the spin foam representation. However, since these are metric dependent 
contributions, their insertion would spoil the very first step of the spin foam quantization - the 
integration over the B field. An attempt to overcome this problem and to get a spin foam model 
along these lines has been made in [29], but the resulting discrete model looks very complicated 
and no concrete results have been obtained in this direction. 

An alternative method is to use the interpretation of areas and volumes in terms of repre- 
sentations and intertwiners, and to replace the original function /i in the path integral measure 
by a function p(A,I), local in the discretized sense, inserted already in the spin foam sum over 



28 Besides this local factor, the measure of the gravity path integral also includes gauge-fixing conditions and 
the usual Faddeev-Popov determinant. 
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representations. This approach was studied in [41] and it was shown that such contributions 
would affect the face and edge amplitudes. On the basis of this observation, it was argued that 
the proper choice of these amplitudes, which must be consistent with the measure favored by 
the canonical analysis, is crucial for the correct implementation of all the symmetries. A wrong 
choice will unavoidably result in an anomalous quantization. On the other hand, the spin foam 
approach does not suggest a unique way to fix the face and edge amplitudes, and therefore the 
link with the canonical quantization seems to be the necessary additional input which might 
allow to fill the gap. 

Let us now return to the 5-function contributions to the measure due to the secondary 
constraints. Although the trick of [95] allows to remove them in the partition function, it does 
not work for correlation functions [9]. This implies that we cannot disregard them and they 
should be included into the measure also in the discretized path integral used in the spin foam 
quantization. It turns out to be possible to derive several important conclusions just from this 
simple fact. 

From the canonical analysis of the Plebanski formulation presented in Section 4.2, we know 
that the secondary second class constraints impose restrictions on the connection variable. Thus, 
their discretized version should give rise to certain conditions on the holonomies. We assume 
that they can be described by a function V'discr (9aT,x T ) where group elements g aT , as in (4.34), 
describe holonomies from the center of a four-simplex a to the center of one of its tetrahedra r, 
and we also included the dependence on the normal x T to the tetrahedron. An example of such 
a discretized constraint is given by the one used in (4.11), which in the Riemannian case can be 
formulated as 

fecr (g<TT,x T ) = {g^ T )~ 1 g~ T u~l = I, (4.44) 

where we have used the decomposition g = (g + ,g~) for Spin(4) group elements and u x was 
defined in (4.36) in terms of g x . 

Relying on these observations, we conclude that the secondary second class constraints should 
be incorporated into the discretized path integral by modifying the integration measure over 
holonomies. We suggest that it takes the following form 

V {Xt) [g aT } ~ 5(V>discr (Sot, X T )) &9<TT (4.45) 

and satisfies the following covariance property 

V^[g aT g}=V^^[ g<7T ], g£ G, (4.46) 

which is necessary in order to ensure the Lorentz invariance of the path integral. 

To understand the implications of this measure, let us consider the discretized path integral 
for a single four-simplex with fixed Bf on the boundary. We then drop the unnecessary index 
a and denote u(f) and d(f) two tetrahedra sharing the f th triangle (as the triangle is oriented, 
one of the tetrahedra is considered as "up" and the other as "down" ) . Since there are no internal 
faces, the only integration involved in the definition of the corresponding quantum amplitude is 
the one over holonomies g T associated with each of the five tetrahedra. The integration is done 
with the measure (4.45), so the amplitude becomes 

A[B f ] = I J] 27^ [g T ] J] exp (iTr [B f g-^g d{f) ]) , (4.47) 

T f 

where in the exponential one can recognize the discretize BF action. The amplitude (4.47) can be 
viewed as a simplex boundary state in the Bf or the so-called flux representation [28, 30]. In the 
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total partition function, these contributions should be multiplied by using the non-commutative 
star-product 29 and integrated over the bivectors: 

Hv[B f ]( + A a [B f ]), (4.48) 
/ 

where the order in the non-commutative product is dictated by the structure of the simplicial 
complex A. The measure T>[Bf], similarly to the measure for holonomies, should contain the 
5- function of the discretized simplicity constraints <pf as well as various volume factors of the 
type (4.43) which can be usually constructed from the bivectors. Note that this representation 
allows to disentangle the issues related to this -By-dependent factors and to the secondary se- 
cond class constraints. It shows that the former become relevant only when one glues different 
simplices together, whereas the latter will affect the contribution of each individual simplex. 
To see how it works, we rewrite ^4[-E>/] as the Fourier transform of another amplitude, 

A W = /il [dgfe^t^A^Xr], (4.49) 
/ 

with 

A\g f ,x T ] = [ n^^^m^v/)^ 1 /))' ( 4 - 5 °) 

T f 

which thus appears as the simplex boundary state in the "connection" representation since it 
depends on the group elements <?/. 30 The representation (4.50) seems to be very similar to 
the one for the partition function of BF theory (4.18) where the S- function implements the 
flatness condition. However, as we will see, the non-trivial measure makes this similarity very 
superficial and the flatness condition does not hold in our case. In [9] it was shown that, using 
the covariance property (4.46) and inserting integrals over the SU(2) subgroup attached to each 
tetrahedron, the amplitude (4.50) can be represented as a linear combination of projected spin 
networks (4.6), 

A[g f ,x T ]= \]l d ^) \]l d Jrf\ M^fJTf,ir)S ira ^ rsSt ^[g f ,x T ], (4.51) 

where T a is the graph dual to the boundary of the four-simplex. The coefficients A(Xf,j T f, i T ), 
given below in (4.54), are naturally interpreted as vertex amplitude of the spin foam quantization. 
The result (4.51) provides a general proof of the observation made in the previous subsections 
that in all spin foam models of four-dimensional gravity, the boundary states coincide with a cer- 
tain class of projected spin networks. Furthermore, below we explain the important implications 
of this result on the implementation of the simplicity constraints in spin foam models and on 
the dynamics of the corresponding quantum theory. 



4.5.2 Quantization of constraints 

All four-dimensional spin foam models are characterized by the presence of simplicity constraints 
which should be implemented in one way or another. In principle, this problem is one of 

29 This product is defined on the plane waves e 9 (a) = e lTr( - 9 ' ' , with a £ g and momentum given by the group 
element g, as e si *e 92 = e 91S2 [28]. 

30 Note that in general the secondary constraints depend on the B field. As a result, the measure (4.45) may 
depend on the bivectors Bj as well. Here we ignore this possible complication. 
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the central topics of the canonical quantization, and therefore we should not expect a drastic 
difference in the way the constraints are implemented in the loop and spin foam approaches. 
Moreover, the Hilbert spaces before and after the imposition of constraints are identical in the 
two approaches. Indeed, on the one hand, the initial state space of the topological BF theory, 
once the closure constraint is relaxed as in (4.5), coincides with the space of projected spin 
networks. On the other hand, the space of the boundary states in a spin foam model of gravity, 
resulting from the imposition of the constraints, must be identical to the kinematical Hilbert 
space of the loop quantization. This shows that it is natural to expect that the constraints should 
be realized in spin foam models in the same way as the one responsible for the reduction from all 
projected spin networks to the kinematical Hilbert space of LQG in its covariant formulation, 
as was presented in Section 4.1.3. 

Let us recall that this reduction is done due to the fact that the holonomy arguments gg of the 
projected spin networks do not span the whole Lorentz group, which in turn can be traced back 
to the constraints satisfied by the connection defining these holonomies. Thus, these are the 
secondary second class constraints that are responsible for the reduction of the Hilbert space in 
the canonical approach, whereas the primary simplicity constraints are automatically satisfied 
due to Dirac's quantization rules. This is in drastic contrast with the usual spin foam strategy 
where the primary simplicity constraints are used to reduce the state space of BF theory, whereas 
the secondary constraints are completely ignored. Of course, the origin of this discrepancy comes 
from the fact that the two approaches quantize different symplectic structures, as was discussed 
in Section 4.4.2: one is given by Poisson brackets and describes the symplectic structure of BF 
theory, and the other is given by Dirac brackets ensuring that we deal with general relativity. 

Remarkably, the inclusion of the secondary constraints into the path integral measure leading 
to (4.51) implies that the constraints are imposed in complete agreement with the canonical ap- 
proach. To understand why this is so, let us return to the initial formula for the simplex boundary 
state in the connection representation (4.50). The <5-function expresses the arguments gj of the 
boundary state wave function in terms of the holonomies g T . But the latter are integrated out 
with the measure (4.45), which is supposed to contain an additional 5- function. This 5- function 
induces certain restrictions on g* , similar to those which hold for g T . For instance, in the simplest 
case (4.44), the constraint induced on gj reads 

9 .f = u XuU) 9 f Ux d(f) • 

Being used in (4.51), these restrictions precisely correspond to the imposition of the secon- 
dary second class constraints on the enlarged Hilbert space of the canonical loop quantiza- 
tion [11, 13]. In particular, these constraints should be reflected in the form of the vertex 
amplitude A(Xf,j T f,i T ), which might, for example, vanish on a certain set of representation 
labels. We consider this issue in more detail in the next subsection. 

Before that, we however discuss another controversial issue of the spin foam approach, which 
is the quantization of the closure constraint (4.20). Usually, it is imposed at the quantum level 
by requiring the invariance of intertwiners associated to the nodes of the boundary states: 

J^T/X T = 0. (4.52) 

fCr 

On the other hand, the canonical quantization leads to intertwiners which carry an additional 
dependence on the normals x T and satisfy instead the relaxed closure constraint (4.5). This 
constraint can be formulated in the infinitesimal form, analogous to (4.52), as [9, 11] 



^ T } Z T (x r ) = f -X T (xt), 

fCr 



(4.53) 
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where f acts on functions of x 1 6 X as a generator of G: 

Tij ■ f(x) = r] IK x K djf - r] JK x K d I f. 

Although (4.53) can be reduced to (4.52) by integrating the intertwiner over x T , the canonical 
quantization clearly shows that this should not be done. 

A moment of reflection shows that the path integral quantization leads to the exact same 
conclusion. Indeed, note that the bivectors Bf appearing as arguments in the quantum am- 
plitude (4.47) should satisfy the linear simplicity constraint (4.26) with a fixed normal x T . As 
a result, this normal should also be kept fixed in the path integral representation for -Afify] 
and, as a consequence, in all the boundary states. Moreover, it should not be integrated even 
when two simplex contributions are glued together. Assume that two simplices u\ and a 2 share 
a tetrahedron. Then the gluing (4.48) can be rewritten in the connection representation, and 
the total boundary state can be obtained from two individual simplex amplitudes as [9] 

where / = (/1, f'2, fn), t = (^l, ?2, T12), the labels 1, 2 refer to non-shared faces and tetrahedra 
of the corresponding glued simplices, and the label 12 marks the shared faces and tetrahedron. 
Using (4.46), it is easy to see that the dependence on x Tl2 drops out. Therefore, an integral 
over the normals would be completely artificial and moreover would lead to a divergence in the 
Lorentzian case. In fact, keeping x T fixed is nothing else but the usual gauge-fixing of the boost 
gauge freedom in the path integral. 

Thus, both the path integral and the canonical approaches require that the closure constraint 
be imposed only in its relaxed form (4.53). This is not in contradiction with the classical closure 
constraint, since the latter arises only as a classical equation of motion and does not need to be 
preserved by quantum fluctuations. The necessity to give up the invariance of intertwiners in 
favor of their covariance has been recently recognized also in [138]. 

4.5.3 Vertex amplitude 

Finally, we consider the implications of our prescription to include secondary second class con- 
straints in the path integral measure for the vertex amplitude, which is the most important 
quantity of the spin foam quantization since it encodes the dynamics of the quantum theory. 
This amplitude is given by the coefficients A(Xf,j T f,i T ) in the decomposition (4.51) and is 
found to be [9] 

A(X f ,j T f,ir) = [ H I)(lr) ^l 5 (r^„A„Mk(/)9 d (/). I 4 ( 4 - 54 ) 

J T 

It is easy to check that the covariance of the measure (4.46) ensures the independence on the 
normals x T . Thus, the vertex amplitude is given by an integral over gauge transformations at 
the nodes of the projected spin network associated with the boundary of a four-simplex. 

This formula is a straightforward generalization of the usual prescription for the evaluation 
of vertex amplitudes in the spin foam models considered in the literature. The only difference 
is that the measure in (4.54) is assumed to be non-trivial and, in particular, to include the 
secondary constraints. Taking instead the measure to be the usual Haar measure on the group, 
T>^ Xt '[g T ] = dg T , one reproduces an analogue of (4.31). Moreover, in this case the integrals over 
the gauge group are equivalent to the integrals over the normals x T . Therefore, as was discussed 
in the previous subsection, they can be viewed as a part of the definition of the boundary states 
ensuring the closure constraint (4.52). Then one gets back the simple recipe that the vertex 
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amplitude is given by the simplex boundary state evaluated on a flat connection, which can be 
traced back to the flatness condition of BF theory. 

However, if one insists on the presence of the secondary constraints in the measure 
this simple prescription is not valid anymore. In particular, flat connections do not play any 
preferred role and the dynamics appears to be more complicated. The concrete form of the 
vertex amplitude depends on a particular choice of the measure which should be dictated by 
agreement with the canonical quantization. 

4.5.4 Choice of connection 

So far our consideration was generic as we did not specify the precise form of the secondary 
constraints. These constraints depend on the choice of the connection which is used to define 
holonomy variables. Thus, we encounter the same issue which already arose in Section 4.1.3 
when we discussed the Lorentz-covariant formulation of LQG. Here we should answer the same 
set of questions: Which classical variable behaving as a connection should be chosen to define 
holonomies? Which constraints are satisfied by this variable? And how can these constraints 
be discretized? 

Usually, in the spin foam approach, no connections other than the standard spin connec- 
tion lo ij appear. Thus, one may expect that this is the connection to be considered. But then 
one runs into problems similar to those which did not allow to use it for covariant formulation of 
LQG. The second class constraints satisfied by the spin connection are given by * ab from (4.17) 
and do not allow a simple geometric interpretation in terms of holonomies constructed from or. 
Given the troubles with the loop quantization in terms of this connection, it is natural to ask 
for something else. 

The next evident candidate is the Lorentz generalization of the AB connection A IJ (4.9). 
Since it was essential in the covariant reformulation of LQG, one might expect that it provides 
a direct link between LQG and spin foam models. This connection satisfies the constraints (4.10) 
which allow a very simple discretization as in (4.44). However, if one plugs these constraints into 
the formula for the vertex amplitude (4.54), one finds as a bizarre result that it is given by the 
vertex amplitude of SU(2) BF theory. This puzzle is easily resolved by noting that the above 
results, and the vertex amplitude (4.54) in particular, are not applicable in this case because 
of two reasons. First, the connection A IJ is not a pull-back of a spacetime connection [7, 140]. 
In other words, it is not possible to extend it to a spacetime connection which transforms 
properly under all gauge transformations including time diffeomorphisms. In the absence of 
such a spacetime interpretation, it becomes impossible to use this connection in the covariant 
path integral quantization. The second reason is even more clear. The connection A IJ differs 
from the spin connection lo ij by a complicated term non-vanishing on the constraint surface. 
As a result, the Plebanski or even the BF action, which is the starting point of our analysis, 
takes a completely different form when it is written in terms of A IJ . Altogether, these issues 
make the use of A IJ in the spin foam quantization very problematic. 

There is another possible choice of the connection which avoids the problems of the AB 
connection A IJ and partially solves the problems of the spin connection lu ij . This is the so- 
called "shifted connection" introduced in [17]. Its name comes from its definition as a shift of 
the spin connection by a term proportional to the Gauss constraint, 

T u, KL 

aIJ — ,JJ< j 7 fST tdPQt 

a ~ a 2(l + 7- 2 ) ■' KL ' PQ ~ a ST ' 

Here ffj N KL are so(3, 1) structure constants, P^ J = g~ 1 g a bP b ' IJ where P a is the canonical 
field of the Plebanski formulation defined under (4.16), g a b = ViJ eI a e b * s the spatial metric 
and g is its determinant, and we defined I 7 = (l — 7 -1 *) / (l — 7- 1 *) with lj/ L the projector 
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introduced under (4.8). As it was shown in [7], A IJ represents the unique spacetime Lorentz 
connection whose holonomies diagonalize the action of the area operator. Due to this last 
property, this connection is suitable for the loop quantization. The corresponding quantization 
has been called covariant loop quantum gravity (CLQG) [15] and is inequivalent to LQG. Its 
nice feature is that it leads to an area spectrum which is independent of the Immirzi parameter. 
This result seems to be much more natural than the standard result of LQG, where 7 appears 
as a proportionality factor, since classically 7 is a free parameter removable by a canonical 
transformation. However, in contrast to the AB connection, the shifted connection is non- 
commutative. This fact leads to serious difficulties in quantizing the corresponding symplectic 
structure and its proper representation in terms of an operator algebra on the Hilbert space 
spanned by projected spin networks has not been found yet. We refer to [16] for further details 
on this approach. 

On the other hand, the fact that the shifted connection has a spacetime interpretation and 
differs from the spin connection by only a term vanishing on the surface of the Gauss constraint 
allows to use the shifted connection in the spin foam approach. But we still need to find a dis- 
cretization of the secondary second class constraints, which in terms of A IJ take the following 



where F a is the SL(2, C) connection compatible with the field P constructed from the B field of 
the Plebanski formulation. The constraints (4.55) are somewhat easier than the corresponding 
constraints ^> ab (4.17) on the spin connection, since they involve a projector on the subgroup. 
However, the dependence of the right-hand side on the B field, which is in fact responsible 
for the non-commutativity of A IJ , introduces additional complications. As a result, a proper 
discretization of (4.55) remains an open question. 

Thus, in both the canonical and spin foam approaches, the use of the shifted connection 
did not lead to a well-defined model for quantum gravity yet. However, the difficulties arising 
on this way seem to be of similar nature and the shifted connection appears to be a unique 
object which is suitable for both quantization schemes. Therefore, from the point of view of 
the agreement between the canonical and spin foam quantizations, it can still be viewed as 
a promising candidate which should be used in the definition of holonomy operators in both 
approaches. 

5 Back to three dimensions — a toy model 

From our discussion in the previous section, it should be clear that all problems and ambiguities 
arising in four- dimensional spin foam models come from the difficulties in imposing the simplicity 
constraints, which are supposed to turn the trivial dynamics of the topological state sum model 
into that of quantum general relativity. Although we have argued that agreement with the 
canonical approach could play the role of the guiding principle, it would be nice to have a simple 
model allowing to test various ways of implementing the simplicity constraints. In order to do 
so, such a model should possess the following features: 

• It should have the form of a Plebanski theory, i.e. be written as 



form: 




(P) 



(4.55) 



St\i2 = Stm + constraints. 



(5.1) 



Both theories, given by the actions Srhi an< ^ ^Th2, should have known spin foam repre- 
sentations. 
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In this case, we could verify which of the methods to impose the constraints reduces the spin 
foam quantization of Stm to that of <Srh2- Those methods which do not satisfy this requirement 
most likely give wrong quantizations also for the original Plebanski theory (4.13) and therefore 
have few chances to describe quantum gravity. 

A model of this kind has been proposed in [9] , where the role of the theories Thl and Th2 is 
played by four-dimensional BF theories with Spin(4) and SU(2) gauge groups, respectively. The 
constraints that reduce Thl to Th2 are given by the analogue of the linear simplicity (4.26), 
whose discretization reads 

BY{x T )j = 0, (5.2) 

and by the secondary constraints (4.10) relevant for the Lorentz generalization of the AB con- 
nection, which after discretization take the form (4.44). It was shown that one recovers the 
known spin foam representation of SU(2) BF theory from the known spin foam quantization of 
Spin (4) BF theory only if one incorporates the secondary constraints in the way suggested in 
Section 4.5. In particular, the vertex amplitude should be modified according to (4.54). 

However, the model of [9] has one important drawback which makes the above conclusions 
uncertain: it has never been shown whether there is an action of the form (5.1) with primary 
constraints (5.2) and secondary constraints (4.44). In other words, the constraint structure was 
imposed by hand and did not follow from the Hamiltonian analysis of any action. 

On the other hand, in [88] it was understood that one finds a well-defined model with 
all the above characteristics if one goes down to three dimensions. In other words, the authors 
proposed a three-dimensional Spin(4) Plebanski-like action in which an analogue of the simplicity 
constraints is used to recover Riemannian gravity, i.e. SU(2) BF theory in three dimensions. 
Then one can take advantage of the well-established results for the spin foam quantization of 
three-dimensional BF theories to test various proposals, such as BC, EPRL 31 and the one from 
Section 4.5, for imposing the simplicity constraints in four dimensions. 

The key observation is that the action 

S[u,B,\]= f d?x(e^PT?{B^F up ) + \^Tr{B^B v )), (5.3) 
Jm 

in which B IJ is an so(4)-valued one- form field, generates the simplicity constraints 



1 



%„ = -eu KL B^B^ = 0, (5.4) 

which behave like in the four-dimensional Plebanski theory. In particular, it can be shown that 
these constraints admit three sectors of solutions. Apart from a degenerate sector in which the 
three-volume V is vanishing, one has the following solutions: 

• topological sector: Bf/ = x 1 ' — x J ej 1 , 

• gravitational sector: Bj/ = e IJ KL x K e^, 

where x 1 E M 4 is a vector analogous to the normal introduced in Section 4. It is easy to see that 
in the topological sector the action (5.3) is vanishing, while in the gravitational sector it reduces 
to the usual SU(2) Riemannian action for three-dimensional gravity. The Hamiltonian analysis 
of this sector shows [88] that the simplicity constraints <3? a f, (5.4) generate secondary constraints 

y ab = Tr(D (a B *B b) ) = 0, (5.5) 



31 Strictly speaking, we are going to use the prescription of the EPR model [67], since we do not consider the 
inclusion of the Immirzi parameter. 
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which do not commute with the primary ones and together form a second class system. As 
their four-dimensional counterpart ^ ab given in (4.17), they can be seen as constraints on the 
components of the spacetime connection uj ij . 

Now we can quantize the theory (5.3) in two ways. First, we impose the constraints (5.4) 
classically and restrict ourselves to the gravitational sector. Its spin foam quantization is well- 
known and provided by the Ponzano-Regge model introduced in Section 3.3.1. It is defined by 
the partition function 

z PR (a*) = y, n (-!) % ( 2 i/ + 1) n ™ ( 5 - 6 ) 

j-+f /eA* »eA« 

with the vertex amplitude given by the SU(2) {6j} symbol. On the other hand, we can start 
from the unconstrained Spin(4) BF theory and try to impose the simplicity constraints at the 
quantum level in order to mimic the situation in four dimensions. It is clear that the state sum 
model for the unconstrained part is given by 

Z Spin(4) (A*) = E II (-l)^ ++2j7 (2j+ + l)(2jj + 1) Jl {6i + }{6r}. (5.7) 

(i+,i-)->//eA* v&a* 

Thus, we come to the following simple question: is there a way to incorporate the constraints 
into the partition function (5.7) which would reduce it to the Ponzano-Regge model (5.6)? 



The usual strategy 

First we would like to test the ideas used to derive the BC and EPRL spin foam models. Accor- 
ding to them, the simplicity constraints are implemented as restrictions on the representations 
and intertwiners labeling the faces and the edges of A*, respectively. These restrictions are in 
turn obtained by imposing the primary simplicity constraints (5.4) on the boundary states of the 
unconstrained theory. As in four dimensions, at the discrete level one can again distinguish the 
diagonal and cross simplicity constraints. The former is always imposed strongly and requires 
that the Spin(4) representations assigned to the faces be always simple, jt = jj = jf. On the 
other hand, the imposition of the cross simplicity depends on the chosen approach and is expected 
to fix the intertwiners labeling the edges of A*. In this way one obtains the following results [88]: 

• BC prescription: After the strong imposition of both the diagonal and cross simplicity 
quadratic constraints, the partition function becomes 

iw(A*) = £ n (d Jf ) 2 n {6j} 2 . (5.8) 

j^//eA* veA* 

• EPR prescription: One replaces the quadratic cross simplicity by its linear version (5.2) 
which is then imposed in the weak sense. The resulting partition function reads 

iw R (A*) = £ n K /2) ) 2 n (m 2 n > ^ 

with a, = 1, . . . , 4 labeling the edges of the tetrahedron dual to the vertex v and the pair 
a/3 labeling the corresponding faces. The fusion coefficients /*+•_ attached to edges are 
given by the evaluation of a spin network and produce the {9j} symbol 



Ji/2 
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Actually, they depend only on three spins jf^ e since all intertwiners appearing in their 
definition are unique. 

As a result, the EPR-like partition function (5.9) differs from the one given by the BC pre- 
scription (5.8) only by a normalization factor assigned to the edges. In fact, this factor appears 
because the intertwiners obtained by imposing the cross simplicity a la EPR have not been 
properly normalized and the operator mapping (5.7) to (5.9) is not a projector [100]. Once the 
normalization is corrected, one finds that 

^3dEPR-norm = ^3dBC- 

This fact has a simple explanation. Note that since we are in three dimensions, for any simplicial 
discretization each edge bounds three faces. Therefore, from the very beginning, up to norma- 
lization, there is a unique gauge invariant intertwiner which can be assigned to the edges of A*. 
As a result, in this case the effect of the cross simplicity is trivial and, after proper normalization, 
the BC and EPR models coincide. However, they do become different on a generic discretization. 

What is important however for our discussion is that neither the BC nor the EPR prescription 
leads to the Ponzano-Regge partition function (5.6). The difference is twofold. First of all, the 
face amplitudes do not reproduce the expected result dj . . This is in fact not surprising because, 
as was discussed in Section 4.5.1, the lower dimensional amplitudes depend crucially on the path 
integral measure, which we did not attempt to fix here properly. The second difference is more 
important: both prescriptions lead to a wrong vertex amplitude. The resulting amplitude is 
nothing else but the original vertex amplitude of Spin(4) BF theory restricted to the allowed set 
of representations, exactly as in four-dimensional spin foam models. However, the correct vertex 
in our example is the SU(2) symbol, which is clearly not equal to its square. Thus, we have 
to conclude that the standard spin foam strategy "first quantize, then constrain" leads to the 
wrong quantum dynamics, as was suspected in Section 4.4.2 on the basis of the canonical analysis. 

Improved quantization 

Since the usual methods to incorporate the simplicity constraints into the spin foam state 
sum failed, we need to try something else. Our basic observation is that the model (5.3), 
just like the four-dimensional Plebanski formulation, possesses the secondary second class con- 
straints fy a b (5.5), which impose certain restrictions on the components of the connection uj. 
As has been argued in Section 4.5, the presence of such constraints requires a modification of 
the integration measure for holonomies. Here we will show that the results presented in that 
section allow to recover all the ingredients of the Ponzano-Regge model and, combining them 
together, one finds precisely the state sum (5.6). It should be noted that the derivation we are 
going to give here works also in the same way for the four-dimensional model proposed in [9] 
with the only difference that the three-dimensional simplicial complex should be replaced by its 
four-dimensional analogue. 

Let us fix the normal field x 1 , which allows to rewrite the simplicity constraints in the 
gravitational sector in the linear form 



This normal breaks the Spin(4) gauge group to an SU(2) subgroup stabilizing x 1 . As a result, 
the boost components of the Gauss constraint become second class and together with (5.5) can 
be written in the same way as the constraints for the Lorentz extension of the AB connec- 
tion (4.10), i.e. 



b; j xj = 0. 




(5.10) 
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In particular, in the time gauge x 1 = 5q, the meaning of (5.10) is simply that the left and 
right chiral parts of the spin connection are the same, oj^ = uj. These constraints are easy 
to discretize and their discrete version is given by (4.44). It is clear that, when inserted in the 
measure, they simply reduce the integration to the SU(2) subgroup. 

Now we are ready to use the prescription of Section 4.5 [9]. First, the vertex amplitude, 
which is given by (4.54) now becomes 32 



= {6i}nV^(/)/M/)/' (5 - n) 
/ 

where the {6j} symbol is constructed out of six representations jf = j u (f)f = 3d(f)f an d we have 
used the property 33 



m,m' ,n,n 



Thus, up to the normalization factor, we have reproduced the vertex amplitude of the Ponzano- 
Regge model! 

However, we should still understand the following issue. The prescription that we are using 
here is based on the use of an enlarged state space spanned by all the projected spin networks. 
This is reflected, for instance, in the formula (4.51) for the simplex boundary state in the 
connection representation where it is represented as a sum over all these spin networks. On the 
other hand, the result (5.11) implies that the label Xf is auxiliary since the vertex amplitude 
does not depend on it. How is this degeneracy removed? Or in other words, how do we reduce 
to the boundary state space of the Ponzano-Regge model which is spanned by the usual SU(2) 
spin networks? 

One could expect that this is done by imposing the primary simplicity constraints which 
relate the labels A/ to jf, as it happens for example in the EPRL model. However, this would 
still leave us with the problem, already discussed in Section 4.4.2, that such restrictions on the 
labels are not enough to reduce the projected spin networks to the usual ones. In fact, as we 
will see now, the mechanism of the reduction is different and is analogous to the one in the 
canonical quantization. 

The idea is to consider the simplex boundary state (4.51) where the projected spin networks 
are weighted by the vertex amplitude. Since the latter does not depend on Xf, the sum over 
this label can be done explicitly. Indeed, it is easy to show the following identity: 



j'+j— m,m',n,n' 

due to which the simplex amplitude can be evaluated to 



A[g f ,x t ] = Y j m 

if 



Ud^SigJixXg+ix))' 1 ] 
f 



S^^gjix)], (5.12) 



where we have introduced <?/(x) = g x ^ (f) gfgx d{f) and 5(r T j) is the usual SU(2) spin network 
associated with the boundary graph of a tetrahedron r. This result shows that the sum over the 



32 Recall that we are in three dimensions so that the edges of A* are dual to triangles of A. This is why the 
index r was replaced by t. 



33 The normalization of the Clebsch-Gordan coefficients C^f^lj^ is chosen as in [75]. 
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auxiliary representation labels restores the secondary second class constraints on holonomies, 
which in turn provide the reduction from the enlarged state space to the kinematical Hilbert 
space of the Ponzano-Regge model. It is worth to note that the primary simplicity constraints 
have not been used so far. It is the secondary constraints that allowed to reproduce both the 
kinematical Hilbert space and the vertex amplitude implementing the quantum dynamics. 

The last step is to glue the amplitudes (5.12) associated to different tetrahedra together to 
get the full partition function. This can be done by using (4.48) with (4.49), and this is the 
only place where the primary simplicity constraints appear explicitly. Due to these constraints 
and the (5-function in (5.12), the integral over the bivectors Bf generates the 5-function on the 
SU(2) subgroup only. It imposes the flatness condition on the product of the holonomies h T f = 
g^f( x ) = 9rf( x ) aroun d each face. Expanding the 5- function in the sum over representations, it 
is easy to see that the full partition function is given by 

U[ms<r T M[hrf]]- 

T 

It is immediate to check that, doing the remaining integration and contracting all indices, one 
reproduces the Ponzano-Regge state sum (5.6) with the same face, edge and vertex amplitudes. 

This derivation explicitly demonstrates how the second class constraints should be incorpo- 
rated into the discretized path integral and how they allow to recover the right quantization for 
the constrained theory. In fact, it is not really surprising that our procedure led to the correct 
partition function because the formula (4.48) we have started with, after integrating out the 
off-diagonal degrees of freedom fixed by the second class constraint, is trivially equivalent to 
the initial expression (3.33) for the Ponzano-Regge partition function. What we have done is 
just comparing the spin foam representations of (4.48) and (3.33) obtained after certain ma- 
nipulations, which of course cannot spoil the initial equality. In more complicated cases, it is 
impossible to solve the second class constraints explicitly and the strategy based on the partition 
function (4.48) and the vertex amplitude (4.54) is the only one which remains at our disposal. 

6 Conclusions 

In this work, we have reviewed the loop and spin foam approaches to quantum gravity in 
three and four spacetime dimensions, focusing on the issue of their mutual consistency. Let us 
summarize the main lessons to be learnt from this study. 

In the three-dimensional Riemannian case, we have a beautiful consistency between the canon- 
ical quantization and the covariant path integral represented by the Ponzano-Regge (for A = 0) 
and the Turaev-Viro (for A > 0) spin foam models. Moreover, both models and their canonical 
counterparts are easily extended to include point-like particles. Nevertheless, even this nice 
picture contains important gaps. The most serious one is that, whereas the canonical quanti- 
zation dual to the Ponzano-Regge model is obtained by using the standard loop quantization 
technique, the Turaev-Viro model has been reproduced only from the so-called combinatorial 
quantization presented in Section 3.2.2. In this case the loop approach has not been realized yet 
to a full extent and its agreement with other approaches, which have proven already to be quite 
successful, remains an open issue. In particular, the crucial question is whether it gives rise to 
the quantum group structure provided in this case by U 9 (so(4))? We have presented some steps 
towards this result in Section 3.3.3. 

It is worth stressing that three-dimensional gravity in the Lorentzian signature represents so 
far almost an unexplored territory. The physical states are in general unknown, and no spin 
foam quantization (at least mathematically well-defined) has been found so far. Only the case 
where A is positive has been studied in great details because it is dual to three-dimensional 
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Riemannian gravity with a negative A. On the other hand, gravity with A < is much more 
interesting because this is the only case where black hole solutions exist in three dimensions. 
One should probably pay more attention to these cases since they are still much simpler than 
the physically interesting Lorentzian four-dimensional case. They can provide new important 
lessons and generate new interesting ideas. 

In four dimensions, the situation is not so clear. On the one hand, the spin foam models that 
have recently been introduced have the advantage of solving some long standing problems in 
the spin foam quantization, and they provide structures reminiscent of those in loop quantum 
gravity. For example, the construction of the EPRL model seems to give rise to the reduction 
from the full gauge symmetry group to its SU(2) subgroup, so that the boundary states of 
this model are claimed to coincide with the kinematical states of LQG. On the other hand, 
a detailed comparison with the canonical approach reveals crucial differences and furthermore 
raises some questions about the validity of the whole quantization procedure used to get to these 
models. Thus, here one cannot claim that a complete agreement has been achieved, even at the 
kinematical level. 

Why is the four-dimensional case so different from the three-dimensional one? The main 
difference between them can be traced back to the existence, in four dimensions, of certain second 
class constraints called simplicity. This becomes especially clear in the Plebanski formulation 
presented in Section 4.2. These constraints are responsible for breaking the topological invariance 
and generating the non-trivial dynamics of general relativity. At the same time, they also 
represent the deviation point between the two quantization approaches: while in the canonical 
approach the simplicity constraints are incorporated at the classical level, either by means of the 
Dirac bracket or by solving them explicitly, in the spin foam approach they are supposed to be 
imposed only at the quantum level. This difference in the quantization strategy is at the origin 
of all the other differences in the results and final constructions of LQG and spin foam models. 

Moreover, as we briefly reviewed in Section 4.4.2, it is this unusual strategy to deal with 
second class constraints used in the spin foam approach that is the reason for various inconsis- 
tencies expected to appear at the dynamical level. In fact, this strategy differs from the usual 
one implied by Dirac's quantization in two aspects. First, as we already said, the simplicity con- 
straints are imposed after quantization, so that one quantizes really a non-physical symplectic 
structure. Second, it completely ignores the presence of the secondary constraints conjugated 
to the primary simplicity. A nice demonstration that such a strategy is indeed erroneous is 
provided by the model introduced in [88] and presented in Section 5. 

But is there an alternative way to get spin foam models in four dimensions and how should 
a model consistent with the canonical quantization look like? We have presented some tentative 
answers to these questions in Section 4.5. In our opinion, it is crucial to start with the path 
integral with the right measure, which must agree with the measure implied by the canonical 
quantization. In particular, besides various volume-dependent factors indispensable for the 
correct implementation of the diffeomorphism symmetry, the measure should contain <5-functions 
of the secondary second class constraints mentioned above. It turns out that the inclusion of 
these constraints has far reaching consequences and strongly affects the dynamics of the theory. 

First of all, the secondary second class constraints ensure the reduction from an enlarged 
Hilbert space consisting of all projected spin networks to the kinematical Hilbert space in exactly 
the same way as the Lorentz covariant formulation of LQG reduces to the usual formulation 
based on the SU(2) gauge group, as recalled in Section 4.1.3. Thus, it becomes possible to 
achieve agreement between the canonical and spin foam quantizations at least at the kinematical 
level. But most importantly, the inclusion of the secondary constraints changes the standard 
prescription for the vertex amplitude: it is still given by an integral of the simplex boundary 
state, but now with a modified measure (see (4.54)). In fact, it is quite natural to expect such 
a modification, since it introduces a new input into the quantity determining the dynamics at 
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the quantum level. Without it, the dynamics turns out to be completely determined by the 
kinematics via kinematical boundary states. In contrast, now it becomes dependent on the 
particular form of the Hamiltonian, since the secondary constraints do depend on its form. 

Again all these conclusions find a nice confirmation in the model of Section 5. Of course, 
whether they can be realized in the full four-dimensional theory of quantum gravity remains to 
be seen. Nevertheless, we hope that the ideas and techniques presented in this review will be 
helpful to go along this direction. 
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